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Summary. [Bounds have been placed upon the 10-90°% rise time and overshoot, 
of the unit step response of phase distortion free networks. Two bounds on the rise time 
ire given, one, in terms of the area under the amplitude function curve, and the other 
in terms of the area under the square of this curve. Two bounds on the overshoot are 
given in terms of these areas and the phase slope or its square root. 

Best possible bounds on the rise time and overshoot are obtained when the amplitude 
function is itself bounded and approaches zero at least as fast as w ", n > O in the high 
frequeney region. These bounds may be evaluated readily from tabulated data. 

Introduction. ‘The transient response of a network is often substantially improved 
by the introduction of phase equalization. However, it should be determined whether 
such phase correction will actually provide sufficient improvement to warrant its use. 
Computations of the transient response are often tedious. If certain key quantities such 
as the 10-90°; rise time and the overshoot, of the unit step response, of the phase cor- 
rected network could be readily estimated, then the advisability of phase correction 
could be determined. Such readily evaluated bounds are presented here. 

Since phase distortion free networks have a linear phase characteristic, it will be 
assumed throughout this paper that the transfer function of the network in question 
is of the form 7(w)e ’** where T(w) > O and k is the constant phase slope. Low pass 
structures will be considered and for convenience 7(w) will be normalized so that 
7(0) = 1. The unit step response will be written as A (t). 

The 10-90‘; rise time of the unit step responses. Bounds on the 10-90% % rise 
time can be obtained by manipulation of the expression for the unit step response 


A(t) = + | TM) ain w(t — k) dw. (1) 
2 tJo 


Utilizing the fact that the unit step response of a linear phase network possesses odd 
symmetry about the 50°% point, which occurs at t = k, the 10-90% rise time r, can be 
obtained by solution of the equation 
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2r, [ T(w) sin [wr, /2] 
Ess anna iy, 


5 wT Jo ww; 

It is assumed here that the 10 and 90% points fall on opposite sides of the 50% 
point. If this is not the case, then the response will fall below 10% after ¢ = k. Such a 
highly oscillatory response is not suitable for most practical applications and will not 
be considered here. 

Equation (2) is difficult to evaluate. In fact, the rise time can, in general be more 
readily evaluated from a plot of the actual transient response. However, Eq. (2) can be 
used to obtain some readily evaluated bounds. 

Since 


smz | 
Ss i, 
M4 


T, = is | | T(w) te | ° (3) 
o Jt 


Often {> 7(w)dw does not exist, which renders inequality (3) useless. Even if the 
integral does exist, the lower bound may be unduly small. 

A somewhat different approach may be used to obtain a bound that exists in all 
physical cases where parasitic capacitance limits the high frequency response 

Apply Schwartz’ inequality to Eq. (2). 


: < taf [T(w)]}° dw [ E = 2) ae} 


i 
“a v0 « QW 


but, 
r sll (Ww, z ’ Tv 
| | } 7 | ‘he os 
‘ wr 17 
Thus, 
i6 fr Te x, 
" = on rT’ [T'(w) * day? ° }) 
7 a | | ) 


In many instances the bound presented by relation (4) will be much stronger than 
that stated in relation (3). Of course, there are also conditions where the converse is true. 

If some relatively weak conditions, which are obtained in practical situations, are 
imposed upon 7'(w), then a very readily evaluated greatest lower bound on the rise time 
may be obtained. 

Consider that 7'(w) < €(w./w)",n > Oforw > w,. That is T(w.) < € and T(w) falls 
off as w ", for w > w,. (For further discussion of such approximations see [1] 

Then the unit step response can be written in the form 


F T(w) sinawl(t — k a 
Pee ae 1 In w ®) deo + 8. 5) 


W 


where 


Us 
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Thus, if | 6 | is sufficiently small, the transient response can be assumed to be unaffected 
by frequency components greater than », . 

If the transfer function is such that frequency components greater than w, can be 
neglected, and, in addition, if T(w) is bounded so that T(w) < M, then the following 
procedure leads to a readily evaluated bound. 

From Eq. (2) 
=2[ rex/r,) ™ az. 


/0 


} 
5 


If w.r,/2 < mw, then the following bound is obtained 


b <2 MSiwo.r,/2) 
5~f 
or 
2 : vi 
r= Inv. Si[4x/(10M)]. (7) 
w, 
Where Inv. S7[x] is defined by the following relations if « = f% [sin z/z] dz = Si(y) 
then y = Inv. Sz(z). 


Note that Inv. S7[4x/(10M)] is single valued since WM > 1 and [4r/(1041)] < Si[2z] 
which is the least minimum of Si{[2]. 
In the derivation of relation (7) it was assumed that w,r,/2 < 7m. If for some T'(w), 
2 > z, then relation (7) will still be valid, since the result is given in the form of a 


W,/ 7 
lower bound and, thus, larger values of r, are allowed. 
It can be shown that values of r, given by the equality sign of relation (7) exist, and 
hence, this is the best possible bound. To do this, use the transfer function 
| lw<e 


T(w) = lim} M,e <0 <u. 


«0 





\ O,w > w, 


Although such a transfer function cannot be realized, it can be approximated arbi- 
trarily closely. 

Note that an increase in the value of M/ will reduce the bound on the rise time. 
However, since 7(0) = 1 < M, large values of M can result in large overshoots. 

The overshoot of the unit step response. The overshoot of the unit step response 
will be defined in the usual way 


0, = max. A(é) — 1. 


If the maximum value of A(¢) occurs at time t,, then manipulation of Eq. (1) yields 


7 : ) sims de. (8) 


I l 


In general, this equation cannot be used to obtain an overshoot since ¢,, the time of 
its occurrence is not known. However, a bound can be obtained. Since | sin 2/x | < 1, 
Eq. (8) leads to 
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0, < as + i —* | T(w) dw. 


y wT de 


Of course, this equation still contains the unknown time ¢,, . The bound attains its 


maximum value when ¢,, — k is a maximum. Since the response of a linear phase network 
is symmetrical about 50% point, and in a physically realizable case is zero for t < 0, 


the maximum value of t,, — k is k. Thus, 
[i a? «= 
3 ET Tw) dw. (9) 


The Schwartz’ inequality provides a second bound which may yield better results 
than the above if {5 7(w)dw is unduly large, or does not converge. 
This is 
( \1/79 
a aera kt 
0,.Ss\—- + 4 | [T(w) | daw? . (10) 
i) aa (2a Jo 


If T(w) is such that | 6 | given by relation (6) is negligible, so that the components 
of T(w) for w > w. can be neglected and, in addition, 7(w) is bounded so that 7'(w) < 1, 
then a bound on the overshoot which requires no integration can be obtained. 
Equation (8) leads to 


l aD abi r( Di es 
= 3 4. : | 1 7 -) isin x/2| dx. 
To bound this expression, consider a comb like T(w) such that 7[x/(t,, — k)] = 0 


if sin x < O and 7[x/(t,, — k)] = M if sin x > 0. The largest value of 0, then occurs at 
t, — k \|max = k. Thus, 


ee s 
0, S —s= t+ — {Su(r) + [S2(8r) — S2(2r)] + ---}. 
7 


The final term in this series depends upon the numerical value of w,k so that the complete 
form of the bound may be stated in the following way. 
If nr < wk < (n+ 1)z,n = 0,1, 2,3, --- and T(w) < M, the following bound on 


the overshoot is obtained. 


l Vx 
0o,<—-=+ > (—1)’ “Si(gr), n= 1,3, 5, 


2 © oo 
4 ’ ] 
P / M J ne ‘ g+lo-: \ . 

3 “34 < Si(w.k) + >» (-1) Si gr) 7 , n= 2. 4, «+. (11) 
= 7 1 / 
l | 

0S -=+ Si(w.k), n = U0. 
2 T 


It can be shown that these are the best possible bounds by using the comb like transfer 
function discussed previously. 

In order to aid in the computation of these bounds and to provide a means of rapidly 
estimating the overshoot the following table is included. 
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TABLE I 


> (—1)¢*! Si(gr) for odd n 

g-l 
n LX (-1)"* Si(gr) 

g=l 

I 1.8519 
3 . 1085 
2504 
3484 
4234 
4840 
5350 
5789 


Wo tw to bh bo bo to 


It should be noted that if the term | 6 | defined in relation (6) is not negligible, the 
bounds given by relation (10) can still be used if they are increased by | 6 | 
Conclusion. ‘Two bounds have been placed upon the 10-90% rise time of phase 


distortion free networks. These are 


} a & m 1 
i= “| | T'(w) ta 
o J0 
16 po ie 
: | Tle)? dof 
Zo \Je 


The following bounds have been placed upon the overshoot. 


and 


IV 


cf 
6, a + - | T(w) dw, 
2 «Jo 
1 k a= z. . 1/2 
4-24 {1 | [T(w)]} aus : 
— (<7 “vt / 


The shape of the transfer characteristic will determine which of these bounds is the 
best and should be used. 
If 7'(w) falls off so that its frequency components can be neglected for w > w, and in 


addition, if 7(#) < .W, then bounds may be obtained which require no integration. 


These bounds are the best possible and are for the rise time 


2 | dao 
Fa > Pe S 
2 in | tt, | 


Ve seiides 
0, < + — 2) (—1)*""Si(gn), n= 1,3,5,--- 
Ww o=1 


and for the overshoot 


bho 


} I Me Pc, a ‘ . 
a, 3 ae { Silo.) + p> (—1) siton} n= 2,4,6,--- 


, a bs 
2S -5 + — Sek), n= 0, 
T 
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where the constant 2 is determined from the relation 


nr <wk < (n + I)r, n= 0,1, 7.3, ++ 
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Appendix A. Application of the work of Zemanian to linear phase networks. 
Zemanian [2] has imposed certain restrictions on the unit step response of a physically 
realizable one port network in terms of the real part of the impedance function R(w). 
These conditions, when modified slightly, apply to the amplitude function of linear 
phase networks. 

Comparison of the expression for the unit step response of a network 


sin wt dw 


2 6° Rio 
A= = | R(w) 


+ 


with Eq. (1) indicates the similarity of the roles played by R(w) and 7T'(w). Some examples 
of these dual roles will now be given. It has been shown by Zemanian that: 
(1) If the unit step response is to be a monotonically increasing function of time, 
it is necessary that 
R(O) > Rw) forall w. 


(2) If R(w) is a monotonically decreasing function of w which approaches zero as w 
approaches infinity, then the overshoot is bounded by 


9 
0. < —1 + — Si(x) = 0.1790. 
T 


These may be written for the case of linear phase networks as 

(1) If the unit step response is to be a monotonically increasing function of time, 
it is necessary that 7(0) > 7'(w) for all w. 

(2) If T(w) is a monotonically decreasing function of w which approaches zero as w 
approaches infinity, then the overshoot is bounded by 

6. = a3 a ! Si(r) = 0.0895. 
ae ae 

In a similar way other restrictions imposed upon the real part of an impedance 

function can be carried over to the amplitude function of linear phase networks. 


REFERENCES 


1, P. M. Chirlian, An investigation of linear phase networks, New York University, Eng. Sc. D. thesis, 
Department of Electrical Engineering, May, 1956. Also, P. M. Chirlian, Bounds on the error in the 
unit step response of a network, Quart. Appl. Math. 16, 432-435 (1959) 

2. A. H. Zemanian, Bounds existing on time and frequency responses of various types of networks, Proc. 
IRE 42, 835-839 (May 1954). Also, A. H. Zemanian, Further bounds existing on the transient re- 
sponses of various types of networks, Proc. IRE 43, 322-326 (Mar. 1955) 























231 


ON THE APPLICATION OF DYNAMIC PROGRAMMING TO A 
CLASS OF IMPLICIT VARIATIONAL PROBLEMS* 


BY 
RICHARD BELLMAN (The Rand Corporation) 


AND 
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1. Introduction. A large and important class of variational problems has the follow- 
ing form. Given a vector equation of the form 


dx ; 

— = ¢(z, y), x(0) =e, (1) 
dl : 

where x is an N-dimensional vector, we wish to determine an m-dimensional vector y 

so as to minimize a given criterion functional 


a7 
J(y) = | h(x, y) dt, (2) 
where A(x, y) is a given sealar function. 
The vector y may be subject to constraints of the form 


r(x,y) < 0, ~= 1,2,--- ,q. (3) 


In problems involving “terminal control,” we meet the problem of minimizing a 


function only of the final state 
I(y) = klx(T)]. (4) 


A problem of this nature occurs when we wish to have the system in some specified 
state a,(7) at time 7’, without caring how the system gets there. This is usually an 
idealization, in the sense that a more realistic problem will involve a combination of 
a criterion of the type appearing in (2) together with some measure of the value of the 
final state. 

\s has been shown in some recent publications, see [1], where further references 
may be found, a variety of problems of this nature arising in economic and engineering 
control processes may be solved computationally by combining the theory of dynamic 


programming with modern digital computers. 
In recent years, problems of less explicit nature have become more frequent. Thus, 


’ control problem requires that 7 be chosen 


for example, what is called the “bang-bang’ 
so that the system tend to a specified equilibrium state as rapidly as possible; see [2]. 
The upper limit of integration is thus not predetermined, but is rather a function 
of the choice of the vector y. In place of a formulation in precise analytic terms of the 
type appearing in (2) or (4), we encounter an implicit criterion of the following type: 
“When x satisfies a set of conditions C, , C, , --: , C, for the first time, we want a 


given scalar function of x to be as small as possible.” 


*Received June 9, 1958. 
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A particular example of a problem of this nature, equivalent to one we shall discuss 
in more detail below, is one in which we require that a preassigned function be a minimum 
for the first value of JT for which x,(7) = a, , a given value. 

A number of quite interesting existence and uniqueness questions arise in conjunc- 
tion with problem statements of the foregoing kind. These will be discussed at some 
time in the future. Here we are interested in describing a technique which can be used 
to obtain computational solutions via the functional equation path of dynamic pro- 
gramming. 

The problem becomes of even more interesting nature if we insert some stochastic 
influences into the process. Let the governing equation be 


dx 


- 5) 
dt ‘ 


= g(r, y, 7), x(O) =, 


where 7 is a random vector. We now wish to minimize an expected deviation, or say 
the probability that the deviation exceeds a given critical value. 

Once again, let us point out that the rigorous groundwork for these questions remains 
to be laid. However, as we shall see below, we have a simple method for postponing 
this type of investigation. 

As is to be expected, certain simplifications are possible if the underlying equations 
are linear, i.e. of the form 


Gacy = ALe t+ Yn Tle; Y= C, (6) 


l 


and the criteria quadratic. We shall discuss these cases in some detail since they are 
of some importance in connection with the application of the method of successive 
approximations. 

Throughout, our aim will be to illustrate the applicability of the functional equation 
technique of dynamic programming to the computational solution of questions of this 
kind which appear in many ways to be outside the domain of the classical calculus of 
variations. 

2. Preliminaries. Since, as mentioned above, we are primarily interested in a 
computational solution of implicit variational problems of the type described in the 
foregoing section, we shall pose our problem in discrete terms. The recurrence relations 
we derive will then be ready for use in a digital computer. 

In place of the differential relation of (1.4), consider the difference equation 


Laser = Ola s Yas Ta) Zo = C, n=0,1,--- ,N. (1) 


One of the advantages of formulating problems in this fashion is that there are now no 
conceptual difficulties concerning the meaning of random functions or the existence 
of minimizing functions. In return, sometime or other we must show that the limit of 
the discrete process exists, and, preferably, yields the continuous process. For a start 
in this direction, see [3]. 

In order to illustrate the method in simple fashion, we shall consider a two-dimen- 
sional process, 

a(n + 1) = 2,(n) — y(n) — r,(n), ti(0) = ¢; , (2) 

xz.(n + 1) = gp[z,(n), x(n), y(n), ro(n)], 2(0) = ¢2 . 
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The aim of the process is to choose y,(7) and y.(n), subject to constraints of the form 
0<a, < yf") Sa, 0<b, < y(n) < bp (3) 


so as to minimize the expected value of [2.(m) — .»|°, where m is the “time’’ at which 
x,(m) 0. The r,;(n) are independent random variables with given distributions. 

The expected value is over the random variables r, and r, , where r, can depend 
upon the choice of y, and y, , but, in any case is subject to the condition that 


y(n) + 7r,(n) > a, > 0. (4) 
It follows that 2,(n) is steadily decreasing as n increases. 
The recurrence relation in (1) is valid until x,(n) = 0. Properly, we should write 
x(n + 1) = max [0, x,(n) — y,(n) — r,(n)]. (5) 


The process ends as soon as x, assumes the value zero. 
3. Functional equations. It is clear that the minimum of the expected value of 
lxv.(n) — xo] depends upon ¢, and c, and only upon these variables assuming all other 


t 


functions and distributions known and fixed. Let us then write 


fle, , 2) = min exp [x,(n) — xo)’. (1) 


We have 
f{(O0, c2) = (2 — 2)’, (2) 
and the principle of optimality, see [1], yields the functional equation 


fe, ,¢2) = min [exp fc, — y; — 17, , gle: , Co , Yo 5 T2)). (3) 
There is no difficulty in treating the case in which the distribution of random effects 
depends upon the decisions that are made. 
4. Probability of deviation. In place of mean-square deviation, let us consider the 
problem of determining y, and y. so as to minimize the probability that | 7. — 2» | > d. 


As above, let 


f(e, ,¢2) = min prob [| xz, — x | > d]J. (1) 
Then 
1(0,c.) = 1, C —2%|>d, (2) 
= (), leo — % | < d, 


») 


while f(¢, , c.) satisfies the same functional equation as in (3.3). 

5. Discussion of computational solution. In order to determine the function f(¢,, c2) 
using a digital computer, we employ a discrete grid in (c, , ¢.)-space. Let c,; assume 
only the sequence of values 0, 6, 26, --- , and c, a sequence of values 0, A, 2A, 
Since ¢, is monotonically decreasing as the process continues, we can use it as a “time” 
variable. Write 


f(k6, ce.) = file). (1) 


Then (3.3) may be written 
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f.(c2) = min [exp f,[g(k6, c2 , Yo , T2)]], (2) 


Yi.U2 Tit 
where p is determined by the condition 

p = [1 — yi: — 71)/4), (3) 
the greatest integer contained in (c; — y, — 7,)/6. 

Since g(k6, co , Yo , To) in general will not be an integral multiple of A, we can either 
take as its value the nearest integer multiple of A, as we did in (3), or we can use inter- 
polation, if more accurate results are desired. 

The value of f,(c.) is determined by the relation 


fo(c2) = (Co — 2)’. (4) 
Consequently (2) furnishes a recurrence relation which enables us to compute the 
function f,(c.) in terms of f,(c.) forn = 0,1, --- ,& — 1. We thus have a feasible compu- 


tational scheme. 

6. Deterministic process. Returning to a purely deterministic process, as specified 
by (1.1), we may wish to determine y so that x is in some desired state at some sub- 
sequent time. One way of attacking this problem is to treat the problem of minimizing 
{7.(T) — zx ]°, where T is the first time at which 2,(7’) has its desired value. The functional 
equations are as above, without the averaging over the random behavior. 

7. Linear equations and quadratic criteria. In general, the application of a straight- 
forward functional equation approach is limited by dimensionality difficulties in the 
sense that functions of three or more variables cannot be readily stored in a fast memory. 
Consequently, the techniques described above must be aided and abetted by successive 
approximations of various types, a subject which has been discussed elsewhere. If, 
however, the guiding equations are linear, and the criteria function quadratic, then the 
sequence of functions {f,(c)} will consist of a sequence of quadratic functions in c. These 
functions are determined once the coefficients are determined. As we shall see, reasonably 
simple recurrence relations exist connecting the coefficients of f,(c) with those of f,_,(c). 

Consider, to begin with, the problem of choosing the y; so as to minimize the ex- 


pected mean-square deviation 
7 
J7(y) = exp | cae —a,2x(T) — a) + z. (yy. Bu» } (1) 


Here 7 assumes the values 0, 1, 2, --- , B is a positive definite matrix, a is a specified 


state vector, x and y are related by means of the linear relations 
w, . = Ag, + Yy; 2 Tn» Lo = C, (2) 


where |r;} is a set of independent, random vectors with identical distributions. 

The process is assumed to proceed in the following fashion. We observe c, the initial 
state, and on this basis and the foregoing information, choose y, , the initial control 
vector. Then a random effect r, occurs, yielding by way of (2) a new state vector 
Ac + yo + 7 . The process then continues in this way, stage-by-stage, a ‘feedback 
control’’ process. 

Although this problem can be, and has been, treated by straightforward variational 
techniques, we shall treat it by functional equation methods. There is some merit in 














ee 
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doing this even in this case, and in addition we shall prepare the way for the following 
section devoted to a process of random duration. 
Define the new sequence of functions {f7(c)} by means of the relation 


fr(c) = min J7(y). (3) 
Then 
f.(c) = (ec — a,e — a), (4) 
and the principle of optimality yields the recurrence relation 
f,(c) = min exp [(y , Byo) + fr-1(Ae + Yo + o)], (5 
forn = 1,2, ---. 
Let us now show inductively that each f,(c) may be written in the form 
f.(c) = (c, M,c) + 2(b, ,c) + uy . (6) 
The result is obviously so for n = 0. 
Substituting in (5), we have 
f,(c) = min exp (Yo , Byo) + [Ae + yo + 10 , M,-1(Ae + Yo + 10)] (7) 
+ 2(b,-; , Ac + Yo + To) + U-s)- 
Taking expected values and using the result that 
min [(y, Cy) + 2(g, y)] = —(g, C™'9), (8) 


whenever C is positive definite, we see that f,(c) has the form stated in (6). Carrying 
through the calculations, we obtain recurrence relations connecting /, , b, and d, with 
M+ « Va4 ONG @. . 

8. Linear process of random duration. Consider now a system specified by the 
equations 


Un+i = Un — Tin 5 Un = Co 
n ( a) (1) 


Sass @ AZ, + BF Tos Lo = C, 
where u, and r,, are scalars, x, , y, and r, vectors. The process ends whenever u,, becomes 
zero or negative. 

The quantity 1,, 


is a uniformly positive random variable, so that the process is 
always finite. The control vectors y, are to be chosen so as to minimize the expected 


value of 
J(y) = [x(m) — a, x(m) — a] + YS (y , By), (2) 
hk 0 


where m is itself a random variable determined by the condition that it is the first integer 
for which w,, is negative or zero. 
Write 


f(co , c) = min exp J(y). (3) 


v 
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Then 


f(0,c) = (€—a,c-— a), (4) 


and 


f(e¢o ,c) = min exp f(¢o — Tio , AC + Yo + 1o)- (5) 


Assume, as previously, that ¢, can assume only a discrete set of values with a similar 
condition on 7:5 . Let, suitably normalized, cy take the values 0, 1, --- , and ry» only 


the range of values d; ,d; + 1, --+ , d.. Then, writing 


f(k,c) = f.o), c=), 1,2, +3>, (6) 


we may write (5) in the form 
\ 


f.(c) = min exp 4 ym Difn-( Ace + Yo + ro)? , (7) 
Vo To \ ty / 

where 

p; = the probability that r,, = 7. (8) 
The function f,(c) is identically zero for k < 0. 

Once again, it is easy to see that each element of the sequence {f,(c)} is a quadratic 

function of c, of the form 

fi(c) = (c, Mic) + 20 , oc) +m. (9) 
The recurrence relations connecting M, . b, , wu, with M,_, , b.-, , u,_; can be obtained 
from (7) in the way indicated above. 
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DETERMINATION OF CHARACTERISTIC VALUES* 


BY 
MARK LOTKIN 


Aveo Mfg. Corp., Wilmington, Mass. 


A method for computing the characteristic values of arbitrary matrices was recently 
proposed by the author in reference [1]. This method utilizes a sequence of unitary 
transformations which are designed to triangularize the original matrix, thereby pro- 
ducing the desired characteristic values along the diagonal of the triangularized form. 
Each of the unitary transformations is constructed in such a way as to reduce the norm 
of the upper-triangular part of the matrix. The basic function underlying this con- 
struction is a certain cubic polynomial whose coefficients depend upon the elements of 
the matrix to be reduced. 

In this paper there is presented a refinement of the cubic polynomial, which is shown 
to possess properties that make it superior to the previous polynomial, on the basis of 
theoretical and practical considerations. Theoretically, the modified approach is seen to 
become identical with the Jacobi method for symmetric matrices, in certain cases; 
practically, the modification has been found to lead to more rapid convergence, at least 
for a considerable number of certain matrices that were subjected to both techniques. 

A brief résumé of the relevant equations inherent in the procedure seems appropriate 
here. Let us assume, then, that the sequence of transformed matrices has reached the 


pth stage A, , p = 0,1, 2, --- , Ag = A, and that it is consequently desired to construct 
As, = T,' A,T, , where the unitary matrix 7, has the ‘‘norm-reducing”’ property, 
1.e., if 


denotes the “upper-triangular’’ norm of A, , and M,,, denotes the corresponding quantity 
for A,., , then 


M,., — M, < 0. (1) 

Let the element b = | b | exp 78, located in the (7, j) the position ¢ < j, of A, be the 
“pivot” for the next transformation 7,,,, , and let the elements in the position (2, 2), 
(j, j), and (j, 7) of A, be denoted by a = | a | exp za, d, and c, respectively; in general, 


let a,?’ = | a,2’ | exp (ta,2’). 
It was shown in [1] that, for R ¥ 0, 


M,.1 — M, < H(R, 0) = RF(R, 6) (2) 


Il 


with 


F(R, 6) = C3 +- CR? + Ck + Co (3) 
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C, = —2\/c| [| d| cos(@+y7— 6) — |a| cos(0@+ y — a)], 


C,=|!d-—al’—2|b c| cos (286 — B+ y) + >» (la? |? + | ajr’ |’), 


ponies 
Co = 2/6! [| d| cos (6 + 6 — 8B) — | a! cos(0+a — B)] 
be [| aie’ | | aj’ | cos (0 + ayy — ay) 


— | a? | | ap? | cos (0 + a; — a;)). 
Any solution (R, 6) of the system 
F(R, 0) = 0 (4) 
dF /00 = 0 (5) 


then guarantees that M/,,., < M, . Having determined R, 6, the elements of the trans- 
formation matrix 7’, are found by means of 


r= (R?+ 1)’, senr = sgn R, t = r exp 78, (6) 


and the new elements of A,., are determined from the relationships (8) through (15) of 
reference [1]. 

In addition to the relationships (16) between the elements of A,,, and A, , stated 
in [1], there exist further identities between these elements, of value in the actual per- 
formance of numerical calculations; some of these are exhibited below. 

A short calculation shows, for example, that 


a, |? +10, 17? +e 1? + {1d (/?=/al/?+/b)] 7+ /e|7*+/d)?, (7) 

by virtue of r°(1 + R°) = 1. For the same reason, 
lax / +tlan ? = lan |? + | a, |’ (8) 
au; | + | a; |? = | a; |? + | a; |’ (9) 


for 1 =< kh <= nk x i, 3. 
Further, whenever 6 = 0, 


b, —c, = b—-e. (10) 


As stated in the introductory paragraph, the determination of R, @ was previously 
based on Eqs. (4) and (5). 

Now a possibly larger decrease in the norm M than indicated by (4) and (5) may be 
obtained by choosing R such that 


M,-—-M< min A(R, 6). 


. 
The values of 2 at which min H occurs must then satisfy 
0H/dR = 4C;R° + 3C,R* + 2C.R + Cy, = 0 (11) 


0H/d6 = 0. (12) 
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If R = R,, satisfies (11), then clearly 
M — M, > R;(3C,R2, + 2C.R, + C,). 


Now for fixed 6, H(0O, 6) = 0, (@H/AR)o.6 = 0, (0°H/AR*),.5 = 2 C, . Thus min 
H < 0 whenever C, # 0. If Cp = 0, but C, < 0, then still min H < 0. Only if Cp = 0, 
Cc? — 4C,C, < 0, is min H equal to zero. Thus in general (11) may be considered superior 
to (4) for the reduction of the super-diagonal norm. If Eq. (11) has, for fixed @, two nega- 
tive minima, then we choose for the transformation z = R exp 76 that root R for which 
min min H is assumed. 

The superiority of (11) may be deduced also from the study of certain second order 
matrices, which are obviously of basic importance in this problem. 


I. Let us consider 


A= [@ b] (13) 


with a = d. If A is skew-hermitian, then a = d = 0. It may be assumed that be ¥ 0. 


Since now a = 4, (11) becomes 

4\c|Rijc|R? —|b}) =0, 
so that R = | b/c |'”, and min H(R, 6) = — | 6 |?.With z = | b/c |'” exp 2™* (8 — y)i, 
and r° = |¢ b| + | c¢ |), the transformation equations lead to 


a, = a+ (be)'” 


d, = a — (be)'” 
b, = 0 
c, = (b/| b |)\(\e| — | b)). 


The expressions for a, , d, are, naturally, the exact roots for the matrix (13), with 
a = d. Thus MV, — M = — |B’. 

II. Again let us consider matrix A as defined by (13), now subject to the following 
conditions: 
(i) a= Fi) 


(il) 8 + ¥Y 2a 


(iii) @ = 2°‘ (8 — y). 


Then 
C, = {el 
C, = —2/c|(|d| —|af) 
C, =({|d|—!a|)*’ —2! be} 
Co =2|b((|/d| —|/a)), 


whence 


SR) = [ic|R? —Cd| —Ja)R—|b|f2le|R—Gad|—|fa)] (14) 
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H(R) =Rile|R?-—(/d| —|a)R—2] bl]ff{e|R—({d|—|a))]. (15) 
If|c| R’ — (|}d| —|a]) R —|b| = 0, then A(R) = — |b |’. If, however, 2 | c | 
R—(|d|—|!a|) = 0, then H(R) = | c| R* (| c| R° + 2|b]) > O. Consequently, 
min H(R) = — | b|° is assumed at the roots of |c| R° — (|d| —|a|)R —|6b| = 0. 


Therefore, b; = 0, so that the characteristic values of A appear immediately as a, , d; . 
III. Now let the nth order matrix A be hermitian, i.e., a;; = d;; , a;; real. 
Then it is seen that the cubic polynomial F(R, 6) of (2) becomes 


C, b 
C, -2!b!|[\d!cos(@+y— 5 —!a!cos(6+y¥7—a)] 
C, = ld—al* — 2|0b/? cos (20 — B+ 7) 
Co = 2/6! [| d| cos(@ + 6 — B) — | a! cos(@ +a — 8B)]. 
for the choice of 6 = 27" (8 — y) above expressions become 

C,=1b/* 

C, = —2|b| (d — a) 

C, = (d — a)’ — 2! b!? 

Ce = 21 1 a = 2), 


and, consequently, 
H(R, 6) = RF(R, 0) = R{| b| RF — (d-—a@dR — 2| b{][| b| R-— d — a)] 
S(R, 0) = 0H/dR = 2[| b| R? — (d—a@R — | b | }[2| b| R — (d — a)]. 
Thus again H(R) = — | 6}? for the roots R of the quadratic factor of S(R, 6). This, 
naturally, implies again that b, = 0. 
IV. The matrix 


rt £& 4 
Be !O0 2 1) (16) 
: @ 3 


has been mentioned as one which defies direct treatment by Greenstadt’s method [2], 
as well as by the method of reference [1]. However, it is stated in [3] that by applying a 
transformation to B which effects a rotation through @ = 7/4, the transformed matrix 
becomes tractable by Greenstadt’s method, and that in twelve cyclically executed 
annihilations of the respective pivot the matrix B becomes triangularized to a sufficient 
degree of accuracy. 

It will be seen that the same preparatory rotation through 0. = 7/4, followed by 
two transforms of the type determined by (11), exactly diagonalizes the matrix (16). 

Since the transpose B’ of B has the lower superdiagonal norm, we subject B’ rather 
than B to the sequence of unitary transformations. The preliminary rotation is effected 
by 

/ cos 6 O —sin 4, 
T=| 0 | 0 (17) 


|sin@, O COS Oo. 
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leading to 
| 3/2 cs 19 
B, = T'B'T = . 1 -2'? 
L-1/22'°7 1/2 


Let us choose next the element a{}? = 2°'”’ as the pivot b. Then case II discussed previ- 
1/2 


ously is found to apply. With z = 2 one obtains 
2 «0 0 | 
B,=|0 1/2 —3'”%/2}- 
1/2 


The superdiagonal norm M(B,) = 5/4 has thus been reduced by the amount 6? = 27' 
to M(B.) = 3/4. 





Next we take a{3) = —3'’/2 as the pivot b. Here case I obtains. Therefore, R = 1, 
2 =7, and 
2 0 0 | 
B, =|0 (1/2) + i/3'”) 0 (18) 
'e 0 (1/2)(1 — i/3'”)_] 


so that B has actually been reduced to diagonal form. 
The diagonalization of B” has thus been achieved by subjecting it to the unitary 
transformation B, = P~'B’P, with 





a7? 9 -2"*1 2/3)" -3"" of1 © 0 | 
P=! 0 1 0 s ow ¢@ | 02'7% -2'? 
ot 9 I 0 0 1| aes” -9"%! 
3°? -(1/2) — (2-37) "i (2-3'?)* + 7/2] 
= |3°'” on"; -3°" . (19) 
3°? (1/2) — (2-:3'”) "i (2-3')" — 7/2] 


According to a theorem of Toeplitz (see, e.g., [4]) a matrix M can be reduced by 
unitary transformations to the diagonal form if and only if the matrix M is normal: 
M‘’ M = MM", where M‘’ denotes the conjugate transpose of 7. Thus the matrix B 
is seen to be normal. The interesting question then arises whether a class of normal 
matrices of which B is a member can be reduced to diagonal form by the general technique 
of this paper. This question can be answered in the affirmative; the results will 
be published elsewhere. 

V. While the choice of a particular value of @ may be appropriate, in special situations, 
in general the condition (12) may have to be considered, for optimum results. In the 
example discussed here the values of 6 = 0, 1/6, 2/3, --- 54/6 were applied to each 
pivot, which was always chosen to be the element of largest modulus. For the trans- 
formation z = R exp 76 that pair (R, @) was selected for which min H(R, @) is assumed. 

The following matrix is taken from [1]: 
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l 0 -2 
salad i —] ae (20) 
2 ] 0_| 


its characteristic values are —2, 1 + 27. Using the method of [1], eight iterations of A 
produced the diagonal terms shown under (33) of reference [1], viz. 


—1.97139 — .07432:, 95014 + 2.119927, 1.02125 — 2.04650:. 
The use of the refinement based on Eq. (11) led to 
— 1.99489 + .01100:, .99161 + 2.001587, 1.00328 — 2.012587, 


clearly a considerably improved result over the previous one. The super-diagonal norm 
at this stage had been reduced from 8.000 to 1.424 X 10°”. 
VI. The matrix 


11ii1i1 41 
6 7 8 9 10 11 


is one of a sequence of non-symmetric matrices of extremely bad ‘“‘condition” [5]. Matrix 
l 


C is nearly singular; the absolute value of its determinant is (31 05 2236 7232 - 10°°)"' = 
3220 3799 - 10°**. It is well known—see, e.g., Todd [6]—that certain calculations with 
these matrices, such as inversions, determination of characteristic values, etc., suffer 
from ‘‘numerical instability’. For machines with twelve decimals, employing floating 
point arithmetic, attempts to invert matrices of even the eighth order have been doomed 
to failure. 

The characteristic values \,, of C, calculated by means of a determinantal method, 


and arranged in order of decreasing magnitude, are listed in Table 1. 


TABLE 1. 


Characteristic values of C. 





n An 

1 .2132 3763 x 10! 

2 ~ 2214 0681 x 10° 

3 — 3184 3305 x 107 
{ — 8983 2330 x 10-8 
5) — .1706 2788 x 10 


1397 4990 « 10-¢ 


| 
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It is seen from this table that 4, - \, --- As = .3220 3799 x 10°"* = det C, to eight 
significant figures. 

The reduction technique described above was programmed for the IBM 704 machine. 
Single precision arithmetic, and floating point with eight significant figures was used. 
At each step, the super-diagonal element of largest absolute value was taken as the 








pivot b = a,; of the next transformation. Some of the results are shown in Table 2. 
TABLE 2. 
Triangularization of Matrix C 

p M, gf a\™ 
0 .494 XK 10° 1.0000 0000 .9090 9090 « 107+ 
10 .236 X 10°? .2124 5647 X 10! .1250 4932 «K 107? 
20 .553 X 10°* .2130 8000 x 10! .2123 4414 x 10° 
30 296 X 10-5 .2132 2888 x 10! — .6031 1972 x 10-4 
10 385 X 1076 .2132 3244 xX 10! — .7919 3310 x 10-* 
50 .616 X 10-7 .2132 3678 x 10! — .2023 3858 x 10-5 
60 544 X 1078 .2132 3763 X 10! — .2646 2287 x 1075 
70 980 < 10-9 .2132 3705 x 10! —.1103 9584 x 1075 
80 149 x 107° .2132 3751 X 10! — .1097 4208 x 10-5 
90 239 K 107% .2132 3772 X 10! —.1057 6060 x 10-5 
100 .448 X 107" .2132 3763 XK 10! — .1451 8099 x 10-6 
110 .647 X 107" .2132 3765 X 10! — .1315 6736 x 10-6 
120 .143 X 10°” . 21382 3766 xX 10! — .1652 8475 x 10-6 
130 .245 XK 10-7" .2132 3765 xX 10! — .1634 9713 « 10-6 
140 581 x 10-" .2132 3765 xX 10! — .1371 2160 x 10-6 
150 834 xk 10° .2132 3765 x 10! — .1409 3837 x 10-6 


In this table p denotes the number of iterations, JJ, the super-diagonal norm of A, , 


as’, a\™ the diagonal elements in A, of largest and smallest absolute value, respectively. 


The diagonal elements at p = 150, arranged in order of decreasing magnitude, are: 
.2132 3765 * 10' 
— .2214 0677 x 10° 
— 3184 3361 « 10°' 
— .8983 2775 x 10° 
—.1705 5897 x 107° 
.1409 3837 x 10°° 


Thus the dominant characteristic value is determined at this stage to about seven 
significant figures, while the smallest one is known to about three. The rate of decrease 
of M, , from .56 X 10° to .8 X 10°"’, would seem to indicate a “linear” type of con- 
vergence. It is obvious that further improvements of the results will be achieved once 
a number of basic routines that are presently in the program have been sharpened. Such 
routines are concerned with the conversion of numbers from the decimal to binary 
system, the calculation of trigonometric functions, the location of roots of polynomials, 
and other operations required in the method. 

VII. Among the many other matrices that have been reduced satisfactorily we 
mention here the Hilbert matrices. These are symmetric matrices whose elements are 
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a;; = (+ j-—1)",2,j = 1, 2, 3, --- . Some results for the eighth order matrix are 
exhibited in Table 3. 
TABLE 3. 


Characteristic values of a Hilbert matrix 


p M, a’ a:; 

0 882 x 10° 1.0000 0000 .6666 6667 « 107 
10 262 X 107! 1693 0662 x 10! .4273 0205 xk 10°? 
20 .388 X 1073 .1695 9118 X 10! 4030 7163 « 107° 
30 543 x 1075 1695 9389 x 10! .3309 0313 «Kk 10~¢ 
40 Mi x iv? .1695 9390 x 10! .1518 9860 « 1074 
50 .123 xX 10-° .1695 9391 x 10! .6202 6355 x 10-5 
60 211 <x 10-8 . 1695 9391 x 10! 3315 4293 x 1076 
70 364 X 107% .1695 9391 x 10! .¢117 3257 X& 10-8 
75 .424 k 10-4 1695 9391 x 10! 6802 6828 « 1075 


The intermediate diagonal elements in A;, are: 

.2981 2524 X 10° 
.2621 2851 XK 107 
.1467 6944 & 10~° 
.5437 2030 x 10°° 
.1297 1307 x 107° 
1589 9581 x 107° 


The trace of the matrix, which is theoretically equal to the sum of the characteristic 
values, is 2.0218 0042. The sum of the diagonal elements at p = 75, on the other hand, 
is found to be 2.0218 006. The well-known Givens method for the characteristic values 


of real symmetric matrices results in a corresponding value of 2.0218 002. 
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SOME PROPERTIES OF RATIONAL TRANSFER FUNCTIONS 
AND THEIR LAPLACE TRANSFORMATIONS* 


BY 
ARMEN H. ZEMANIAN 


College of Engineering, New York University 


Introduction. The behavior of a fixed, linear system at a pair of output terminals 
due to a signal impressed at a pair of input terminals is determined by the transfer 
function between these two terminals. More precisely, this transfer function Z(s) is 
the ratio of the Laplace transform of the output function of time to the Laplace trans- 
form of the input function of time. Moreover, if the system is lumped and finite, the 
transfer function will be a rational function of the complex frequency variable s = ¢ + je 
and may be written as the ratio of two polynomials in s 


_” 1 “*-e € y 7 
s" + 4,-18 +---+a, , N(s) (1) 


OF SF ee ok ie 
The constant multiplier K and the coefficients a; and 6, are all real numbers. It will be 
assumed still further that the system is a stable one. That is, if the input signal is removed, 
the output signal becomes arbitrarily small with increasing time. This restricts the 
denominator D(s) to being a Hurwitz polynomial (a polynomial whose roots all have 
negative real parts). Finally, the transfer function of any physically realizable system 
must become arbitrarily small as the magnitude of the complex frequency s becomes 
arbitrarily large. This is due to the stray capacities and parameter dissipation that 
must exist in such a system. Thus, this discussion will be restricted, in addition, to those 
transfer functions that have more poles than zeros (m > n). This insures that the Fourier 
transform, relating the transfer function Z(jw) to its corresponding unit impulse response 
W(t), exists. Only those transfer functions which satisfy the aforementioned conditions 
will be considered in this paper. Therefore, all unit impulse responses will be a finite 


sum of terms of the form, 
At’e *' cos (Bt + y), 


where A, a, 8, and y are constants with a > 0 and uz is either zero or a positive integer. 

In a previous paper [1], certain classes of these transfer functions have been in- 
vestigated and many restrictions on the transient responses corresponding to such 
transfer functions were shown to exist. These classes are related to the class of positive 
real functions for which m — n equals one. However, the sum of all these classes does 
not include every transfer function of the type considered here. A transfer function that 
belongs to one of these classes is called a class k function and the integer k identifies 
the particular class. Throughout this paper, the positive integer k will represent the 
number of poles in excess of the number of zeros. 

The purpose of this paper is to determine other properties of the class k functions 
and then to extend the restrictions on their corresponding transient responses to the 


*Received June 25, 1958. 
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transient responses of any transfer function. That is, a means will be developed for 
determining bounds on the transient responses of any transfer function, whether or 
not that transfer function is a class k function, by relating it to a class k function. 

The class k functions. For this investigation, the form of the definition of class k 
functions will be altered somewhat from that given in [1]. It should be emphasized, 
however, that these two definitions are the same in that they define exactly the same 
classes of functions. The new definition leads to simpler expressions in this paper. 

Consider the function Z,(s) which is obtained by successively integrating the transfer 
function Z(s) according to (2) 


Z.(s) = | diy | d 8-2 °° | Z(8) d 8 « (2) 
The integer g is restricted to being less than or equal to the integer k — 1. The right- 


hand side of (2) is, in general, a multivalued function and, if all the poles of Z(s) are 
simple, it may be written as 

oe (s + p;)* | — 1 

(—1)°K 2 i, ° (q—1 yt in (s + Ge) = =] 
where the p; are the poles of Z(s) and the k,; are the residues of these poles. [Z,(s) will 
be single-valued when Z(s) has only one pole of multiplicity m.] For the purposes of 
this paper, only one branch of this multivalued function is needed. Therefore, given 
any Z(s) having at least two more poles than zeros, branch cuts consisting of straight 
line segments connecting all the poles will be assumed. In addition, these straight line 
segments will be so chosen that any pair of poles can be connected by a path traced 
exclusively along these segments. That is, it will be assumed that these branch cuts 
form a tree with the distinct poles at its vertices. The number of these branch cuts 
will be one less than the number of distinct poles of Z(s). Finally, it will be assumed 
throughout this paper that the paths of integration in (2) never cross any of these 
branch cuts. Under these conditions, (2) defines a single-valued function Z,(s). 

Representing the real and imaginary parts of s, and Z,(jw) by s, = o, + jw, and 

Z,(jw) = R,(w) + jl,(w), the real part of Z,(jw) is given by the following expressions. 
For q odd, 


a «x 


R(w) = (—1)°*"! | dw, | dws. *** | I (wo) dw (3) 


yw Jwa-1 Wi 


and for q even, 


a & 


Rw) = (=9"J dw, | dw. *** | R(wo) dary , (4) 


q@-1 J @i 
where PR(w) and 7(w) are the real and imaginary parts of Z(jw), respectively. 
The class k functions may now be defined as follows. 
Definition. The transfer function Z(s) will be called a class k function, where k = m — n, 


Ss 


if Z,.-,(s) is a positive real function in the half plane o > 0. 

An immediate consequence of this definition is easily obtained. The necessary and 
sufficient condition for Z(s) to be a class k function is that R,_,(w) > 0 for all w. This 
follows from the definition since Z,_,(s) is an analytic function and, moreover, Z,_,(c) 


is a real function for o > 0. 
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Furthermore, class k functions may be generated from the class 1 functions by 
differentiation. More generally, if Z(s) is class k, then (— 1)"d"Z/ds" is class (k + h). 
Another elementary property of these functions is that aZ,(s) + 6Z,(s) is class k if 
Z.(s) and Z,(s) are each at least class k and a and 8 are positive numbers. 

Realizability conditions on the impulse response. It will be assumed throughout 
this paper that all driving functions of time are applied at ¢ = 0. The unit impulse 
response W(t) is, therefore, related to the transfer function Z(s) through the Laplace 


transform, 


Z(s) = [ Wie" dt, (5) 
0 
where o is greater than the real part of any pole of Z(s). Moreover, the function Z,(s) 
may be obtained by integrating (5) under the integral sign qg times. The resulting ex- 
pression (6) is valid for g < k — 1 [2, Theorem 17, p. 273). 
ad ” Ee xe 
Zs) = | pee" a (6) 
A well known result of circuit theory is that, if Z(s) is a positive real function, then 
it is expressible in terms of three positive definite quadratic forms which are related 
to the total stored energy and power dissipation in the network [3, Chap. 4]. There is 
a somewhat similar result in the theory of Laplace transform which will place restrictions 
on the unit impulse response corresponding to a positive real Z(s) [4, Theorem 6]. More 
precisely, Z(s) is a class 1 function if and only if W,(¢) is a positive definite function 
where W’,(é) is the even function given by 


(wit) for «>0, 
W(-?t) for t< 0. 


Wid) = 


By definition, W,(t) is positive definite if the Hermitian form 


n 


je > W At; — t)x,a; (7) 


a=1 y3=1 


is non-negative for all values of the real numbers ¢, , t2 , --- , t, and 2% , %, +++, 2X, 


and for all n. 
Applying these known results to expression (6), the following realizability theorem 


is immediately obtained. 
Theorem 1. The necessary and sufficient condition for Z(s) to be a class k function is 
that W.(t)/ | t |"~' be a positive definite function. That is, Z(s) is class k if and only if 


.<. Wil, — 
yy Pe aa, 20 (8) 


t=] j=] 
for all values of real numbers t, , tz , +--+ t, and 2, , %2, +++ , x, and for all n. 
There are several necessary and sufficient conditions for a function to be positive 
definite [4 through 8] which could similarly be used for realizability conditions. However, 


all of these are, in general, too complicated to be useful as a practical test. 
A condition which is sufficient (though not necessary) and at the same time quite 
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to 
= 
¥ 2) 


simple to apply may be obtained by applying Theorem 5 of [9] to the function W(t)/t*~* 


for the case where the integer m, defined in that theorem, equals one. This yields the 
following sufficiency test for a class k function. 
= ‘ ‘ ‘ - 22% k—1 » ‘ a ‘ 

Theorem 2. If the second derivative of W(t)/t is non-negative for t > 0, then Z(s) 
will be a class k function. 

A necessary but not sufficient condition on the impulse response, in order that. it 
be the Laplace transform of a class k function was given in [1, Theorem 9]. A new proof, 
which is considerably simpler than the one given previously, may be constructed by 
using an operation transform pair of the Laplace transform. 


Theorem 3. If Z(s) is class k, then, for t > 0, 


1 


wi | < —K! (9) 


- ( — J) 
Proof. By the definition of class k functions, the function Z,_,(s) is positive real and 


the first term of its inverse series expansion is 


A (10 
(kk — 1)fs’ 


where k = m — n. It has been shown that such a positive real function has a unit impulse 
response which is bounded by the coefficient of the first term in its inverse series ex- 
pansion [9, Theorem 1]. Moreover, each integration in expression (2) corresponds to 
the division of W(t) by the time variable ¢. This transformation of the operations of 
complex integration and division by ¢ is a well-known property of the Laplace trans- 
formation [2, Theorem 17, p. 273]. Thus, for ¢ > 0, 


| Wie) K 
es < ” ; 
t — (k — 1)! 
which is the desired conclusion. 

By use of the superposition integral, the previous result may be extended to those 
transfer functions that are products of class k functions. This is applicable, for instance, 
to a system of cascaded amplifiers each of whose interstage networks is class k. It may 
also be applied to the transfer functions of ladder networks [10]. 

Theorem 4. If the functions A(s), B(s), --- , H(s) are class a, class b, --- , and class h, 
respectively, and their constant multipliers are K, , K,, --- , and K, , respectively, then the 
unit impulse response W(t) corresponding to the product of these functions, A(s) B(s) 
H(s), 1s bounded by 
» AK, biden K +h+cootf—}) P 
W(t) | < ——— * e rT (11) 

eee +. hw 1 
fort > 0. 

Proof. Let W,,(t) be the inverse Laplace transform of the product A(s)B(s). By the 

superposition integral, 


W,,,(2) < | | W.(r) || W(t — 2) | dr. 


70 
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But, since A(s) is class a and B(s) is class b, 


er K.K, — 1 
WAO1S Gd — mJ a ae 


Integrating this expression by parts b — 1 times, 


= t > = 

KK, a+b-2 K.K, (a+b-1) 

1 T dr l ; 
- #0 


. < los ee 
Wald) —~(a+b—-2 (a+ b— 1)! 


Applying this process repeatedly to the sequence of products A(s)B(s), A(s)B(s)C(s), 

- , A(s)B(s) --- H(s), expression (11) will be obtained. 

A relationship between an arbitrary transfer function and the class / functions. 
One of the objects of this paper is to develop a method of determining bounds on the 
transient responses of any transfer function satisfying the restrictions enumerated in 
the introduction. This will be accomplished by relating such a transfer function to a 


class k function as follows. 
Theorem 5. For any given Z(s), a real number c can be found such that Z(s + c) is 
class k. Moreover, if Z(s + c) is class k, then Z(s + d) ts class k for all d greater than c. 
Proof. By the definition of class k functions, Z(s) will be class k only if Z,_,(s) is a 
positive real funetion. Now consider the loci of Re [Z,_,(s)] = 0. Since Z,_,(s) behaves 


as 1/s in the neighborhood of s = ©, only one such locus exists in a sufficiently small 
neighborhood of s = ©. Moreover, this locus is asymptotic to a line parallel to the 


imaginary axis. Thus, for c sufficiently large, the line s = c + jw, which is parallel to the 
imaginary axis, will not intersect any locus of Re [Z,_,(s)] = 0. Furthermore, Re [1/s] > 0 
for o > 0. Therefore, Re [Z,_,(s)] > 0 foro > c. 

The second statement of this theorem follows from the fact that a positive real 
function of a positive real function is positive real. 

Theorem 6. If the real number c is such that Z(s + c) is class k, then c cannot be less 
than (a,-; — Bm—,)/k, where a,., and b,,., are the coefficients indicated in (1). Moreover. 
c cannot be less than the real parts of all the poles of Z(s). 


Proof. Integrating the inverse series expansion of Z(s) according to (2), a series 
expansion for Z,_,(s) may be found 
- K [{k a,-, — b,, 
Z,-\(8) = = fe EE oes J 
k! \s s 
It can be seen from the first two terms that a locus of Re [Z,_,(s)] = 0 is asymptotic 
to the line, 


$= Gn hoe Pmt + jw, 


for w sufficiently large. Moreover, such loci are continuous curves which pass through 
the singularities of Z,_,(s). Sinee all the singularities of a positive real function are 
excluded from the right half s plane, the constant c cannot be less than (@,-,; — Dm-1)/k 
if the line s = c + jw is not to intersect a locus of Re [Z,_,(s)] = 0 and if Z,_:(s + c¢) is to 
be analytic in the right half s plane. Moreover, the singularities of Z,_,(s) occur at the 
same points as the poles of Z(s) so that ¢ cannot be less than the real part of any such 


pole. 
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One physical interpretation of the addition of a positive constant ¢ to the complex 
frequency variable s is that dissipation is added uniformly to all the reactive elements 
of any network which realizes Z(s). In particular, every inductance L, has a resistance 
R,; = L, ¢ inserted in series with it and every condenser C;; has a conductance G; = C; 
¢ connected in parallel with it [3, pp. 706-708]. 

A class k function Z(s + c) can be obtained from a given Z(s) in the following way. 
Consider the polynomial in the numerator of Z(s) 


N(s) = s*" + a,-:8" '+---+as+a. (12) 


The polynomial N(s + c) may be determined by setting u equal to s + c and expanding 
N(u) into a Taylor series around the point » = c 
n—-1 m—-lar ry 
a ‘ 8 d" N(s) dN(s ae ' 
N(is+ec)=s8 toy | + -°- tee + N(c). (13) 
(n— 1)! ds _— a +... 
A simple method exists for obtaining the coefficients of N(s + c) from the coefficients 
of N(s) without having to perform the differentiations indicated in (13) [11, pp. 52-54]. 
This procedure is especially useful when a numerical value for c is used. Divide N(s) 
by s — c until a constant remainder r, is obtained. 


N(s) ; . 
= = Q,_.(8s) + —-- (14) 
s—c s—ec 
The value of N(c) is r, . Dividing the quotient Q,-,(s) in (14) by s — c until another 
constant remainder r, is obtained, the value of the first derivative of N(s) at s = c will 


be found. 


Q,-1(8) = Q,_.(s) +—, 


=-¢ eee 


Ni) " 
pm ” =. 


Repeating this process, the value of the 7th derivative of N(s) at s = c divided by factorial 
4 is found to be equal to the 7th constant remainder r, 


1 d'N(s) _ 
ii ds |. 


The verification of these relations may be found in [11, Chap. 3]. 
Denoting the remainder, obtained from the denominator D(s) through the procedure 
indicated above, by r’ , the expression for Z(s + c) may be written as 
iid, - tris" |) tes +78 + 1% . 
Zis+ec=K oe eae m—-1 | en eee ee ee ee (15) 
$s 7 Tans8 T °** +18 + % 
Therefore, for any given Z(s), a c may be determined which produces a class k function 
Z(s + c) by applying the tests for a class k function [1, Theorems 4 and 6] to expression 
(15). 
If Z(s) has only one finite zero, an application of [1, Theorem 5] yields values for c 
which insure that Z(s + c) is class k. In this case, Z(s + c) has the form 
Stet a -s+et+a 


Zs+cec)=K ™ + ?_3” ry ine oe Ws re = K “esta : (16) 
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where r’,_, = mc + b,,-,, rf = dD(s)/ds at s = c, and ri = D(c). For m even, Z(s + c) 
will be class m — 1 for all values of c greater than or. equal to (a) — bn_,)/(m — 1). 
For m odd, this theorem places two restrictions on c; it must be greater than or equal 
to (ap — bn-:)/(m — 1) and greater than (rj/r{ — ao). 

A simple graphical method for determining this constant ¢ to insure a class k function 
Z(s + c) is based on the fact that, if the phase angle of Z(jw) equals an odd multiple of 
x/2 exactly k — 1 times for k odd or if this angle equals a multiple of + or zero exactly 
k — 1 times for k even and if dg/dw < 0 at w = 0, then Z(s) is class k [1, p. 282]. The 
procedure is quite simple if N(s) and D(s) are factored 


(s — p,)(s — wo) «++ (8 — w,) 
—— 2: ~ . 17 
(s — pi)(s — fx) -°* (6 — 9p.) ( 





Zs) = K 


Following the convention that the phase angle of any factor in (17) remains within 
the interval (37/2, — 2/2), the phase angle of (s — y;) will be denoted by y, and the 
phase angle of (s — p;) will be denoted by 6; . The phase angle ¢ for Z(s + c) is then given 
by 


e= Lu LH. (18) 


Now the phase angle of the factor (s + ¢ — u,) is the angle of the vector from the point 
u; to the point s + c. Consequently, the variations of g along a series of lines parallel 
to the imaginary axis in the s plane may be obtained graphically by applying (18). For 
some sufficiently large value of c, the aforementioned sufficiency criterion on the values 
of ¢ for s = jw will be fulfilled and, for this value of c, Z(s + c) will be class k. 

Bounds on the impulse response of any transfer function. Now that any transfer 
function Z(s) may be related to a class k function Z(s + c) by an appropriate choice 
of c, the restrictions on the impulse responses of class k functions may be extended to 
the impulse response of Z(s). 


Theorem 7. The unit impulse response W(t) corresponding to any transfer function 
Z(s) ts bounded by 


[wo|s =; (19) 


where the constant c is such that Z(s + c) ts a class k function. 


Proof. Consider the inverse Laplace transformation between Z(s) and W(t) where 
the path of integration is along the line s = ¢ + jw 


+j@ 


: [ Z(s)e*' ds. (20) 
7) Jc-—j@ 


— 


W(t) = 


Making the change of variable s = v + c, (20) becomes 


W(t) = e'WA0, (21) 
where 
WA) = 5 [ Zw + oe" de. (22) 
47) J-jo 











ARMEN H. ZEMANIAN (Vol. XVII, No. 3 


no 
or 
to 


Since Z(v + c) is a class k function of v, Theorem 3 may be invoked to complete the 
proof. It should be noted that ¢ need not be a positive number if Z(s) is class k. 

An application of [1, Theorem 10] to the quantity W(t)/e' yields still another re- 
striction on the unit impulse response of any transfer function. 


Theorem 8. Let | W(t)e~°* | be less than or equal to the positive constant M for t > r. 
For any transfer function Z(s), let c be such that the tmaginary part of Z,-2(¢ + jw) is 
non-positive for w > 0. Then, forO < y <1, 

: | K(yr) sin ry , 2My‘sinzry < 1 | 
Wiy7r) | < e&” — , i px | 
Y =" | mk — 1)! T 1 p> p (p —y) 
Appendix. As an illustration, the bounds on the unit impulse response of the trans- 


fer function, 


Zs) = ate a 2 ' 
(s + p) 


will be obtained. The only restriction imposed upon the real numbers u and p is that p 
is positive. The unit impulse response corresponding to Z(s) is 


W(t) = onl +E . e e|. (23) 





Now, the lowest real number c which makes Z(s + c) a class 2 function will be found. 
According to the definition, Z(s + c) will be class 2 when Z,(s + c) is positive real. But, 


e+ H#te 

7 _ ° _So + Ke —_— gee 2 - 

2,(8) = ; das, = a 
2 (Sp + p) (s + p) 


Calculating the real part of Z,(s + c), it can be seen that this real part is non-negative 
only when ¢ satisfies the following inequalities 


: 7 es 
c2- 9 


lA 
> 


for MK 


u — 3p 
c> 9 


IV 
> 


for pz 
Thus, by Theorem 7 

W(t) | < t exp (- c + ‘ for wp 
and 

| W(t) | < t exp (! ad 7 for p> p. 


By comparison with expression (23), it is easy to show that these inequalities are indeed 


true. 
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BOOK REVIEWS 


Lineare algebra. By Werner Graeub. Springer-Verlag, Berlin, Gottingen, Heidelberg, 
1958. x + 219 pp. $9.35. 


A remarkably complete text on the subject of the title. By linear algebra is meant the theory of 
finite dimensional vector spaces and the linear and multilinear functions on them. This theory is con- 
veniently subdivided into three parts, the general theory which applies to vector spaces over an arbitrary 
field, the theory of real vector spaces, and the theory of complex vector spaces. The principal results 
in all three branches are presented here. The book is divided into eleven chapters. The first two on 
“Lineare Riume’”’ and “‘Lineare Abbildungen und Gleichungssysteme”’ are prerequisites for all the chap- 
ters which follow. Most of these latter, however, can be read independently of each other. Chapter IIT 
“Determinanten”’ is needed in order to read Chapter IV “Orientierte Lineare Raume’”’ which deals with 
real vector spaces and is perhaps the most unusual part of the book. The main result of the chapter 
states that the basis z, can be deformed continuously into the basis y, if, and only if, the determinant 
of the mapping z, — y, is positive. Chapter V “Multilineare Algebra’’ is the longest in the book and 
treats the algebraic theory of tensors in an invarient (coordinate free) manner a la Bourbaki. The next 
four chapters are concerned with real vector spaces. Chapter VI “Der Euklidische Raum’’ introduces 
the scalar product. Chapter VII “Lineare Abbildung Euklidischer Raume’’ is concerned with Eigen 
values and the spectral theorem for self adjoint transformations. Chapter VIII ‘“Symmetrische Bilinear- 
funktionen”’ treats quadratic forms and the inertia index. Chapter IX “Flachen zweiter Ordnung”’ 
presents both the affine and Euclidean classification of quadric surfaces. Chapter X “‘Unitare Riume”’ 
deals with the spectral theorem for Hermitian and normal transformations. The final chapter XI “Invari- 
ante Unterraume”’ returns to the general theory giving an invarient treatment of canonical forms of 
transformation and the Cayley-Hamilton Theorem. 

The book is suitable for a text, probably on the graduate level, and has many exercises in each 
chapter. I suspect it will be even more useful as a reference book in the subject for practicing mathe- 


maticians who are not specialists in this particular field. 
Davin GALE 


Readings in linear programming. By 8S. Vajda. John Wiley & Sons, New York, 1958. 

vii + 99 pp. $3.00. 

In this introduction to linear programming typical applications are illustrated by numerical examples 
which suggest the appropriate concepts and techniques. Only very elementary algebra is required for 
the understanding of the book, and the exposition is clear and easy to follow. The contents of the book 
are indicated by the following (incomplete) list of chapter headings: Transportation Problem, Caterer 
Problem, Production Scheduling, Transhipment, Bid Evaluation, Flow through a Network, Ship Sche- 
duling, Personnel Assignment, Routing Aircraft, Investment, The Simplex Tableau, Nutrition Problem, 
Airlift, Blending of Aviation Gasolines, Smooth Patterns of Production, Duality, Selection of Products, 


Train Loss Reduction, Attendant’s Rota, Warehousing, Games, Bibliography. 
W. PRAGER 


Calculus of variations and its applications. L. M. Graves (Editor). McGraw-Hill Book 
Company, Inc., New York, 1958. v + 153 pp. $7.50. 


This important book contains the papers which were presented at the Eighth Symposium in Applied 
Mathematics sponsored by the American Mathematical Society and the Office of Ordnance Research, 
and held in April, 1956. The papers are noteworthy for their clarity of exposition and avoidance of 
an excessively cramped style, so that they should be of interest to physicists, engineers and operations 
analysts, as well as to mathematicians. They cover a wide range of topics in the calculus of variations 
and its applications to elasticity, plasticity, electromagnetic theory, mathematical economics and 
other fields, as the following list of authors and titles shows: 


(Continued on p. 262) 











HEAT CONDUCTION FROM A CYLINDRICAL SOURCE 
WITH INCREASING RADIUS* 


BY 
H. R. BAILEY 
The Ohio Oil Co., Littleton, Colorado 


Introduction. The theory of heat flow due to conduction from a moving heat source 
is of interest in a number of applications; for example, welding [8]**, heat exchangers 
[2, 3] and progressive freezing of a liquid [2, 3]. In this paper we consider heat conduction 
in an infinite homogeneous medium from the surface of a cylinder of finite length whose 
radius is increasing with time. This problem arises in connection with secondary oil 
recovery by an underground combustion process, e.g., [4, 7]. The restriction of the 
cylindrical source to a finite length corresponds to considering an oil reservoir of finite 
thickness and including vertical heat losses to media bounding the reservoir. 

The differential equation describing this problem is written and the Greens function 
method is applied to obtain a solution in the form of an integral. A limiting value of 
this integral is then obtained for the case of the source moving at a constant velocity 
with no vertical losses. The problem is to evaluate a limit of the form lim,.. fo f(t, r)dr 
when /(¢, 7) has an asymptotic representation, for (¢ — 7)/r sufficiently large, which 
can be integrated explicitly. For the integrand considered in Sec. 2, it is shown that the 
integral can be divided into two parts, namely 


~ 


at/N 


[ w t) dr = i f(t, 7) a+ fe f(t, 7) dr, 


where the last integral goes to zero as t + @ and the integrand in the range [0, ¢/N] 
can be replaced by its asymptotic expression and evaluated in terms of N for t > ~. 
Finally, the desired limit is obtained by passing to the limit as N —> ~. 

In See. 3 an explicit evaluation of the integral solution is obtained for the case of 
the radius of the cylindrical source increasing at a variable velocity, namely rz = V/rr . 
This result is obtained by showing that the solution in a transformed coordinate system 
is independent of time and thus the partial differential equation in the new coordinate 
system reduces to an ordinary differential equation which is solved explicitly. 

1. A heat conduction problem. The partial differential equation in cylindrical 
coordinates for the temperature in a homogeneous conducting medium is given by 


aT ,10T , &T 20T 1 
or rar az "oat &k 2,1, 0), (1) 


where r is the radius, ¢ is time, a’ is reciprocal diffusivity, k is the conductivity,  (z, r, ¢) 
is the source function and T' is the temperature which is independent of the angular 
position, 6. The source is assumed to be at the surface, r = r,-(t), of a cylinder with a 
fixed axis at r = 0 and height A as shown in figure below. 





*Received July 7, 1958. 
**Numbers in brackets refer to the References at the end of the paper. 
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The intensity of the source is given by L(¢) in units of heat per unit of surface area per 
unit time. The source function, (2, r, t), for this type of source may be expressed in 
terms of the Dirac delta function, 6(r — rp), by the equation 


&(z,r, t) = L()Bz, h/2) o(r — rr), (2) 
where B(z, h/2) = 1 for | z| < h/2, = O for | z| > h/2. 
If we assume 7(z, r, 0) = 0, then Eq. (1) has a solution of the form [see 1,7], 


t a © alr 


T(z,r, t) = (4xk)™ | dt, | | [ B(2 5 To , bo) GTo Aro AO, d2o , (3) 


70 Jro=0 « 


where @ is given by Eq. (2) and G is the Green’s function corresponding to the left side 
of Eq. (1). G is given by the formula [see 1, 7], 


G = Gle,7r, 0, t | 20 , To » Oo ; bo) 
= 2°ar'(t — t)*” exp {—a’[R® + (2 — 29)")/4(t — t)}, 
where R? = r’ + r? — 2rr,cos(@ — 0). 


After performing the indicated integrations in Eq. (3) with respect to z) , % and ro 
and making the substitution ¢ — ¢, = rt we obtain 


T(z,r, ) = (4k)”’ | drr ‘Le exp |= - = re) |, (2) 
Jo 7 or 


.a(h/2 z) -a(h/2 — 2) 
|e avn + Ps ont avn i 
2Vvr 2V rT 


where LZ and rr are evaluated at ¢ — 17, J,(z) is a modified Bessel function defined in 
Sec. 2.1 and erf x = 24r'” f5 exp (— a)da. 

A reasonable assumption in an underground combustion process is that L = qr; 
where q is a constant. Two cases of this problem are of particular interest in underground 
combustion, namely: Case I=constant radial velocity, rr = Ut, and Case II1=frontal 
radius given by r; = 2V¢t. The solutions for these two cases may be obtained from the 
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above equation; and for no vertical losses, i.e. h — ©, are given by Eqs. (4) and (5) 


respect ively. 


. t ihe T2744 2 T(t — 


T(r. Py ee (Qh) 'qV dr - exp (a + (= oy fereve— ory, (5) 


In Sec. 2 a steady state solution of (4) is obtained, that is, the lim,... T(r, t), assuming 
that s = r — Ut remains finite as t ~ o. Thus a steady state solution for r positions 
near the source is determined; this has been called quasi-stationary state [see 6]. In 
Sec. 3 an explicit evaluation of the solution for Case IT is obtained. 

2. A steady state solution for the constant velocity case. 2.1. Preliminary Consider- 
ations. In this section we prove a lemma, state some known properties of J,(z), and 
discuss an explicit evaluation of an integral. These results will be needed in Sec. 2.2. 

Lemma. Let f(t) be a real function of the real variable / satisfying the conditions 
(a) 0 < f,(N) < lim,... f(t) < 7,(N) for all finite N, and 

(8) lim f,(N) = lim f.(N) = P, then lim f(t) = P. 
Ws. “—— mee 

Proof. By (8) there exists an N, such that | f,(N) — P| < € and! f.(N) -— P| <.e 
for all N > N, . Combining this with (a) we have P — ¢ < lim,.. f(t) < P + € and 
by choosing NV, sufficiently large we can make ¢ small and the lemma is proved. 


The modified Bessel function, /,(z), is defined by J,(z) = (2r)7* fi e*°°""d0. I,(z), 


for real z, can be represented for large values of z by its asymptotic series 
T(z) ~ (2nz)7'? e* [1 + O(1/z)], and thus 

lim (2rz)'’’e"*J,(z) = 1. (6) 
Since e *J,(z) is a continuous function of z and lim,.... e *J)(z) = 0, then e~*J,(z) is 


bounded for z > 2» for any 2p . 
We define a function g(z) by the equation 


e ‘I,(z) = (2nz)'’*[1 + g(2)], 2 = 0. (7) 
And if we define g(0) = — 1, then g(z) is a continuous function of z. Finally from (6) 
and (7) we have that lim,... g(z) = 0. 
An explicit evaluation for the integral, 


v(u,v, t) = [ x”? exp (—u’/x) exp (—v’x) dx, u #0, 


has been given by W. Horenstein [5]. The result is 
Yiu, v, t) = (Qv)'e'?[—2 sinh (2 | uw |v) +e?” erf Wt’? + [ul t'”) 


+e?" erf vt’? — Jul t'”)). ” 


The result of Horenstein does not include the absolute values signs on the u’s and 
these must be added to make the formula correct for u < 0, thus we must have (wu, v, t) = 
¥(— u, v, t). In the application of this formula in Sec. 2.2, the case u < 0 corresponds 
to positions inside the cylinder, i.e. s < 0. Passing to the limit as > @ in (8) we obtain 
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lim Y(u,v, t) = va’ e7?"™"", (9) 
t-@ 
2.2. Evaluation of a limit. Equation (4) may be written in the form 
T = (2%)'qU'l, (10) 
where 


- 


[= | dr(t — r)r' exp {-—a’[r — U(t — 7)]?/4r} exp [—a’rU(t — 7)/27] 


-T,farlU(t — 7)/2r] 
and putting s = r — Ut, U > 0, we obtain 
J = | dz(t — r)7r' exp [—a’°(s + U7r)’/4r] exp [—2(t, 7) Jo[e(t, 1], 11) 
where 
2(t, r) = (27) a U(s + Ud(t — 7). 12) 


We shall temporarily assume s ~ 0 and in this case the integrand in (11) is bounded 
in 0 <7 < ¢, since for 7 # O the integrand is a product of continuous functions of r 


and the integrand approaches zero as tT > O 


I may be divided into two parts, ‘ fg 7 qe. where J, So” f(t, r)dr and I Pb y 
f(t, 7)dr with f(t, 7) written for the integrand in (11), where N is chosen > 1. For 0 < t/N 
<7< it, we haveO < (t — 7r)/r < (t — t/N)/(t/N N — 1. Then it follows from 


z(t, r) is bounded below for0 < tN <7<t 


Eq. (12), with U > 0 and s finite that 
*“""T, [2(t, 7)] is bounded. 


and N > 1, say z(t, 7) > 2 and thus, by See. 2.1, that « 


Hence there exists an JJ such that 0 < « ky : lz 5. T)] < M for 0 : t/N < 
< t, N > 1. The left side of the above inequality is true since [)(z) > 1 for all z. Using 


7 


the above evaluations we can now evaluate /, . J, can be written in the form 


i. = | drr (t—T exp (—a's'/4r exp (— asl7/2 
-exp (—a°U~r/4) exp [—2(t, 7) ]Jofz(t, 7)] 


and thus 


0< /, < (NV — 1)M exp (—a’sl’/2) | exp (—a°l~ 1/4) dr 
V — 1)M exp (—a’sl’/2)(—4)(aU’ exp aU" t/4) — exp (—a@ U't/4N)] 
and thus for t — @ and for all finite N we have 
lim J.(t, N) = 0. 13) 


If we replace « [,{z(t, r)] in J, by its equivalent given in Eq. (7) we obtain 


I, = | dr(t — r)r exp [—a” s+ Ur)'/4r]}1 + glz(t, 7)]} [2xe(t, 7)] 
: F ig iia 14) 
= fa Unr(s + UD] drt — 7) ’"r ** 


-exp [—a°(s + Ur)*/4r] {1 + glz(t, 7], 
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where z(t, 7) in the last factor of the integrand of the first equation of (14) has been 
replaced by its expression in Eq. (12). If we define 
h(t, 7) = (t — 7)'?r"'? exp [—a’°(s + Us)’/47], h(t,0) = 0 and gfz(t, 0)] = 0, 
then the integrand of the second equation of (14) can be written as the product of 
h(t, r) and 1 + g [z(t, r)], where both factors are continuous functions of 7; and A(t, 7) > 0 
for 0 < r+ < t¢/N and for all t > 0. Hence we may apply the first mean value theorem for 
integrals to the integral in the second equation of (14) and obtain 
t/N 

I, = la’ Urls + UO)? {1 + gle(t, &)} [ dr(t — 7)'*r"'” exp [—a°(s + Ur)’/47], 
where 0 < ¢ < t/N. For 0 < +r < t/N we have (t)'?(1 — 1/N)'” < (t — 7r)'7 < (®)'” 
and thus 

I, < OU? [a Un(s + UOS'? EL + gle(t, &]} I(t, N) (15) 

I, > ¢?°U — 1/N)'?[e?Ur(s + UOT? {1 + gle(t, Dt, N), 


where 


1M 
J(t, N) = [ dr r'” exp [—a’(s + Ur)’/4r]. 


“0 


For 0 < — < t/N, we have (t — £)/ & > (t — t/N)/(t/N) = N — 1; then 


> =_ Ss 


t, &) 2'a’U(s + Ud(t — &)/E > 2'a°U(s + UY(N — 1), 
and thus lim ,.. 2(t, ) = © and finally we have lim,... g[z(t, £)] = lim,.. g(z) = 0. 


If we now pass to the limit as ¢ > © in (15) we obtain 


(1 L/N)'?(aU) x '? lim J(t, N) < lim J, < (alU’)"'e'” lim J(t, N) (16) 
and evaluating lim,_. J(t, N), we have 


JN 


lim J(t, N) = lim [ dr 7” exp [—a’°(s + Uz) 


t—o v0 


4r] 


= exp (—a’slU’/2) | dr 7 '” exp (—a’s’/4r) exp (—a°U*1/4) 


and from Eq. (9) we have 


lim J(t, N) = exp (—a’sU /2)2(aU)"'(r)'” exp (—a’ | s | U/2). (17) 


Combining (13), (16) and (17), we have 


(1 — 1/N)'?2(aU)~* exp [—a?U(s + | s |)/2] < lim (7, + J.) 


< 2(al’)~* exp [—a°U(s + | 8 |)/2] 
and the above inequality is of the form discussed in the lemma of Sec. 2.1 and thus 
we may pass to the limit as N — © and obtain 


lim J = 2(aU)~* exp [—a°U(s + | s |)/2]. (18) 
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In the above calculations we have assumed s + 0; if s 0, we may replace the 
lower limit of the integral in Eq. (11) by e, repeat the calculations and then pass to the 
> 0; and we obtain (18) for the case s = 0. 
(8) and (18) we obtain fort — « 


limit as 


If we combine Eqs. 
lim T = 2k"ql? lim | (t — Ar" 
‘exp f-—a [re + Ut — 7) )/4r} Ja rU(t — 7) /2r7| dr (19) 
gh a exp — a [ (s ae Ss Zi. 

The above result is the same as the corresponding known result, [see 6], for a plane 
source moving with constant velocity. This would be expected since the eylindrical 
source approaches a plane source for large radii. 

3. An explicit solution for Case II. Substituting a rrp and r Af in Eq. (5) 
gives 

a ax —a Via +1—X a Va(1 — X) 
rae 2h) gn | exp | (20) 
a: r 2d Dn ‘I 
Thus 7, when written as a function of x, is independent of ¢. The partial differential 
equation, Iq. (1), for the corresponding case may be written in the form 
oT ‘ _ \ aT eS i 
2 +- ( ——- a Ver) : — 2a” | [= 0, 21) 
or bs OX ot 
where we have made the substitution Ti Py >. The term 0°7 /d2" = 0, since, for the 
case h T is independent of z. The source term is replaced by its equivalent con- 
dition, namely 
oT oT : 
— - —(k) repLit —qV/k, 22) 
or Ol ; 
where 
oT oT 
Ox Ox 
>Uasr- x 


are the right and left hand derivatives, respectively, at 2 
Due to symmetry no temperature gradient exists at 7 0, also T 
23) 


The corresponding boundary conditions on T(x, t) are 
Q. 


oT : ae 
— (0, 2) = 0; lim T(x, ¢ 
Ox ‘ 
solution of Eq. (21) subject to conditions (22) and (23). Since the 
0 and Eq. (21) reduces to the ordinary 


Thus Eq. (20) is 
20) is independent of ¢, then d7'/dt 


(24) 





solution 
differential equation 


aT ; oe 
X—3 + (1 + avr) == (), 
dx dx 
can be integrated giving the two solutions 


Equation (24 
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T = constant 
* dz ree 
7 = C, [ _ exp |[—a’ V2"/2]. 


To satisfy the boundary condition, d7’/dx = 0 for « = 0, we must choose the solution 


T = constant for 0 < x < 1; thus the boundary condition corresponding to (22) becomes 
dT —qV 
AL \-n14 k 


and the solution in the range x > 1 is given by 


T —qV(k)~' exp (a’V /2) | d(z)(z)~' exp (—a’ V2? /2) ° 
{ ] | ] (25) 
—(2k)~'gV exp (a°V/2) ki(—a’ V2x"/2), x>I1 
where C, and C, have been chosen so that (25) will satisfy the boundary conditions at 
x land x ©. Finally, the value of the constant temperature in the range 0 < x < 1 
is obtained by requiring that the solution be continuous at x = 1 thus 


T = —(2k)~'qV exp (aV/2)Ei(—a?V/2), O<2<1 (26) 

where Ei is the exponential integral, — Hi(—x) = {2 e* z ‘dz. Returning to the r 

coordinate, by putting r = arr(t) in the solution (25) and (26), gives the desired explicit 
solution for Case II. 

Since this solution is unique, we may equate the solution (20) to the solution (25), 


(26) and obtain the following interesting evaluation, with a°V//2 = R, 


Pp” aa —R(a + ) 2Px ~— J Fase M(—R),|al< , 
| : exp | / ; + | zal - nT = 4 ki(—k r|/sl (27) 
‘ ‘ ; ‘ \—Ei(—Rr’), |x| > 1. 


a 
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BOOK REVIEWS 


(Continued from p. 254) 


1. E. Reissner, “On Variational Principles in Elasticity” 

1, D. C. Drucker, ‘Variational Principles in the Mathematical Theory of Plasticity’’ 

3. P. G. Hodge, Jr., ‘Discussion of D. C. Drucker’s Paper ‘Varviational Principles in the Mathe- 
matical Theory of Plasticity’ ”’ 

4. J. B. Keller, “A Geometric Theory of Diffraction”’ 

5. J. B. Diaz, “Upper and Lower Bounds for Eigenvalues’”’ 

6. J. L. Synge, “Stationary Principles for Forced Vibrations in Elasticity and Electromagnetism” 

7. H. F. Weinberger, “A Variational Computation Method for Forced-Vibration Problems’’ 

M. M. Schiffer, ‘Applications of Variational Methods in the Theory of Conformal Mapping”’ 

9. R. Bellman, “Dynamic Programming and Its Application to Variational Problems in Mathe- 

matical Economics” 
10. S. Chandrasekhar, “Variational Methods in Hydrodynamics” 
11. E. H. Rothe, “Some Applications of Functional Analysis to the Calculus of Variations” 
Unfortunately, only a one-page index is provided. 

The personal tastes of the reviewer governed the choice of three of the longer papers for further 
discussion. 

Keller introduces a generalized theory of geometrical optics designed to include the effects of dif- 
fraction. This is accomplished through the introduction of new rays, diffracted rays, which account 
for the appearance of light in shadow regions and which alter the light found in illuminated regions. 
These diffracted rays arise when a ray impinges on an edge or a vertex, or grazes an interface. First 
the diffracted rays are specified explicitly by describing the rays which arise under various circumstances. 
Then the rays are characterized by an extension of Fermat’s principle, with the equivalence of the 
two descriptions following from standard variational considerations, Following this diffracted wave- 
fronts, eiconal functions, and imaginary rays (which account for the light found in the region on that 
side of a caustic through which no ordinary rays pass) are introduced. Finally, to achieve a quanti- 
tative theory, amplitude and phase functions are associated with the rays. The author states, justi- 
fiably, “if this ray theory had been available at the time of the controversy between ray and wave 
theory, it might have forestalled the acceptance of the latter.’’ The paper concludes with an interesting 
section relating this theory to earlier works in diffraction theory, along with an extensive bibliography. 

The determination of approximate values for the characteristic numbers of self-adjoint differential 
operators is of great practical importance. The Rayleigh-Ritz method furnishes upper bounds and the 
method of Weinstein lower bounds. Using the equation for the vibrating clamped plate as an example, 
Diaz states Courant’s minimax characterization of the eigenvalues and notes that by contracting the 
class of admissible functions upper bounds are obtained, and by enlarging it, lower bounds, which leads 
to a unified view of the basic ideas underlying these two methods of approximation. Then each is discussed 
in detail, along with Aronszajn’s extensions. The paper is rounded out with some remarks concerning 
estimates of the errors and relations between the eigenvalues for the vibrating plate and vibrating 
membrane problems. Though no numerical examples are given, many references to the literature are 
provided which still further enhance the value of the paper. 

The purpose of Bellman’s paper is the discussion of a variety of optimization problems arising 
in mathematical economics. The approach is via the functional equation technique of dynamic program- 
ming, a field in which the author has pioneered and which has been extensively explored in recent years. 
Interest centers on the analytical and computational determination of optimal policies for use in multi- 
stage decision processes. The key to the solution of these problems is provided by the principle of opti- 
mality: an optimal policy has the property that whatever the initial state of the system is, and what- 
ever the state is that results from the initial decision, the remaining decisions must constitute an optimal 
policy with regard to the state resulting from the initial decision. In turn, this is but a special case of 
the more general principle of invariant imbedding (Proc. Nat. Acad. Sci., v. 42 (1956), p. 629), which 
can be used in treating a number of problems in mathematical physics. Following a brief outline of 
dynamic programming and the relationships between continuous decision processes and the calculus 
of variations, several multi-stage allocation processes and smoothing processes (of both deterministic 


(Continued on p. 270) 











IMPROVING THE CONVERGENCE IN AN EXPANSION 
OF SPHEROIDAL WAVE FUNCTIONS* 


BY 
J. MEIXNER (Technische Hochschule Aachen) 


AND 
C. P. WELLS (Michigan State University) 


1. Introduction. Spheroidal configurations have proved to be of considerable value 
in various types of diffraction and radiation problems. Of these we mention only one 
example, the prolate spheroidal radiating antenna. This has been discussed recently by 
Myers [1] who studied the radiation patterns and by Wells [2] who studied the near 
field of the antenna. Other examples and references to the literature can be found in 
the books of Meixner and Schafke [3] and of Flammer [4]. 

The advantages of spheroidal models both in a mathematical and in a physical 
sense are well known and will not be discussed here. The disadvantages are the com- 
plexity of the spheroidal functions and the lack of sufficient numerical values of the 
functions. Further, the expansions in spheroidal functions which represent the near 
field, e.g. the current on the antenna, converge slowly, a fact that is also true of ex- 
pansions for other geometric models. It is this slow convergence which will concern 
us in this note and we shall show that the convergence can be improved in the sense 
that fewer terms of the expansion are needed to obtain the desired accuracy. The method 
is not new; it has been used, e.g. by Meixner and Kloepfer [5] in the problem of the 
ring shaped antenna. However in the case of the spheroidal functions some device 
such as this is particularly useful and we present it here with the idea that it may be of 
use to others. 

By improving convergence we mean the following: Let z. a,y, be an expansion of 
eigen functions y, and consider the asymptotic behaviour of a,y, for large n. Suppose 
that to orders of 1 /n, a,, ~ b,g, . Then in the expansion 


2. (A,Wn — Dngen) (1) 


the terms of order 1/n are absent and this sum can be expected to converge faster than 
the original expansion. If in addition it is possible to write the sum )> b,¢, in closed 
form, then 

b a,v, = 7 (a,W, — Dign) + f, (2) 
where f is the sum >> b,g, . The results of this procedure, based on the prolate spheroidal 
antenna, will be given in what follows. Obviously the method can be applied to other 
types of eigen function expansions. 

2. Notation and expansions. A variety of notations have been used in spheroidal 
functions most of which have been summarized in [4]. For theoretical work the notation 
used in [3] is preferable while for numerical work other notations have some advantage. 
We shall use the notation of [3] where also the relevant information concerning spheroidal 


functions can be found. 


*Received September 12, 1958. 
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The spheroidal coordinates are £, n, ¢ where é is the radial variable and 7 and ¢ are 
angular variables. The ranges of the variables are > 1, | »| < 1,0 < g < 2m. The 
angular functions for the symmetrically driven antenna are ps} (n; y°), where y = 27a/), 
a is the semi-focal length of the spheroid £ = const. and \ is the wave length. The radial 
functions are S\‘*’ (£; y) and are functions of the third kind in that they provide the 
proper wave function behaviour for large &. The ¢g-component of the magnetic field 


is given by 


where , is given by 


aN 


2((—& — 1)°°°S, (Eos y)]/b,/(2n + 1) = iwea | EE; n)(& — 1°)' "ps, (ns y°) dn 
and where & = const. represents the antenna, /, is the tangential component of the 


electric field and the prime indicates differentiation with respect to & . 

Equation (3) is our basic expansion. An indication of how slowly it converges for 
£ = £. is given in Table 4. Since the spheroidal functions are difficult at best to compute 
and existing tables so far are limited in range, some method which cuts down this com- 
puting and speeds up convergence is desirable. 

3. Comparison series. Ignoring the physical constants and the integral over the 
gap in (3), the problem is the convergence of 


VY = > psi'(n’; -y’)psi(n; y(n + 1/2)V,(), (4) 
n=1 
where 
1(4 2 1/2 ol (4 , 
Vig) = Sa Es M/E — IS. yD]. 


We now consider some comparison series formed by replacing the terms of (4) by 
their asymptotic expressions for large n. The details of obtaining these expressions are 
outlined in the Appendix. According to Lemma I in the Appendix we can choose as 


a& comparison series 


W = > Po'(n)PUa)(n + 1/2)W€) (5) 
with 
W,(€) = (& — 1)'°Q2(E)/n(n + 1)Q,(8), 


and where P) , P,' , Q, are Legendre functions. Numerically the comparison of V,, 
with W,, depends on y and the smaller the y the better the comparison. Table 1 shows 
this comparison for y = 1 and y = 2 and for two different values of ~. In this table 
only the the real part of V, is given; the imaginary part falls off rapidly as n increases. 

One concludes from Table 1 that (5) is, in fact, a very good comparison series if 
only y is not too large. The difficulty in using it lies in the fact that W cannot be summed 
in closed form. But even so there is some practical value in replacing the terms of (4) 
with those of (5) with some accuracy after, say, n = 6, especially in view of the limited 
tables of spheroidal functions. The functions Q,(&) at least for values of — near 1 as well 
as the functions P'(n) can easily be computed. 
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TABLE 1. 














—& = 1.00001 
y=1 y=2 
n Re(V,.) Wr Re(V,.) Wa 
] —37.16 — 21.923 29 396 —21.923 
2 —8.922 —~8.098 — 13.199 —8.098 
] —4.548 —4.362 —5.235 — 4.362 
4 —2.845 —2.774 —3.065 —2.774 
5 —1.979 — 1.946 — 2.058 — 1.946 
6 —1.491 —1.477 —1.521 —1.477 
7 — 1.163 } —1.141 —1.163 
¢ = 1.001 
= | y=2 
n te( Vn) W, Re(V,) | Wr 
| 
1 —6.763 —3.962 | 4.aee —3.962 
2 —1.744 —1.586 | —2.593 —1.586 
3 — 947 — 911 — 1.080 | — 911 
4 — .626 — .613 — .669 | — .613 
5 — .456 — .450 — .476 — .450 
6 — 353 — .351 — .364 — 351 
7 — ,292 — . 286 


From the comparison series W, with the help of Lemmas II and III in the appendix, 


we obtain as an approximation to (5) 
X = SOP. (n)PUn(n + 1/2)X,8)/n(n + 1), (6) 
n=1 
where 
X,,(é —} eoth u — (n + 1/2)[K,[(m + 1/2)u]/Ko[(n + 1/2)u] 
— 1 —1/(2n + lhu — n(n + 1)/(n + 1/2) — 1/(4n + 2)] 
with ¢ = cosh u. Let us put XY = X, + X, + X; + X, and 


X, = DO Po'\(n) PU) X,.8), t= 1,2,3,4 


1 


where 
X —4i(n + 1/2)/n(n + 1)] coth u, 
or —|l, 
x —1/4n(n 4 1), 


: —l(n + 1/2)? /n(n + 1))[K,(u,)/Ko(u,) — 1 — 1/(2n + Iu} 
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and u, = (n + 1/2)u. The coefficients are of the order X,, = 0 (n7"), Xi3 = 0 (n7?), 
X,, = 0 (n~*). Thus it would appear that it would be sufficient to take only X, asa 
comparison series in order to improve convergence. However this improvement becomes 
effective only after the coefficients of (6) with X,, in place of X, approximate those of 
(5) closely enough. And this is the case only for rather large values of n if £ is very close 
to 1. In order to approximate the terms for smaller 7 also it is useful to take X, + X, + 
X, + X, as the comparison series. 

As Table 2 shows X,,3 is never very essential and can be omitted. For & = 1.081, the 
main contribution is given by X,, although the additional consideration of X,, and 
X,, is useful. For £ very close to 1, X,, and X,, are of much greater importance than 


An2 . 











TABLE 2. 
é = 1.081 t = 1.02 
t = 0.4, coth u = 2.632 u = 0.2, coth u = 5.067 

n Kaa Xn2 X br = Xne Xn3 Xne 
1 —0.988 —1 —0.125 0.181 —1.900 —1 —0.125 0.495 
2 —0.548 —1 —0.042 0.077 — 1.056 —1 —0.042 0.220 
3 —0.385 —1 —0.021 0.041 — .739 | —-1 —0.021 0.122 
4 —0.296 —1 —0.013 0.025 — .570 —1 —0.013 0.081 
5 —0.241 | —-1 —0.008 0.017 — .465 —1 —0.008 0.061 
6 —0.204 —1 —0.006 0.012 —.392 | -1 —0.006 0.050 

¢ = 1.001 £ = 1.00001 

t = 0.0447, coth u = 22.38 u = 0.0047, coth u = 223.6 

n X X X Raa to + Xn3 Ras 
1 -8.40 -] —0.125 3.64 —84.0 —1 —0.125 16.5 
2 —4.67 -1 —0.042 1.76 —46.7 —1 —0.042 26.3 
3 —3.27 —1 —0.021 1.11 —32.7 —1 | -0.021 18.0 
4 —2.52 -] —0.013 0.76 — 25.2 —1 —0.013 13.8 
5 —2.05 1 —0.008 0.58 —20.5 —1 —0.008 10.8 
6 —1.73 -] —0.006 0.47 —17.3 —1 —0.006 9.0 


In the next section we shall show that X, and X, can be written in closed form. The 
series X, can be neglected in comparison with X, . The series X, is more difficult to 
handle since it cannot be summed in closed form. However as Lemma IV shows, good 
approximations exist for both large and small values of &. Some values of K,/K»y and 
approximations for large and small argument are given in Table 3. 
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TABLE 3. 








| | 
Un Ky(un)/K (un) | un [y + In up/2] | 1 + 1/2u, | 1+ 1 — Ki/Ko 


| | | Qu, 
0.02 12.4 —12.4 26 13.6 
0.04 7.5 —7.5 13.5 6.0 
0.10 $05 —4.12 6.0 1.95 
0.20 2.73 —3.0 3.50 0.77 
0.40 1.97 2.25 0.28 
1.00 1.43 1.50 0.07 
2.00 1.23 1.25 0.02 


teturning to the original expansion (4), it can be seen from the coefficients V,,(é) 
that the convergence improves as £ approaches 1. However it does not follow that the 
comparison series give the best results for — in the neighborhood of 1. If we study these 
comparison series for various values of ~ we find that X, can be neglected in comparison 
with X, and X, if (n + 1/2)u > 0.5. This holds for all n if > 1.54. For smaller values 
of &, X, plays an essential role. Thus for § = 1.00001, X, and X, give the main contri- 
bution for n < 0.5/u or n < 100, approximately. This can be seen from studying the 
coefficients X,, and X,, in comparison with X,, . Hence if we are interested in values 
of & < 1.54, the series XY, is important and since this series cannot be summed we must 
compute or approximate the sum numerically. This, however, is not difficult in view of 
the simple analytic form of the coefficients X,, and their approximations as given in 
Table 3. 

4. Closed form expressions for XY, and X,. We consider first X, and notice that 
this series is the Green’s function for Legendre’s differential equation 


(1 — 9 )y’’ — 2ny’ + [n(n + 1) — m’/(1 — Jy = 0 

with n = 0 and m = 1. If we label this Green’s function as G(n, ’) then 
Gin, n’) = YA + 0 -9)/A-n)d+ 7)]!7, 2a’ 
=I —nd+n)/1+n1—-)*%, ave” 


° . , *1/2 . . ° ° 

when expressed in terms of the solutions [(1 + )/(1 — »)] of the differential equation. 
In order to sum X, we make use of a result of Watson [6] on summation of the 

Gegenbauer functions. This result, specialized for our purposes, states that 


; sin’ w dw 
sin @sin¢ | . 


np! >-1 ona 
a3 —_7 > os 6) 2S A 
(1 — 2¢ cos Q + ry n s UP nes (CO! Fees (cos ¢) (7) 


where cos 2 = cos 6 cos¢g + sin @ sin ¢g cos w. If we make the substitution w = x — 2y, 
the integral on the left becomes 

* [" , sin 2ydy __ 

“Jo (1 — 2tcos2 + ¢)*” 
For ¢t = 1 this becomes 


(1 — k’ sin’ y)” ( 


. 2 [ 2 ; _sin® 2y dy i 8_ (2 = r)K — 2E] 
ge" /0 r ke | -_ 
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where a = 2 — 2 cos (0+ ¢), k = 2sin 6 sin g/{1 — cos (@ + ¢)] and K and £ are 
the complete elliptic integrals. The integrals fail to exist when 0 = ¢. 

Now the expansion on the right of (7) converges for | ¢ | < 1 and since P' (eos 8) 
P;' (cos ¢) is 0(1/n) for large n and for @ and ¢ between 0 and z, the expansion con- 
verges also for ¢ = 1. Then by Abel’s theorem [7] the equality in (7) holds'for ¢ = 1 


except for 6 = ¢. Hence we have 
X, = — >> P! (cos 0)P,,' (cos g) = (1/rk’) ese [(6 + ¢)/2][(2 — )K — 2E] 


with 7 = cos 0, n’ = cos¢. 

5. Some numerical results on convergence. In Table 4 we give some numerical 
results which offer some indication of how well the convergence can be improved. The 
first six terms of the series for V and for X as well as for their difference are given for 
y = 2 and for —€ = 1.001 and é = 1.00001; for y = 1 only the first five terms are given 
with € = 1.02. Such a table is limited by the amount of information available on numerical 
values of the spheroidal functions and the above values of £ were chosen since for these 
we have sufficient information to make a reasonable comparison of V and X. No attempt 
has been made here to tabulate X, and X, from their closed form expressions although 
such computations can readily be done. 


TABLE 4. 


5 io 2 a.” l 
= 1.00001 ¢ = 1.001 t = 1.02 

J X V-X J X V-X J X V-X 
17.05 —16.6 33.65 4.25 —2.54 6.79 —1.695 —1.09 — .605 
—2.15 —1.38 —0.77 — .423 — .255 —Q.168 — .136 —.121 — .015 

.916 603 313 .189 122 .067 085 063 .022 
— .538 — 602 064 —.118 | —.137 019 || —.0599 — .074 014 

.995 1.107 —.112 230 . 206 — .026 . 126 146 — .020 

.156 135 021 .037 .033 004 

For all considered values of ~ the values of V — X are fairly small after n = 6 and 


can be neglected if two or three-place accuracy is wanted. If the comparison series is 
not used many more terms would be needed to obtain the same accuracy. 
Appendix. Lemma I. For large n, the coefficients V,,(&) are given approximately by 
W,(€) = (F — 1)'°Q2(6)/n(n + 1)Q,(€). 
From reference [3] we find that for large n 
psn; y°) & P\(n), 
SEs y) ee —iQh(e)/Kily)n(n + 1). 
Then 


— 1)'7 Si (Es y)]’ 
;y)] (8) 
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Now Q'(&) (e — 1)'” dQ,/dé and from the differential equation satisfied by the 
Q, we find [(# — 1)'* Q'(é)]’ = n(n + 1)Q, . If we make these substitutions in (8) we 
get the desired result. 
Lemma ITI. For large n, 
Q),(&) ~& (w/sinh u)'?Ky((n + 1/2)u). 


with — cosh u and K,(w) the modified Bessel function of the third kind or the Hankel function 
of imaginary argument, H, (iw). This result follows immediately from the results of reference 


8}. 
Lemma III. For large n 
& 1) QUE) /Q,(E) & —(n + 1/2)K,(w)/Ky(w) + 1/2u — § coth u 
where w (n + 1/2)u and K,(w) ts again a Hankel function of imaginary argument. 


This result follows directly from Lemma IT. 


— w' fy + In w/2]"', approximately, 


Lemma lV. For small w, K,(w) K,(w) = 
= 2w — 1/8w’*, 


where ¥ 0.5772 7s Euler's constant; for large w, K,(w)/K,(w) = 1+ 1 
approximately. 

The first result follows from the definition of A,(w) and K,(w) in terms of the modified 
Bessel functions J)(w) and J,(w) together with the approximations, J,)(w) ~ 1, J,(w) + 
w/2. The second result follows from the asymptotic expansion of K,(w) which can be 


found, e.g. in reference [9]. 
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BOOK REVIEWS 


(Continued from p. 262) 


and stochastic type) are discussed in some detail. The paper closes with the solution of a problem arising 
from consideration of a miniature mathematical model of the interdependent steel and automobile 
industries. This leads to a bottleneck problem involving the allocating of steel resources over a given 
period of time to the production of (a) additional steel itself, (b) additional steel-producing facilities, 
and (c) automobiles, where the objective is to maximize the total number of automobiles produced 
during the process. 

Even a casual perusal of this volume will impress the reader with the role of the omnipresent varia- 


tional principle in applied mathematics. 
Rospert KALaBa 


Boundary layer effects in aerodynamics. International Symposium held at National 
Physical Laboratory, England, in 1955. Philosophical Library, New York, 1957. 
$12.00. 


This volume contains a complete record of an International Symposium held at the National Physi- 
cal Laboratory, England, in the Spring of 1955. Nine papers were presented, covering most aspects 
of the subject of current interest. The Symposium was opened by Professor L. Howarth, with a broad 
and enlightening survey of boundary layer theory as it stood in 1955. 

The first two papers are essentially of a theoretical character. A paper on three dimensional boundary 
layers by Timman analyses the flow on a general surface in terms of intrinsic coordinates with application 
to thin wings and yawed elliptic cylinders. Lighthill and Glauert investigate laminar boundary layer 
flow on a long thin cylinder and obtain solutions valid near the nose and a long way downstream. In 
a paper on the stability of boundary layer flow on a rotating disk experimental work by Gregory and 
Walker is compared with theoretical analysis by Stuart. In many respects the two approaches show 
good agreement and discrepancies which do arise are satisfactorily explained. 

An analysis of boundary layer transition by Schubauer and Klebanoff lends strong support to the 
turbulence spot theory of Emmons. Kitichemann considers viscosity effects on swept-back wings with 
emphasis on separation and the growth of tip vortex sheets. Pankhurst describes various methods 
of boundary layer control including the use of suction on thick wings. 

Young and Kirkby calculate the profile drag on supersonic wings of biconvex and double wedge 
sections. The symposium ends with two long papers on Shock Wave Boundary Layer Interaction. 
Holder and Gadd review the work done on this problem on the flat plate and estimate the influence 
of the effect on the calculation of base drag. The second paper by Pearcey is concerned with Turbulent 
Boundary Layers on Transonic Airfoils. 

Most of these papers have already been published separately but the collected proceedings also 
contain accounts of the discussion following each paper. This is invariably interesting and in many 


cases, quite extended. 


M. Hour 


Integral equations and their applications to certain problems in mechanics, mathematical 
physics and technology. By S. G. Mikhlin. Pergamon Press, New York, London, 
Paris, Los Angeles, 1957. xii + 338 pp. $12.50. 

This is a useful addition to the literature in English on integral equations and their applications to 
continuum mechanics. The first part of the book, which comprises approximately two-fifths of the whole, 
provides a readable account of the fundamental theory of integral equations. The classical theory due 
to Fredholm and the Hilbert-Schmidt theory for equations with symmetric kernels are treated in the 
first two chapters, and the third chapter is devoted to singular integral equations. Methods for approxi- 
mate solution of the equations are given together with numerical examples. For the most part detailed 
proofs are given for the theorems used to develop the theory. 


(Continued on p. 284) 
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APPROXIMATE DISTRIBUTIONS OF NOISE POWER MEASUREMENTS* 


BY 
WALTER FREIBERGER anno ULF GRENANDER 


Brown University 


Summary. The frequency functions of certain spectral estimates are studied analyti- 
cally and numerically. An approximation is obtained for the case of a Poisson weight 
function and compared to the true distribution. The eigenvalues of products of Toeplitz 
matrices play a crucial role in the sampling theory of quadratic forms; an approxi- 
mation to their distribution is discussed and its accuracy studied numerically. This 
leads to approximate probability densities which are thought to be valid for moderate 


or even small sample sizes. 


1. Introduction. Let x, be stationary Gaussian noise with the power spectral den- 
sity {/(A). The parameter ¢ takes on integral values only, ¢ = — 1,0,1, ---.Givena 


finite sample x, , 2 , , , we want to estimate f(A), using a quadratic form Q of the 


type 
l n 
Q=- ; R Ww, 0%, (1) 
n y,u=l 
where 
l ~ A 
wv. =— [ e wd) dd. (2) 
2nr J_s 
Here the weight function w(A) defines the spectral window and satisfies 
1) w(rA) > 0 
Bu (3) 
2) ~ | w(A) dA = 1. 
2n J_, 


If the sample size n is large, it is known that Q has an asymptotically Gaussian dis- 
tribution with parameters that can be expressed simply in terms of f(A) and w(A); see 
Grenander and Rosenblatt (1957), p. 134 [1]. On this basis a choice between various 
a priori possible estimates can be rationally made, in accordance with the well-known 
relation between the band-width and variance of a spectral estimate. 

A considerable number of particular weight functions has been suggested in the 
statistical literature. The usual procedure is to choose an even, periodic weight function 
v(A), subject to the above conditions (3), and to form 

w(r) = 4[v(A — Ao) + ofA + A)], (4) 
where A, is the frequency at which we want to determine f(A,). [f the Fourier coefficients 
of »(A) are denoted by v, then clearly w, = vy, cos kro . 


*Received October 22, 1958. The work described here has been partly supported by the U. 8. Army 
Signal Corps, under Contract DA-SC-78130. 
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to 
“I 
Li) 


We will be concerned with two particular cases. The first is the rectangular weight 


function 


aT <n 
v(A) = 

. yi SS %. 
It is of course desirable to keep the bandwidth parameter h small, but it must not be 


chosen so small that the variance becomes excessively large. 


The other is the Poisson weight function 


(A) = Se == ll a oO = p < th, (6) 
1 — 2p cos + p 
where values of p close to one correspond to narrow spectral windows. 

The Fourier coefficients v, are in the first case sin kh/kh and in the second p~'. Both 
these windows are convenient to deal with and there are reasons to believe they have 
good sampling properties. 

The distribution theory for spectral estimates is fairly complicated. It is true that 
one can write down immediately the characteristic function ¢(z) of Q as 


7 —a " 
g(z) = [] l= 2 ; (7) 


y=l n 
where the X, are the eigenvalues (always real) of the matrix product RW. Here R is 


the covariance matrix 


he 7, * De 1.2, +> 
( at ) , 
)1 a eee ; { (8) 
= i= | ee T(N) da; yp=1,2,--- ne 
|\Qr é ) 
\at J-—r / 
and 
W = {w,_, ; yp = £2, ey. (9) 


The expression (7) is of little immediate practical use since (a) the eigenvalues are in 
general very difficult to obtain numerically and (b) it is very hard to get the frequency 
function g(x) of Q by numerical Fourier inversion. Therefore, one is almost forced to 
resort to approximations, of which the most widely discussed is the Gaussian one 
mentioned above. 

tecently attention has been focused on the need for sharper approximations, valid 
for moderate sample sizes. One such approximation can be obtained by applying known 
results from the theory of Toeplitz forms. Indeed, one knows that the A, behave dis- 
tributionwise approximately as the values of the function w(A) f(A), viz. 


| 
inumber of A, < w} = 5- meas [A | f(A)w(A) < 4}; (10) 


i) ~ 


lim 


see Grenander and Szegé (1958), Chap. 8 [2]. This implies that we have approximately 


~ 2Qizr os © . Qizh(rA)wlrA 
log ¢(z) = —i p> loe( _ ) ~— | og 1 —_ aE | a. (11) 


if] 
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Before entering into a discussion of how (11) can be used to determine the approxi- 
mate frequency function of Q, one important observation will be made. It is clear that 
a spectral estimate will be useful only if it resolves the spectrum well enough so that 
w(d) is peaked around the frequency of interest \, . If the spectral window is narrow, 
measured in terms of a typical frequency of the stochastic process, then the values of 
the distributions of the two functions f(A) w(x) and f(A)) w(A) are almost identical, 
so that little is changed in (11) if we substitute the second function for the first one. 

Hence the stochastic variable Q/f(\,) will have approximately the characteristic 


function 


ar » »( 
ID exp [7zQ f(Xo) | =~ ¢,(z) = exp b 4 | log [ a am »] an}. (12) 


| @rJo. n 


This distribution apparently does not depend upon the (unknown) spectral density 
f(\), so that we can use it to construct approximate confidence limits f*(A,) and f*(A»). 
Indeed, f*(X Q/x,; with 


| g(x) dx = p, (13) 
where p is the confidence coefficient and g,() is the frequency function corresponding 
to ¢,| 

A parallel investigation, Grenander, Pollak and Slepian (1959) [3], discussed the 
approximation (11) for f(A) = 1. The approximation mentioned in the last paragraph 
does not seem to have been studied numerically in the literature. This will be done 
for a few cases in the following section. 

2. Numerical studies. Let us introduce the eigenvalues Xj (k = 1,2, --- , nm) of 
the matrix f(\,)1V; as mentioned above, these eigenvalues can be expected to have 
about the same distribution as \, . In studying this approximation numerically we 
are constrained by the limitations of our computing facilities. To avoid excessive com- 
puting effort, we have chosen n = 20. This number may not correspond to situations 
met frequently in practice, but it is hoped that the information gained will be of some 
general interest. 

The spectral density was chosen as 


f(\) = 102 + 22 cos 2d + 20 cos 4X, (14) 
a graph of which is shown in Fig. 1(1). The value of f(A) at A» = O is f(0) = 144; the 
Poisson weight function (6) was used with p = .7, .9 and .95. The eigenvalues \, and 


df are given in Table 1. 

The statistical distribution of the noise-power estimates will depend mainly upon 
the largest eigenvalues and the table shows that the error in these is of the order of a 
few per cent. We have also computed the exact frequency function g(x) and compared 
it with the frequency function corresponding to the eigenvalues Xj ; see Figs. 2-4. 

The frequency functions have been computed from the expression 


10 n(2hor) ~? 20 


q(x) = : >> (-1)"" | ee I] (1 — 2yr,)7'? dy (15) 
T > a 


y= / (2hon7—1)~ pel 


and the corresponding formula with the A, replaced by the \/ . This expression was 
given by Slepian (1958) [4]. To eliminate singularities at the ends of the intervals of 
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TABLE 1 











k p= .95 | p= .9 p= .7 
Ny ri’ Ne rx’ Ar | i’ 

l 2078.9 2114.6 1596.0 1618.3 || 723 .3 729.6 
2 386.9 405.4 552.1 576.0 || 529.4 | 547.3 
3 124.6 134.1 220.2 7.6 (| 353.4 379.0 
1 57.1 64.2 108.8 122.4 232.0 | 260.9 
5 31.6 37.6 62.1 73.9 155.2 | 185.0 
6 19.6 | 24.8 39.0 49.5 107.0 | 136.4 
7 13.0 17.7 26.2 35.7 76.3 104.3 
8 9.2 13.4 18.6 | 27.2 56.2 | 82.7 
9 6.8 10.7 | 13.8 21.6 42.8 | 67.5 
10 5.2 88 | 10.7 | 17.7 33.9 | 56.6 
11 13 7.3 «| 8.7 | 15.0 27.9 48.5 
12 3.6 6.3 T- 13.0 23.9 42.5 
13 £3 5.6 6.5 11.5 a. «(Cl 37.9 
14 2.9 5.0 || 5.9 10.4 || 19.5 | 34.3 
15 2.7 4.6 || 5.6 9.5 18.5 31.5 
16 2.6 4.3 | 5.4 8.9 18.0 | 29.5 
17 2.6 4.0 5.3 8.3 17.7 27.9 
18 2.6 3.9 5.3 | 8.1 17.6 26.8 
19 2.6 3.7 5.3 | 7.8 || 17.6 26.1 
20 2.6 3.7 5.3 | 7.6 6 5.3 





| 


integration, a change of variable 
y =c,cos7mx + d,, 


with 


; 1 
C, = ¥(Yorss — Yoo)s d, = 3 (Yous + Yor), Paes Qn, 





is effected and the resulting integrals are evaluated by a Newton-Coates five-point 


formula. 


g(x):10 * 


ad —-RW, p=.9 


, 











) “a 2 3 4 ‘ 6 7 5 9 10 x-10 
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-4 RW 
g(x) -10 Pa ae 


x-l0 





Fia. 4 


Because of the proximity of the eigenvalues, which is characteristic for matrices 
of this type, the classical iteration procedure broke down. Instead, the Jacobi diagonali- 
sation method was modified to apply to (non-symmetric) matrices which are products 
of two symmetric matrices; see Kalker (1958) [5]. The computations were carried out 
on an IBM 650 computer and took of the order of a few hours for the determination 
of the eigenvalues and of each frequency function. 

It is obvious from the graphs given below that the approximation works well for the 
narrow windows p = .95, .90, .70. If the weight function is flat one will not expect the 
approximation to be adequate and in the cases p = .5 and .3 we give only the X{ in Table 
2 and the corresponding frequency functions in Figs. 5 and 6. This table may be of 
interest to a reader who wishes to investigate wide spectral windows. 

TABLE 2 


} p= oO p= 3 
he’ dA 

1 416.4 264.1 
2 375.5 254.6 
3 322.6 240.3 
} 269.1 223.1 
5 221.9 204 .6 
6 183.3 186.2 
7 152.6 168.9 
8 128.7 153.4 
9 110.2 139.4 
10 95.6 127.4 
11 84.2 ii7 2 
12 75.3 108.4 
13 68.3 101.1 
14 62.8 95.0 
15 58.3 90.0 
16 54.9 86.0 
17 52.3 82.8 
18 50.4 80.5 
19 49.0 78.8 


bo 
“J 
ay 
_— 
. 


20 48 . 3 
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S/- 
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°F 
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2 | 
—Toeplitz approximation (26) 
ig 
O | 1 —— ~ 1 ; Se a ae 
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Fie. 5 
g(x) io 4 
5 L 
41 
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2 e 


— Toeplitz 
opproximation (2@ 





_ 4 i = i 
0 2 3 4 5 6 7 8 9 x-10° 


Fic. 6 


For the rectangular window (5) the analogous quantities have been computed for 
the spectral density 


f(\) = 19 + 6 cosd — 10 cos 2A + 2 cos 3A + 8 cos 4A (16) 
the graph of which is shown in Fig. 1(2). The value of f(A) at the center A, = 3/10 x of 
the windows is f(A)) = 17.2. The weight function (5) was used with 

he see Le * oe bad . 
‘=40° 20° 340° 8 4 


For the first three of these values we have Table 3. 
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TABLE 3 








k h = 7/40 h = x/20 h = 34/40 
Ax NK Xx r;’ Xe r;’ 

1 176.2 161.2 ] 147.5 135.3 118.9 107.2 
2 166.6 160.9 140.9 134.7 153.7 106.7 
3 13.8 10.9 14.0 35.2 | 68.2 56.5 
1 11.6 10.7 37.2 | 34.9 58.5 56.2 
5 175 1 2.7 1.9 i2.3 8.7 
6 126 1 2.0 Ly 9.0 8.0 
7 0 0 m | 0 6 0 

8 0 0 0 0 3 0 
9 | | 0 0 0 0 
10 | 0 0 
11 
12 ¥ ¥ 





While the approximation still seems adequate it is interesting to note that the relative 
errors of the largest eigenvalues are now somewhat larger than for the Poisson weight 
function. The explanation for this may be found in the local behavior of the spectral 
densities at the frequencies of interest. 

The resulting frequency functions are represented in Figs. 7-9. We also give the values 
of \{ for h = w/8 and w/4, in Table 4. 











I h = 7/8 | h = 2/4 
i’ Ak 
1 46.0 30.9 
2 45.7 26.6 
3 | 29.7 23.7 
4 27.8 21.1 
5 | 17.0 19.6 
6 16.0 18.9 
7 14.9 18.4 
8 14.9 18.4 
9 14.9 18.2 
10 14.9 18.2 
11 14.9 16.6 
12 14.9 16.6 
13 14.8 16.5 
14 | 14.6 16.4 
15 14.4 16.2 
16 14.3 15.8 
17 11.6 14.3 
18 10.6 10.4 
19 | 4.8 3.9 
20 3.3 3.5 
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Similar observations to those made in the case of the Poisson weight function seem 


to apply here. 
3. Some closed-form approximations. We have seen in the last section that if the 
spectral window is sufficiently narrow, little accuracy is lost if we regard f(A) constant 
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throughout the interval (— 7, 7). If this is done, and if f(A) is normed to 1, we can use 
Eq. (12) to give us the approximate expression for the characteristic function 


Ga(Z) = exp |-8 ve) | 


“ 2izw(X) 
¥y(z) = = | log E _ Bie) | dx. (17) 
ra n 


70 


where 


For the case of a rectangular window this equation gives us the characteristic function 
of the well-known type III distribution with 


y\ h | Qriz 
= og —_ 
” is nh 
(2) I Qriz |~"’”** 
9.(2) = -_-— P 
nh 


This corresponds to the frequency function 


so that 


PN idl ] oe : 
Ja(t) = & x" exp [—(anh/2z)], 18) 
ee" Tha 


\ Qn) 


~ t 
an approximation which may in essence be contained in the work of S. O. Rice. 

The Toeplitz approximation for a Poisson spectral window can also be expressed 
in closed form. Indeed, the function 


SS oe 2iz Il—p ' 
exp ¥(z) = exp log| 1 —- — P dx (19) 
Y or 5S n | P h 12 
Jo ; 


~ 


is analytic in z in the plane cut from 


n1il— ni + ‘ +: 
Z —- ei to z = — i, see Fig. 10. 
21+ p 21-—p 


Noticing that 


= , oe, 
Re Y(z) | < 3 on | 1 + i (! — | : (20) 
eae ne n l—p 


so that | ¢,(z) | is dominated by a power of z, we can change that contour of integration 
¢ A 5 


in the Fourier inversion formula 


g(t) = = | e '**y(z) dz (21) 
2r . 
to give us 
l i z 99 
Ja(t) =; e€ p(z) dz (22) 
oR 


However 
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Fic. 10. z-plane. 
= — w we get 


rhe last term does not contribute anything to the integral and putting z 


[tog (B | 1 — Re® |?) dd, 


] 
on ~ 





where 
B\1— Re” |? =|1— pe”)? —=(1 — 9p’) 
which determines the constants 
p 
B= = 
7 R 
and 
Qn ; ; v are (23) 
+p-— a — pe tae — [14 0° - a- |} 
n ? 
R — — ees aaa Caen 
2p 
Introducing the angle @ by 
; 2v 
[ee ele oF 
R cos @— isind@=e", cos 6 = ReR = ———— > 
2p 
we get (z) log B so that 
: R n/2 e in6/2 
¢g.(2 =B"’? = (2 = (24) 
p p 
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Inserting (24) into (22) we get 


t => ame 
IAL) = -; = $ e” dv, 
2rp” J 
where we have supposed the sample size n to be even, n = 2m. The above expression 
reduces to 
2m tl i 
9.(t) = 9 ; | € sin m@sin 6 dé 
mp (1— p) 
2me “ Ea ee . 
= — or | ee * sin m@sin 6 dé 
xp (1l— pep) Jo 
with 
l+p 
a m5 =; 
r (25) 
2mp 
B= ——, 
loot 
Hence 
4 me i = \ “ 
g(x) = +7; —3~ | e” [cos (m — 1)6 — cos (m + 1)6] dé 
wp (1— 9p) Jo 
me %74"*" 
= = [Jn -1 (182) + J 1 (882)] 
p (1— op) 
me gr" 
= —-—7, x {t”~—,,.-: (82) + o”""T.+1(82)] aot 
p (l — p) (26) 
= me “* 
= (—1)* ———;-———_ > ([I...:(82) — I,.-:(82)] 
p (l— op) 
2m*e °7(—1)""" 1 
So ———— = J (pz) 
p (l—p)é x 


- m 
= fd T,,(B2). 
Xp 

As a numerical test of the accuracy of this Toeplitz approximation we have plotted 
(26) and the corresponding frequency functions of the exact distributions for p = .5 
and p = .3 in Figs. 5 and 6. Because of the relatively small sample size used in the com- 
putations, the approximation (26) was thought to be adequate only for small values 
p (wide spectral windows); it turned out, however, to be surprisingly close, as seen in the 
figures, so that its validity may well extend to somewhat narrower windows than expected. 

4. Discussion of the approximations. We have examined the behavior of the eigen- 
values of the matrix product RW and compared it to that of f(A.) W. Because of the 
limited scope of this numerical study, our conclusions are necessarily of a tentative 
nature. It seems that the approximation is better the larger the sample size, the narrower 
the spectral window and the flatter the spectral density. Obviously, the spectral window 
has to be taken narrow enough not to smooth out too much of the large scale structure 


1959] DISTRIBUTIONS OF NOISE POWER MEASUREMENTS 283 


of the spectrum. If the time constant corresponding to this spectral window is only a 
fraction of nm we would expect this approximation to be satisfactory. 

The Toeplitz approximation (26) seems to be valid if the time constant p/l — p 
is small compared to n. 

Acknowledgment. We would like to thank Miss Cornelia M. Kalkman for her 
assistance with the numerical work. 
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BOOK REVIEWS 


(Continued from p. 270) 


The second part of the book contains a variety of applications, almost all of which are two-dimen- 
sis is on the method 





sional problems. As is to be expected and desired in a work of this nature, the emphz 
of application of integral equations to obtain the solution, so that the working is carried only to the 
point where detailed results can be found by straightforward calculations. Discussions are included 
of the integral equation approach to the solution of the wave equation and the heat conduction equation 
but most of the examples are taken from Russian work in the plane theory of elasticity and hydrody- 
namics, 

The translator deserves thanks for a translation which very seldom jars. However, the realer may 
be somewhat startled to find that Hooke has lost his final “e’’ due to the strain of being re-transliterated 
into English, and one may seek in vain for references numbered higher than thirty-five. 


R. T. SHreip 


Vibration and impact. By Ralph Burton. Addison-Wesley Publishing Co., Reading, 


Mass. x + 310 pp. $8.50. 


The author states that this book is intended for senior undergraduate, and for beginning graduate 
students, and the objectives of the text are listed as treatments of (1) “the natural frequency” (free 
vibrations without damping), (2) ‘“forees which tend to suppress vibrations’’ (free vibrations with 
damping), (3) ‘‘periodic forced and transient forced vibrations’, and (4) “self-excited vibrations’’, 
There’ are chapters on free vibration, vibratory systems commonly found in machinery, 


impact, nonlinear vibrations, measuring instruments and analogs, numerical computation of natural 


damping, 


frequencies (of multi-modal systems), waves, vibrating beams and related subjects, (the related subjects 
including such items as plate and membrane vibrations, and damping in a system having two lumped 
masses, two springs and one viscous damper), and chapters on the analysis of control systems, and 
fatigue. Nearly every chapter contains some examples and is followed by a number of problems 

In view of the title of the text, the reader may expect to find a book devoted in equal parts to vi- 
brations and to impact; he will be disappointed. There is one chapter of approximately twenty pages 
dealing with elementary aspects of the response of rigid bodies, elastically supported, to non-periodic, 
time-dependent forces. In another chapter, there is a short qualitative discussion of wave propagation, 
wave interaction, and wave reflection from discontinuities, consisting of largely verbal results deduced 
from heuristic arguments; however, this latter material does not add greatly to the study of impact 
on elastic bodies. 

The instructor who intends to use this book as a text in senior undergraduate, or in graduate, 
courses must expect students who have no background whatever in ordinary differential equations 
but a fair preparation in partial differential equations, for there are quite lengthy derivations of the 
solutions of linear, homogeneous equations with constant coefficients, while the solution to the one- 
dimensional wave equation (the only partial differential equation in the book) is given without any 
derivation whatever. 

To the reviewer, the book appears disconnected, badly arranged and replete with errors of commission 
and omission. To list some of the former, the author says of the equation —kr = mx that “unfortunately, 
we cannot solve this equation by direct integration’’. The “perturbation method” is illustrated by the 
example of the pendulum, and the author believes that this method consists of annihilating the non- 
linear terms by amputation; thus the “perturbation solution’”’ of 6 + (g/l) sin @ = 0 is obtained by 
6 = 0. The author notes that ‘Fourier analysis cannot be applied to the solution 


solving 6 + (g/l) 
.. Any solution based upon superposition 


of nonlinear systems under periodic nonharmonic forcing. 
’’. The criterion for stability of vibrating systems is said to be the same as 
that for static stability—i.e., that the system “‘... should return to its initial configuration after a 
small disturbance’’. Of Raleigh’s method, it is said that “the more boundary conditions the assumed 
solution meets, the better will be the estimated frequency’’. It is also said that “It can be shown that 


of solutions is invalid. . 


(Continued on p. 298) 
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FORCED HEAT CONVECTION IN LAMINAR FLOW 
THROUGH RECTANGULAR DUCTS* 


BY 
C. R. DENNIS, A. McD. MERCER anp G. POOTS 
The Queen’s University of Belfast 


Introduction. In this paper we consider the problem of finding the heat transfer to 
the wall of a duct which has a rectangular cross-section and through which a hot viscous 
fluid passes in steady established laminar motion. We shall make the usual simplifying 
assumptions that the thermal properties of the fluid are independent of temperature, 
that liquids are incompressible and that gases obey the perfect gas law. The first is 
strictly true only for small heat input and, of course, the assumption of established 
motion ignores the hydrodynamical boundary layer in the inlet. The problem is of 
engineering interest since, in many applications of gas-flow heat exchangers, flow passages 
are used which have small cross-section and a high ratio of surface area to core volume, 
so that the Reynolds number is often small enough for laminar flow to exist. Practical 
cross-sections can often be approximated by simple geometrical shapes and theoretical 
correlations of heat transfer with cross-section are of value in reducing the amount of 
practical test data required. 

The rectangular cross-section gives rise to an essentially three-dimensional temper- 
ature distribution and has therefore received less attention than those involving two- 
dimensional! distributions, such as the circle and the case of infinite parallel walls. Some 
results have been given by Clark and Kays (1953) [1] by considering conditions far 
enough from the thermal inlet to assume a fully developed temperature profile, but 
this gives only asymptotic results and no information is obtained on the variation of 
heat transfer in the thermal entry region. On the other hand experimental data are 
given regarding this variation and it is therefore of interest to obtain theoretical results 
taking into account the undeveloped temperature profile. This is the object of this 
paper although it is hoped that the numerical analysis of the governing partial differential 
equation, which occurs in wider fields, will also be of interest. 

Governing equations and basic thermal quantities. We consider a duct whose axis 
is the ¢-axis of rectangular coordinates (£, 7, ¢) and whose cross-section perimeter is, 
in general, the curve C(é, n) = 0. The constant cross-sectional area is A and the length 
of the perimeter is S. In customary notation the velocity field is (u, v, w), but for steady 
established laminar motion under a constant pressure gradient P/L we have u = v = 0 
and w = w(é, 7) where 

ow, ew P 
——  ——— = —— (1) 
OF On pL ’ 


= 0 on C. The energy equation governing the temperature 7'(é, 7, £) of the 


with w 
fluid is, subject to the stated assumptions, 
fs Pe 4 BE r) aT 
1 —|/=w—, 2) 
ee ary ( 
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where «x is the thermometric conductivity, and second order terms, such as that due to 
internal heat generation, are neglected. The origin is at the thermal inlet and we suppose 
that the fluid enters ¢ > 0 with constant temperature 7, . We introduce dimensionless 
coordinates x = §£/d, y = n/d and z = {/d Pé where d = 4A/S and Pé is the Péclet 
number, equal to the product of the Reynolds and Prandtl numbers, that is, dwy/k. 
The Reynolds number dwop/u is based on the mean velocity w, , that is, the ratio of 
total flow to cross-sectional area. It now appears that the ratio of axial conduction 
x07 /d¢* to the axial convection term wdT7'/d¢ is of order (1/Pé)’, so that axial con- 
duction may be neglected for large enough Pé. This can lead to discrepancies in the 
special case of low Reynolds number flow of low Prandtl number fluids, such as certain 
liquid metals, but for water, air and high Prandtl number oils it is justified. Equations 
(1) and (2) then become 


FP 


Viw= - (3) 
eI ) 
and 
; woe 
Vid = ae t) 
W,O0z2 
where 
a a T —T, 
Vi=-3 —; and @=- — 
arta y = a 


T, (< T,) being a representative temperature associated with the duct wall in the 
region ¢ > 0. The boundary condition for Eq. (3) is that w = 0 on the transformed 
boundary C’(zx, y) = 0 while that for Eq. (4) depends upon the assumptions that are 
made. If 7’, is taken as the wall temperature, assumed constant, we have 6 = | within 
C’ when z = 0 and 6 = 0 on C’ when z > O. There is also another boundary condition 
in which we interpret 7’, as the temperature of the medium just outside the duct wall, 
again assumed constant. It has been shown by Hampton [2] that, when heat losses by 
radiation and natural convection take place from a body at temperature 7’ into surround- 
ings at temperature 7’, , the flux of heat H(cal/cm*/sec) is well represented for tempera- 
tures from 0 — 400°C. by the empirical formula 


H = A(T —7,)+ BIT — T,)°. (5) 


The constants A and B depend upon the emissivity F of the body, A varying from 
1.96 X 10°* when EF = 1 to 1.33 XK 10°* when FE = O and B varying correspondingly 
from 1.71 X 10°° to 0.25 X 10°°. If we identify 7 with the temperature of the duct 
wall (assumed ideally to be of negligible thickness so that 7 is the temperature of the 
fluid in contact with it) the heat flux to the wall from the fluid is H = — kdT/dv, where 
k is the thermal conductivity of the fluid and » is the outward drawn normal from the 
duct wall. Substituting in Eq. (5) and introducing dimensionless quantities we obtain 
pe? Ao+ BIT, — 7,6, (6) 


dav’ 
where v’ = v/d. Now @ < 1, tending to zero for large z, while B is of order 10 ° A so 


that for small 7’, — 7, (which is implied in the basic assumptions) we may neglect the 





1959] HEAT CONVECTION THROUGH RECTANGULAR DUCTS 287 
second term on the right hand side of Eq. (6). Putting N = Ad/k, the complete boundary 
conditions for @ in this case may therefore be stated as 

6 = 1 within C’ when z = 0, 00/dv’ = -- N@ on C’ when z > 0. (7) 


If N is infinite we get 7 = 7, at the wall, so that this special case gives the constant wall 
temperature condition. If N = 0 we have the trivial solution 7 = T, . In practice N 
must lie somewhere in between these limits. 

The solution of Eq. (4) may be written 


6= > @,0,(x, y) exp (—2,2), (n = 1,2,3---), (8) 


where 9, and i, are eigenfunctions and eigenvalues of the membrane equation 


Aw(x, ¥) 
Wo 


Viet+ 6 = 0, (9) 
subject to the boundary condition, deduced from the second of the conditions (7), that 
00/dv' = — N® on C’. The theory of this equation is well known and is dealt with, 
for example, by Courant and Hilbert [3]. Since w(z, y) is positive within C’ the eigen- 
values X,, are real and positive and the eigenfunctions form a complete orthogonal set 
satisfying the property 


l| w0,0, dx dy=0, (msn), (10) 
J* Dp’ 


where D’ is the domain bounded by C’. Each 9, has arbitrary amplitude which we 
choose, for convenience, so that 


II 


D’ 


w0; dx dy = [| w dx dy. (11) 
Putting z = 0 in Eq. (8) we have from first of the conditions (7) that 
l= > 2,0, (12) 
aot 
so that from Eqs. (10) and (11) 


a, = [| we,, dx ay/ || w dx dy. (13) 
7D’ v/D’ 


The temperature @ (2, y, 2) is therefore known to any desired accuracy once sufficient 
0, have been found. Two further thermal quantities are of interest. Experimental 
measurements are made on the basis of a mean mixed temperature of the fluid, that is, 
6 (x, y, 2) averaged with respect to the local fluid velocity over any section of the duct. 
This temperature is a function of z only and its difference between any two sections 
gives a measure of the heat transferred across the wall between them. Denoting it by 
Ty then Oy = (Ty — 7T,)/(T. — 7;) and is given by 
6y(z) = I| wé dx dy I| w dx dy, 
J! Dp’ JJp* (14) 


= ps a: exp (—A,z). 
n=1 
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The remaining quantity to be considered is the rate of heat transfer per unit area to 
the wall of the duct, defined by means of a heat transfer coefficient. If H is the heat 
flux to a given area of the duct wall, the fundamental equation defining such a coefficient, 
h, is H = hAT, where AT is a representative temperature difference. We consider two 
coefficients, of somewhat different natures, obtained by different choices of AT’. First 
taking AT = T,, — T, we obtain the local coefficient of heat transfer h, , which measures 
the local average rate of heat transfer to the duct wall as a function of longitudinal 
distance down the duct. The total heat transferred to the area between the sections at 
¢ and ¢ + dt is dq = h,Sd¢(Ty — T,) = h_Sdé6y(To T,) and if s is the distance 
measured along the perimeter of the boundary curve C in an anti-clockwise direction 


we must also have 


rae - a ° 00 
dq = —kdt | —- ds = —K(T, — T,) dt | <— ds. (15) 
c Ov “JcOv 


We equate these two and introduce the appropriate dimensionless heat transfer coefficient 
or Nusselt number, defined as h,d/k, so that we obtain for the local Nusselt number 


06 ae ; 
——; as'/ SOx , (16) 
Ov 


( 


Nu(z) = —d [ 


where s’ s/d. We now eliminate @ in terms of the 0, by Eq. (8) and apply Green’s 


theorem to Eq. (9), so that 


F @ 6. . A. ie i 
| —= ds’ = —-— I we,, dx dy. (17) 

Jc: OV Wy JJ p’ ; 
Using Eq. (13) and since {fp wdz dy = A’w,, where A’ is the dimensionless area within 

C’, we finally obtain 
' l 
Nu(z) = } A,8, exp (—A,Z), (18) 
LO am 


where ®, = @; . At large distances down the duct Nu(z) approaches a limiting minimum 
value. If \, is the smallest of the \’s we have, as z —~ ©, that 40,Nu(z) ~ ,@, exp 
(—X,z) and 6,,(z) ~ B, exp (—A,z) so that Nu(o) = d,/4. For experimental measure- 
ments a mean coefficient is generally more useful than the local coefficient. This is based 
on the total heat, g, transferred to the wall between the thermal inlet and the section 
¢. Definition of this coefficient again depends upon the choice of AT and the one most 


commonly used is the logarithmic mean temperature difference 


i A T max — AT min _ (7, — T,) — (Ty — T;) . 

In (A 7 max) — In (A T' min) In (7, — 7,) — In (Ty — 7);) 
Adopting this definition in the fundamental equation we have gq = h,,SfAT,, . Now q 
can either be obtained by integrating Eq. (15) from zero to ¢ or, alternatively, it is the 
heat given up by the fluid in cooling from T, to Ty , so that gq = AwopC, (T,) — Ty). 
Equating these two and introducing dimensionless quantities we have the mean log- 


arithmic Nusselt number, h;,d/k, given by 


] ] 
Nu'(z) = *s In (+). (19) 
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The advantage of basing the mean Nusselt number on the logarithmic mean temperature 
difference is that Nu’(z) tends to the same limiting value as the local coefficient Nu(z). 
For as z— ©, 04 ~ @, exp (— A, 2) and hence 


Mee 1 1 1)\t. 
Nu'(z) ~ i, + In (2h (20) 


The foregoing results are based on fundamental definitions given by Jakob [4] and are 
true for a duct of any cross-section. 

Basis of the method of solution. We consider the general domain D’. The following 
is similar in principle to the method of Galerkin [5]. Let {¢,} be the complete set of 
eigenfunctions of the equation 

Vid + Ad = 0, (21) 
with 0¢/dv’ = — N¢onC’. Any arbitrary function O(z, y) which satisfies these boundary 
conditions and which possesses continuous partial derivatives up to the second order 
can be expanded in an absolutely and uniformly convergent series in the form 


ox 


O(z, y) = be AnPmn(X, Y), (22) 
where 
1 : [| O¢,, dx di (23) 
Go. = ~* ae m d * 
‘ 6,,(m) JJ p- y 
and 
5,(m) = I ¢,o, dx dy, (24) 
vs D’ 
so that 6,(m) = O unless m = p. We can make @ the solution of Eq. (9) so that multiply- 


ing this equation by ¢,, and integrating over D’ we have, by Eqs. (21) and (23), 


5,,.(m)A,»Qm = Xr [| r(x, y)o,,0 dx dy, (m = 0,1, 2, ---), (25) 
vs DP’ 


where r(x, y) w(x, y)/we . If we substitute for @ by Eq. (22) then Eq. (25) is reduced 
to an infinite set of homogeneous linear algebraic equations 


be {5,(m)A,, — Ab,(m) fa, = 0, (m = 0,1, 2, ---), (26) 
where 
b.(m) = I rond, dx dy. (27) 
vs D’ 


The matrix associated with Eqs. (26) is symmetrical since b,(m) = b,,(p) and the elimi- 
nant for a non-trivial solution gives an infinite determinantal equation A(A) = 0 whose 
latent roots are the eigenvalues of Eq. (9). Dividing each row of A(A) by its leading 
diagonal element, the resulting determinant converges [6] if the off-diagonal sum is 
absolutely convergent and \ has no value which makes a leading diagonal element 
zero. If this condition is satisfied the convergence of >—*., | a, | and, moreover, 
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of >°*., 5,(p) A, | a, | follows. The eigenvectors {a‘”} corresponding to a given root 
\ = X,, ean then be obtained, theoretically, in terms of any arbitrary coefficient but in 


practice the determination of a given eignensolution is a problem in numerical analysis. 


One special point concerning the above formulae may be noted. It will be necessary in 
the rectangular case to specify solutions of Eq. (21) by number pairs, that is {¢,,.,} 


rather than {¢,,}, and the expansion (22) is now a sum over all number pairs from 


m,n = 0, 1, 2, --- . Thus 6,(m) in Eq. (26) is then written b,.,(m, n) and is associated 
with a coefficient a,., in a double sum over number pairs from p, g = 0, 1, 2, --- . The 
equations hold for m, n = 0, 1, 2, , and 6,.,(m, n), written for 6,(m) in Eq. (24), is 
non-zero only if both m = pandn = q. 


The rectangular cross-section. In this case the boundary conditions become 


00 0 00 z 0 
- +NO when r=,, ; = +NO when y=—,,, (28) 

02 l Oy l 
where / = (1 + a@)/2, l’ = (1 + 1/a)/2 and a is the aspect ratio. Now the functions 
X,,(«) = sin (6,2 + 8%), where tan 6% = 6,,/N, (0 < BY < 2/2), satisfy the first of 
Eqs. (28) if 8,,(m = 0, 1, 2, ---), is a positive root (the negative roots only repeat the 


functions) of the equation 
tan Bl = 2NB/(B° — N°). (29) 


The roots of Eq. (29) form two separate sets which satisfy respectively the equations 


Btan 361 = N 
and —s (30) 
N tan 36/1 + 6 = O, 


the corresponding functions being respectively symmetrical and anti-symmetrical 


about x = 1/2. The root 6 = O of the second equation does not contribute a function 


X,,(x) but the root 8 = 0 of the first in the case N = 0 contributes a function X,,(z) = 1, 
which must be included for completeness. A similar set of functions Y,(y) = sin (yay + 
+’), where tan y, = yv,,/N, (0 < yi < 2/2), satisfy the second of Eqs. (28) if y, (n = 
0, 1, 2, ---), is a root of Eq. (29) with /’ for 7. Adopting the double suffix notation 


defined above, ¢, X,,(x) Y,(y) satisfies Eq. (21) with the results 


Aen = Ba + Ya (31) 
and 
2 _ Té 9 7 2 ofits T2 9 T / , 
as ae a su (2 N? +2N t\(% + N TaN iL), (39) 
m7 $ . 72 2 72 OL 
Bn +N %¥ +N 
We can now obtain a formula for b,,,(m, n), given by Eq. (27), in the form 
b, Qin, Nn) = z{Cc m—p|,\n—a| — €\m—pl .n+e + Cm+p.nta Cosateents (33) 
where 
Cuan ave | r cos }(6,, Bx , 
pinse | | + B, (34) 
+ (B87, + Bj)} cos {lm + vay + (vi + vi)} dx dy, 
and a suffix | m — p |, say, involves a change of sign between elements with suffixes 


m and p on the right hand side of Eq. (34). Note that this notation is defined only for 
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compound suffixes. A term ¢,,,,, has no meaning in the general case, although exceptionally 
it has in the following limiting cases. Putting N = ©, the constant wall temperature 
case, we have B/ = 0, vy, = 0,8, = mr/l, y, = nw/l' and bn..(m, nN) Ann = (mr /4a) 
(m® + a’n’). The expansion (22) is now in the form of a double Fourier sine series in 
(0<2<10< y < I’). The coefficients given by Eq. (34) can now be identified with 
members of the set whose general term is 


l zg? 
d..,= | [| rcos(émx/l) cos (jny/l’) dx dy, (6,7 = 0,1,2 +++), (35) 


that is, they may be associated with the coefficients of the double Fourier cosine ex- 
pansion of r(z, y) in (0 < x < 1,0 < y < I’). Equation (33) still holds identically with 
d for c. On the other hand, if N = 0 the only difference is that 6, = y, = w/2 and we 
find that we can write d in place of c in Eq. (33) provided we change the negative signs 
in this equation. This case is of no interest in the present problem but may be so in 
other applications since the above formulae are true for arbitrary r(x, y). In the present 


problem r(z, y) is found from Eq. (3) to be 


r(x, y) = ; 3 rs i*7 (2? + aj?) sin (irx/l) sin (jry/l’), (36) 
where 
fo = (4/n”’) : > i 777°? + aj’) 
and, because r(x, y) is symmetrical about both z = 1/2 and y = l'/2 , 7 and j are 


restricted to be odd integers only. Substitution into Eq. (34) yields a formula for c,,.,..., 
which is expressed as a double summation with respect to 7 and j but which can be 
summed with respect to one of these variables of summation to give a rapidly convergent 
single series. We also find that, because of the symmetry present in r(z, y), Cnepneg 1 
zero under certain circumstances. Let us associate odd integers with the roots of the 
first of Eqs. (30) and even integers with the roots of the second. Then it is readily shown 
that Cnsp.n+¢ iS non-zero only if m + p and n + q are both even integers, and that the 
same applies to each of the other three coefficients in Eq. (33). The formula for the 


nm 


non-zero coefficients is 


> raj tan {4(6,, + 8,)l} — (8, + 8,)l tanh (sajr) (37 
¢. — ae <M Ts a Sarees OF eo a: eee Oo" od) 
' OL jaja + Pb. + By yiare — Um +) } 


where ¢,45..+¢ = 7 l’ cos (8, + B/) cos (y, + ¥/)/afo (8m + 8,), and jis odd. The other 
three coefficients required in Eq. (33) are obtained from Eq. (37) by appropriate changes 
of sign. Since, for all values of N, the roots of Eqs. (30) all approach values which are 
integral multiples of z/l, it is clear that the quadruple sum of b, ,(m, n) with respect 
to p, g, m, and n converges absolutely so that A(A) converges. 

We may now consider the special nature of the eigenfunctions derived from the 
solutions of the algebraic equations. Since b,,,(m, n), which is the coefficient of a,., in 
the (m, n)th equation, is non-zero only if m + p and n + q are both even it follows at 
once that the equations break up into four independent sets. Since also a,,, is the co- 
efficient of X,(x)Y,(y) in the expansion for © there are four corresponding independent 
sets of eigenfunctions which exhibit the alternative properties (7) symmetry about both 
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x = 1 2and y = I’/2, (ii) anti-symmetry about x 1/2 and y = I'/2, (ili) symmetry 
about x = //2 with anti-symmetry about y = /’/2 and (iv) the opposite of the last 
ease. Of these only (i) concerns us in this problem since by Iq. (13) only these solu- 
tions give non-zero coefficients in Eq. (12). The special case a = 1| needs further con- 
sideration since here, in addition to the matrix symmetry b, ,,(m, n) b,. »(p, 7) present 
in all cases, we also have b,.,.(m, n) = b,,,(n, m). It follows that if an ordered set of 
coefficients {a,.,} with a given eigenvalue \ = X, satisfies the algebraic equations then 


so, with the same eigenvalue, does the set {a,,,} obtained by interchange of a,,, with 
a,.¢ - that is, if Of(7, y) = 2. be 1 Ay gX,(x)Y,(y) is a solution then a linearly 
independent solution is Of(y, 2x) Ye, Doe, a,,X,(x)¥,(y) and, since the eigen- 
values are equal, these solutions may not satisfy the orthogonality property given 
by Eq. (10), which would invalidate Eq. (13). On the other hand the sum and difference 
of these solutions are both themselves solutions and we can write their contribution to 


the right hand side of Eq. (12) as 


a, {O/(x, y) + Of(y, x)! + AZ{O%(a, y) — O*(y, x)}. 
Multiplying Eq. (12) by w{/O/(x, y) — Of(y, x)! and integrating over D’ we find at 


once that @’ 0 and in the remaining term, considered as a single eigenfunction with 
eigenvalue \,, , the terms X,(x)Y,(y), X,(x) Y,(y) oceur with equal coefficients. It follows 
that in the case a = | we can ab initio puta,., = a,,,(p,q = 1, 3, 5, ---), in the algebraic 
equations and that each eigenfunction derived from this reduced set of equations corre- 
sponds to a unique term in the expansion (12) with @, given as usual by Eq. (13). We 
have, of course, assumed that the X, of the reduced set of equations are themselves 
distinct. It follows also in this case that the expansion (8) consists only of functions 
for which @(z, y) = 9O(y, x), which we would expect physically. 

It remains only, in the general case, to evaluate @, from each computed 0, . Sub- 


stituting in Eq. (13) from Eq. (17) we have 


l F ee. , 
SS = = | 7- ds 
A, 


Jc Ov 


and since 
) 


F fs ) r° ff. 
AE he es , "Jd On p ‘ 0 9, \ | as 
| —== ds’ = -2 | oe f ms dx + | {a . is ; dy | (38) 


01 LO 
then 
l6aN~ : ae , oe ' 
Qs 3 z } sin 8B, sin y((6,, + Yq )@n.a + (39) 
(1 . Qa) A, p=1 q=1 
A more rapidly convergent formula is found by substituting directly for 0, into Eq. 
(13) but it is more complicated except in the special case N = ©, in which it becomes 
I = - - 4 2 .2\—1 n 
ws ale tf, zs » Pp ‘Gg '(p +a’) 'a,, . (40) 


In these formulae f{a”’,.,} are the particular set of coefficients which refer to 9, and 
which satisfy Eq. (11). In practice, solutions of the algebraic equations have been com- 
is such a solution 


n ’ 
pas 


puted by arbitrarily putting one coefficient equal to unity. If {A 
and we put A,") = %, a,",, then %, is found from Eq. (11). From Eq. (9) we 
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have | fo w2? dx dy = (Wo /rn) yf So », Ve, dx dy so that 


> ¥ 5... g(Be + y[AL?P. (41) 
ae ae oe ee ted me 
Computational results. For a given value of a the computation from the algebraic 
equations of the first few eigensolutions, arranged in ascending order of \, is a standard 
problem. We have used relaxation methods in conjunction with Rayleigh’s principle, 
\, for a given eigensolution being estimated by the Rayleigh quotient 


I 


. r lL ip 
+ ~ | 0,70, dx dy * | rO, dx dy, 
z a Jo /0 ( $2) 


> d 6,..(p, 0+ ¥ p> p> > x i. £20 


Since the diagonal elements strongly dominate the matrix, good initial approximations 
to the eigenvalues are found by equating each successive diagonal to zero giving estimates 
i 5,,.,(m, n)(BZ + ¥2)/bm..(m, n) where, for obvious reasons, we temporarily adopt 
double-suffix notation for the \’s. We obtain an initial estimate of the (m, n)th eigen- 
function using this approximation for \,,,, , arbitrarily putting A‘”;"’ = 1 and using the 
other equations except the (m, n)th to estimate the remaining coefficients in the form 


eS M(B, + 7, 4) ‘ 
: ~ bh. (p, Fi —~ ba.n(m,n) — b, .(p, . : 
1, (Ps ) E (m, n\(82 + 4) Ym.niM, nN al Ps Q) (43) 


This approximation is computed until the coefficients attain some prescribed order of 
smallness acceptable as error. It is then improved by the relaxation process until a final 
accepted solution is reached upon which we impose the criteria that no coefficient larger 
than the agreed order of smallness be omitted, nor must any omitted coefficient affect, 
to the agreed working accuracy, the value of any included coefficient. The eigenfunctions 
themselves are of little physical interest but to illustrate the numerical — we give 
in Table 1, correct to four decimal places, the first two in the case a = 1,N = o. The 
physical results of interest depend on the expansions (14), (18) and (19). The @’s converge 
fairly rapidly and we need only the first few terms to describe the physical domain, that 





2 > 2 (say). It is pointless to extend the solution too near to z = 0 for since z = £/dPé 
TABLE 1 
az=1,N=o 
q ] 3 5 7 9 
Pp 

l 1.0000 —0.0369 —0.0022 —0.0004 —0.0001 —0.0001 17 
3} -—0.0369 0.0043 0.0003 —0.0002 15 
5} —0.0022 0.0003 —0.0003 13 
7} —0.0004 Or —0.0011 ll 
9} —0.0001 Q2 0.0001 0.0001 —0.0021 9 


0.0001 0.0004 0.0001 —0.0068 
ne 0.0001 0.0004 —0.0002 0.0334 —0. 1583 
0.0001 0.0001 0.0334 —0.2967 1.0000 


ee 
=().0001 —0.0002 —0.0003 —0.0011 —0.0021 —0.0068 —0.1583 1.0000 0.3711 








zr Kewo 





17 15 13 11 9 7 5 3 l q 
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then, for fixed ¢ and a given fluid, this would be equivalent to increasing the Reynolds 
number beyond the laminar flow range. In our results we have treated the constant 
wall temperature case in considerable detail since this is the case dealt with by Clark 
and Kays (loc. cit.). In Table 2 the first three \, and associated @, are given in each of 
the cases a = 1, 2/3, 1/2, 1/4, and 1/8. For the more general radiation boundary con- 
dition we have considered only the square duct. The first three terms for 6,, in the case 
N = 2 are 

64% = 0.972 exp (—4.8lz) + 0.023 exp (—47.6z) + 0.003 exp (—127z)- + 
while the first only for the cases N = 10, 20 are respectively 

64 = 0.893 exp (—9.18z) + --- and 6, = 0.860 exp (—10.37z) + 


The eigenvalues are well separated for the square duct so that in fact even the first 


term describes well the physical domain. Beyond N = 20, X, [and hence the important 
quantity Nu(@)] can be calculated to good accuracy from the formula 
Qi)v_ ne Bi (6 = 2 y+ Ny [bi 1 » 1)]y. : (44) 
(Ay) v= r B; + V* [b, (1, 1)]y 


This formula is based on the assumption that the correct value of \, for a given N bears 
the same ratio to \, in the case N = © as do the corresponding diagonal estimates of 
these eigenvalues in the initial approximation given previously since, by Rayleigh’s 
principle, these latter are always over-estimates. For example when N = 20, Eq. (44) 
gives \, = 10.42 against the correct value \, = 10.37. 

Comparison of thermal results. The theory used by Clark and Kays is based on 
the assumption, previously used by Seban and Shimazaki [7] for turbulent flow in 
cylinders, that far enough from the thermal inlet 


a(T-T, ) 
om | = —1} = ( 
2 (; mre 7 Reale 


and this is borne out in our work since, for large enough z, 0/04 ~ 9,(x, y)/@, 

independent of z. Calculated results for the limiting Nusselt number in the cases a = 1 
and 0.5 are given respectively as Nu’(~) = 2.89 and 3.39 and these compare well with 
No theoretical information on the variation of Nusselt 


, and is 


our values of 2.98 and 3.39. 
number in the thermal inlet length is given, but experimental data have been obtained 


by the authors for the aspect ratios a = 1 and 0.382. These largely confirm the theoretical 





B 0.014 0.043 0.020 0.028 0.034 


TABLE 2 
AN 

| 
1.000 0.667 0.500 0.250 0.125 | 
My 11.91 12.49 13.57 17.76 22.38 | 
71.07 51.58 $1.17 28.17 25.61 
157.9 99.71 94.93 17.82 31.81 | 
Bi 0.804 0.802 0.789 0.756 0.737 | 
BR 0.104 0.064 0.071 0.107 0.091 | 
| 
| 
| 
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Fig. 1. K against aspect ratio for the case of constant wall temperature. 


values of Nu(@) and also it is found that for small values of dPé/¢ the variation of the 
logarithmic mean Nusselt number is linear according to the formula 


Taal, , 
ae r(#Pe"). (45) 
Nu () 

This linear law follows theoretically from Eq. (20) which gives the theoretical value 
K = (1/A,) In (1/@,). In the case a = 0.382 the theoretical value A = 0.016 agrees well 
with the experimental value 0.017, but this is not so for a = 1 where we find K = 0.018 
against the experimental value 0.042. The theoretical curve for K against a is compared 
with Clark and Kays tentative linear correlation in Fig. 1. The disagreement is serious, 
but we must point out that the experimental curve is determined by only two observed 
points, the result for a = 0 being theoretical, so that an error in an observed point could 
give a very different curve. On account of the discrepancy for the square duct we have 
investigated the radiation boundary condition fully in this case and the results for K 
against N~’ are given in Fig. 2. Clearly the values of K are always lower than those in 
the constant wall temperature case, so no possible explanation is forthcoming from 
these results. On the other hand, comparison of experimental and theoretical values of 
K for the circular cross-section [8] suggests that experimental values may be considerably 
higher. This may lessen the discrepancy in the square case but, for consistency, the 
experimental value for a = 0.382 should also be higher. This is possible since Clark and 
Kays state that in this case the ratio of duct-length to mean hydraulic depth used in 
their apparatus could possibly be higher than the assumed value by as much as 100%; 
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Fic. 2. K against N~ for the square duct. 


this would lead to a larger value of K. In the absence of more detailed experimental 
results, however, it is impossible to state precisely the cause of the disagreement. Finally, 
we should notice that there is some doubt regarding the theoretical value for A near 


the limiting case a 0. We consider only N = © but the general case is similar. If we 
keep the side of the duct parallel to the é-axis fixed and let the other become large then 
0°0/dy? — 0, w/w, — 24 x(1/2 — x) in Eq. (4). The solution for 6 may now be written 
A(x, 2) = >-2_, @Jd,,(x) exp (—AJz) where 3” + 24 N’x(1/2 — x) ¥ = Owith 3) = 
3(1/2) = 0. Now each #,,(z) can be written as 
Tore) lw = 1,8, 8 «+, 
T if) 


7 = 
that is, it can be considered as a sum of functions 9,,.,(2, y) with identical eigenvalues 
d/ and these latter functions can, for varying m and n, be identified with the limiting 
solutions, here written in double-suffix notation, of our previous algebraic equations. 


It is therefore clear that when a = 0 the value of @, to be used in Eq. (20) should be 


the sum > Bin», (n 1, 3, 5, ---), of the @’s associated with 0, ,(2, 1) and since 
it is easily shown that @,,, = nm ~ ®,,, this sum is x @,,,/8. When a ¥ O and the X’s 
are all distinct we should use @,., which, in double-suffix notation, corresponds to the 
smallest \. There is therefore some doubt as to the correct procedure near a 0 but, 
in practice, it can make very little difference to Fig. 1 since K is so small at this end of 


the curve. 
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BOOK REVIEWS 


(Continued from p. 284) 


estimated first mode frequencies [by the Raleigh method] will ordinarily be too high”. One could add 
many more examples to this list of incorrect or misleading statements. 

In some respect, the book lacks order. Aside from questions of sequence of subject matter which 
appears to the reviewer often unmotivated, one finds a portion of a discussion of the nonlinear pendulum 
problem, as well as that of cases of nonlinear damping, in Chapter 2 and 4, while nonlinear vibrations 
are discussed in Chapter 7. In Chapter 12, dealing with vibrating beams and related subjects, one 
finds that “mention of damped vibration with many degrees of freedom has been avoided up to this 
point .. .”’ while the dynamic vibration absorber with damping, and the Lanchester damper are treated 
in Chapter 9. 

Perhaps more serious than the errors of commission are those of omission. One of these is the virtual 
absence of references. Instead, each chapter is followed by a list of “suggested reading’’. These lists 
show a serious disregard, both for important work done in the field, and for the relation of the elementary 
level of the book under review to the advanced level of some of the suggested readings. As one example 
for each may serve (1) the omission of McLachlan’s book on “Ordinary non-linear differential equations” 
in the section on nonlinear vibrations, and (2) the inclusion of Mindlin’s paper on the “Influence of 
rotary inertia and shear on flexural motions of isotropic elastic plates’’. This latter is suggested reading 
for students from whom even the elementary equation of motion of the one-dimensional vibrating beam 
was withheld in the body of the book, and who have been left completely in the dark on the existence 
of rotary inertia and transverse shear effects in beams. 

No mention is made of Lagrange’s equations of motion or of generalized coordinates: in fact, there 
is no treatment of theoretical mechanics. The section on nonlinear vibrations contains no phase-plane 
considerations, no perturbation method, no iteration method and, not a single differential equation. 
Duffing’s name or equation are not mentioned and, although there is a paragraph on relaxation oscil- 
lations, Van der Pol’s name or equation are also overlooked. The section on vibrating beams, plates 


ranes (in that order) fails to derive or even to state the equations of motion of any of these 





and mem 
systems. In the treatment of plates, Rayleigh’s method is applied to an approximate energy expression 
without stated restrictions on plate thickness or vibrational amplitudes. In a discussion of the limitations 
of beam theory as given in the book, reference is made to “massive extensions . . . attached to the beam’’, 
and to gyroscopic action, but none to rotary inertia or transverse shear. A chapter on control systems 
contains neither operational methods nor transfer functions, and the suggested reading includes “W einer, 
Norbert, Cyberneties’’. A final chapter on fatigue presents some largely descriptive paragraphs; phenom- 


enological or solid state physics theories and statistical approaches are ignored. 


R. M. RosenswerG 


Handbool of SUPeErSONIC aE rodynamics. Volumes hgcey aye De Supt. of Documents, U. S. 

Government Printing Office, Washington 25, D. C 

The purpose of this Handbook is to provide a compilation of formulae, tables and other related 
information which would be of use to ‘designers of supersonic vehicles.’ The volumes are issued in loose- 
leaf form and all are not as yet complete, provision having been made for the insertion of additional 
material as it becomes available. When completed, this series will consist of six volumes cont ining 
21 sections. Some idea of their content may be gathered from the following list of titles of sections which 
have thus far appeared: 

Vol. 1: Section 1—Symbols and nomenclature; Section 2—Fundamental equations and formulae; 
Section 3—General atmospheric data; Section 4—The mechanics and thermodynamics of steady one- 
dimensional gas flow. Vol. 2: Section 5—Compressible flow tables and graphs. Vol. 3: Section 6—Two- 
dimensional airfoils. Vol. 4: Section 12—Aeroelastic phenomena. Vol. 5: Section 15—Properties of gases. 


W. H. Rei 


(Continued on p. 321) 





THE EXACT SOLUTION OF BORDA’S MOUTHPIECE IN 
TWO DIMENSIONS* 


BY 
CHARLES A. HACHEMEISTER 
(Pelytechnic Institute of Brooklyn) 


The solution of the semi-infinite two-dimensional Borda mouthpiece is well known. 
Its coefficient of contraction o is 0.5. The solution of the finite mouthpiece seems to 
have been avoided, apparently because it involves elliptic integrals. 

The finite Borda mouthpiece, Fig. 1, is formed by two walls DC and D’C of length b 


A’ 


/ Ie 


Pat —_ tt 
D c f 
B 
Freuvre / 


projecting into a semi-infinite reservoir which is bounded by the semi-infinite walls 
A'B’'C’ and ABC having a gap of width 2a between them. Inviscid incompressible fluid 
flows out of the reservoir through the mouthpiece to form a jet which is bounded by the 
free stream lines D’E’ and DE. The jet contracts to the width 2ca at E’E far from the 
mouthpiece where the speed of the fluid is uniform and of value U’. The total efflux from 
the reservoir is therefore 2caU. The speed of the fluid along the free stream lines D’E’ 
and DE is U’. Along the wet side of the mouthpiece and reservoir walls, A’B’C’D’ and 
ABCD, the surface speed varies. It is zero at the corners C’ and C and at the infinite 
points A’ and A. At two places, B’ and B, the surface speed has a maximum. This follows 
because B is between A and C where the surface speed is zero. 

The values of the coefficient of contraction o and the maximum surface speed with 








its location are of interest. 

Because the flow is in the general category of ‘‘potential flow’, the techniques employ- 
ing conformal mapping are applicable. This principle of analysis is not new. 

Starting with the lower half of the symmetrical geometry of Fig. 1 redrawn on the 
z-plane of Fig. 2a (the mapping figures are grouped together on a later page) define: 
the complex potential P = ¢ + yi, where ¢ is the potential function, and y is the stream 
function, so that the z-plane velocity V, = V, + V,i = — (dP/dz)* where the * means 


conjugate. 





*Received August 12, 1958. 
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A pair of mapping sequences, one set involving the potential P, the other the potential 
derivative dP /dz, is formed to coalesce finally, establishing the desired formulas. This 
procedure is indicated by the diagram, Figs. 2a through 2e, indicating the sequential 
transformations. 

The potential sequence is started by assigning to the potential function ¢ and to the 
stream function y values which are consistent with the existing conditions at specific 
places in the z-plane: 

Along the left infinite semi-circle ¢ © (including points A and G), at the point 
D = 0, at points £, F ¢ = — o, along the stream line ABCDE y = caU, along the 
stream line FG y 0. 

These values are then plotted to form the P-plane map of Fig. 2b. The P-plane strip 
is then mapped onto the upper half of the w-plane of Fig. 2c by means of the Schwarz- 


Christoffel transformation 
dP/dw = cal’ /rw: P (cal’/r) In (w/k’). 1) 


The choice of | and — k” as the w-plane locations of C and D respectively is made 


in the light of subsequent developments. 





§y DIAGRAM SHOWING THE SEQUENCE OF 
MAPPINGS USED TO ESTABLISH THE 
DES/RED FORMULAS. 
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The potential derivative (velocity) sequence of mappings is started by utilizing 
certain known facts about the velocity in the z-plane. Rather than mapping the potential 
derivative, the reciprocal of the conjugate velocity normalized with respect to the 
terminal velocity U of the jet is mapped. This function maps into a simple geometry. 
Because 
re —dP/dz = q exp (— %) 

where g is the speed and @ is the direction of the velocity, 
Sys = (C7 q) exp (02). 

\long the rigid walls and the stream line of symmetry the direction, @, is dictated 
U/Vve 


by these boundaries. Along the free stream line the speed q is constant. The values of 
V* along the boundaries can then be tabulated 


\long ABC CD DE FG 

U/Vz U’/q) exp(i/2) (U/q) exp(1) exp( 92) U/q 
vith > U/q > U/ae > U/q>1 r>O>0 o> U/¢>1 

These v: 


ilues are then plotted to form the U'/V*-plane map of Fig. 2d 
The upper half of the ( 


) 


onto the semi-infinite strip in the Q-plane of Fig. 2e 


V*-plane excluding the unit semi-circle is then mapped 


by means of the transformation 
C= 


f In (U/V*). (2) 
Finally the two mapping sequences are coalesced by mapping the Q-plane strip onto 
the upper half of the w-plane of Fig. 2c applying the Schwarz-Christoffel transformation 

dQ/dw = 


. } fe 4 
N(w + f)/[w(w + k”)*(w + 1)). 


N=} fs 


30th V and f are evaluated by integrating between G to A and C to C in the Q-plane 
and along the corresponding semi-circles in the w-plane. The results are: 
I+k with k 


=1—k". 


Putting these into the expression for dQ/dw, integrating and making appropriate 
idjustment for correspondence of points, there results the Q to w-plane transformation 
[wi + (w +k?) W[(w +k?! + kw?) 
Q = In - 3 . (3) 
k*(w + 1)? 

The connection between the w- and 2-plane is effected by substituting the potential 
derivative for the velocity in Eq. (2) and using the potential derivative from Eq. (1) 


l , dz . dz dw rw dz 
Q In= > = In{-l = In|{—-l =} = In | —————]- (4) 
y* dP dw dP oa dw 
Iquating the arguments of the logarithms in Eqs. (3) and (4) establishes, after 
rearrangement, the z to w-plane transformation in derivative form 
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2g = Wt wt HMw + b+ bw! - 
oa k’*w(w + 1)? 


Equation (5) can be integrated to form the z to w-plane transformation. The result 
involves elliptic integrals of modulus k or k’ with arguments that are complex, real or 
imaginary depending on the locations of the points in the z and w-plane. It is much 
easier to integrate Eq. (5) between limits which correspond to specific points, arranging 
the integrand in each case so that the resulting elliptic integrals have real arguments. 
Three such integrations suffice to establish the relationships among the parameter /, 
the mouthpiece length (b/a), the coefficient of contraction a, and the location of the 
maximum surface speed (g/a). 

The relation between the mouthpiece length (b/a) and the parameter k is obtained 
by integrating Eq. (5) between limits corresponding to the locations of points C and D 


in the z- and w-planes 


vr f 1+k sf" w' dw 
—- dz = 3 — 3 Ct ; 
oa ] ; I [((w + k°)(w + ppt“ + #) 
i dw eRe f dw aii dw 
‘| lwlw + k’°)(w + 1)}} ¥ k”? | . (w+) + | , ww t+ 1)} 


After performing the indicated integrations the results are arranged to 


rh ke _ E(k) — k”K(k) a i +s. (6) 
Q2ca 1—k l—k ~ £8 | 


Integrating between limits corresponding to the locations of points D and E£ establishes 
the relation between the coefficient of contraction « and the parameter k. Here the 
infinite value of x is avoided by approaching point £ in the z-plane: 


; of" Ii+k (w+ 1)! dw kd + k) 
——_ dz 79> =: 72 , = jo 
lim om | g RE | wa [ww +k } 
; dw Ii+k er dw dw 
e | : 2 art 72 - J — #2 | spar j -+- | : ne F 
: [wiw + k~)(w + 1)] k Jpn (w+ 1) : » wiw + )d) 


The real parts are obviously infinite and so are ignored. The imaginary parts yield 


the desired formula: 


© E(k’) — K(k’), x : 
, ere K? Ve 2 a 


Eliminating o from Eqs. (6) and (7) establishes the relationship between (b/a) and 


the parameter kf 


b_ k+ K(k) — k” K(k) — [0 — *)/2] In (1 + 4/0 — k)] (2) 

a E(k’) — k’K(k’) + r(1 — k)/2 
It is interesting to note that k = 1, k’ = 0 corresponds to b/a = ©, the semi-infinite 
Borda mouthpiece. Putting these values in Eq. (7) yields o = 1/2, the well known 
coefficient for this mouthpiece. Similarly k = 0, k’ = 1 corresponds to b/a = 0, the 
mouthpiece vanishes and discharge is through the remaining slit in the reservoir wall. 
Equation (7) now yields ¢ = 1/(1 + 2/m), the known value for discharge through a 


slit. 
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The location of the maximum surface speed along the reservoir wall is found by 
integrating (5) between limits corresponding to the locations of points B and C 





a di —(1+ kj) [ (—w)! dw se 
_ dz i= =, ( - 
oa | r k . e 1+k [(—k _—_ w)(—1 = w)}? + : + ke 
.[ es __ a+hi 
J-usy (—0l—k* — o(—-1 — oY 
[ — [ © dw 
d 14k) (— 1 — w)? ie 1+k w(—1 — w)? 
and 
mg — @) _ (1 + k) {ki K(k )— K(k )] se je + 1+ ke (9) 
2ca i 2(1 — k’) 


Iquation (9) with (7) determines the value of g/a in terms of the parameter / and 
so the relative location of the maximum surface speed. 

The value of the maximum surface speed is found using the potential derivative in 
terms of the transformation derivatives of Eqs. (1) and (5) evaluated forw = — (1 +h), 
the w-plane location of point B: 


ve dP _ _dP dw 
' dz dw dz 
Substituting from Eqs. (1) and (5) and making w = — (1 + hk) 





k 


Fig. 3. Graphical representation of Equations (7), (8), (9), and (10) showing the relations of o, the 

coefficient of contraction; b/a, the mouthpiece length; g/a, the location of the maximum surface speed; 

and Qmax/U, the maximum surface speed, to the parameter k. 2a is the width of the mouthpiece. U is 
the terminal speed of the jet. The values on any ordinate correspond. See Figures 1 and La. 
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1—k 


Vi, (10 
eS | 


V,= 
the velocity maximum at point B. 
Eqs. (6), (7), (9) and (10) constitute the desired relationships. The results are shown 
graphically in Fig. 3. 
For computation purposes it is convenient to introduce into Eqs. (6), (7) and (9) 
E(k) — k’°K(k) 


Bs 73 


which is tabulated in Tables of functions by Jahnke and Emde, Dover Publications, 1945. 
The following formulas from Handbook of elliptic integrals for Engineers and Physicists, 
Byrd and Friedman, Springer, 1945, were used for evaluating the various elliptic integrals: 


Equation Formula number 


(6) 233.00, 233.01, 110.06, 1110.07 
(7) 236.00, 236.01, 110.06, 110.07 
(9) 232.06, 321.02, 232.00, 111.03, 122.10 
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ON A FREE BOUNDARY VALUE PROBLEM FOR THE HEAT EQUATION 


BY 
WALTER T. KYNER* 


University of Southern California 


1. Introduction. W. L. Miranker [1] recently published an existence theorem for a 
free boundary value problem for the heat equation. Using a method due to I. Kolodner 
[2], he obtained a functional equation for the free boundary function R(t) and showed 
that the existence of a solution to the functional equation implied the existence of the 
solution to the free boundary problem. He then solved the functional equation by an 
iterative method. 

The mathematical problem which Miranker solved represents the heating of a long 
insulated metal rod which has begun to melt at one end (x = 0) and after a layer of 
liquid metal A units thick has formed, heat is applied at « = 0. The layer of liquid metal 
is assumed to have an initial temperature distribution f(x). It is essential for Miranker’s 
formalism that A be positive and that df(A)/dz be negative. Physically, this means that 
the front separating the liquid and solid metal must be moving before the mathematical 
model applies. The purpose of this paper is to present a constructive existence and 
uniqueness theorem which is not subject to this restriction. 

The problem is to determine two functions u(z, t) and R(é) satisfying the following: 


= U,, 0<2< Rid, 0 <t, 


Ues 
= — q(t) 
u,(0, t) gt), O<4, (1.1) 
u(R(t), t) = —dR(t)/dt, o<¢ 
R(O) = A, 


u(x, 0) = f(x), O< a <A, 


where g is a positive continuous function, and f is a continuous function such that, for 
some constant b, 


0 < f(x) < D(A — 2), O<2<A. (1.2) 


Miranker required that f be continuously differentiable, non-negative, and that 
df(0)/dx = —4g(0), df(A)/dx < 0, (1.3) 


the latter condition being essential for his proof. Although it was not stated explicitly, 
it follows from his conclusion that u, is continuous on the boundary** that f must vanish 
at x = A. Hence Miranker’s initial value function satisfies (1.2). 

In 1951, G. W. Evans [3] published an existence theorem for this problem with 
A = 0 and g constant. His proof consisted of an iterative argument applied to a heat 


*Received November 3, 1958. The research for this paper was done while the author was a Tempo- 
rary Member of the Institute of Mathematical Sciences, New York University. 


**See the proof of lemma 2 in [1]. 
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balance equation. He proved the existence of a solution for ¢ restricted to the interval 
(0, 1/4]. J. Douglas and T. M. Gallie [4], A. Datzeff [5], G. Sestini [6], A. Friedman [7], 
and the present author [8], have proved existence theorems for similar problems. 

2. The existence theorem: There exists a unique solution to the free boundary 


problem. 
Proof: Following Evans, we derive a heat balance equation by evaluating 


’ aR(t’) 
[ | (u,. — u.) dxdt’ = 0, (2.1) 
70 0 


using the boundary conditions (1.1). We obtain 


aR(t 


t A 
R(t) = A+ [ g(s) ds — | u(x, t) dx +f f(x) dx. (2.2) 


be) 


70 ) 


We use this equation to define a transformation S = F(R) by taking R to be a given 
differentiable monotonic function such that R(0) = A, and taking u to be the solution 


of the reduced problem: 
Un. =U, S< so <= Kip). Oo < g, 


u(0,)=—-g(), O<¢4, a 
u(R(t), ) = 0, o <.. «, 
u(x, 0) = f(z), O<2< A. 


If we can find a differentiable monotonic function which is left invariant by the 
transformation F, then it, together with the corresponding solution to the reduced 
problem, satisfies (1.1). In this paper, we show that the boundary function we seek is 
the limit function of a sequence of iterates, Rp = A, R,., = F(R,). 

The sequence of iterates is well defined, for if R is differentiable, and if S = F(R), 
then S is monotonic and differentiable. In fact, 


0 < dS(t)/dt = —u,(R(0), 0). (2.4) 


The equality follows from (2.2) and (2.1). The inequality from the following argument: 
if u, were positive on x = R(t), u would be negative nearby. But then, by the maximum 
principle, « would attain its negative minimum on x = 0. This cannot happen since w, 
is negative there. 

Our goal is to prove the existence of a solution for an arbitrary time interval. The 
first iterative process which we carry out will converge if the time interval is small. 
Then, taking as the initial function the solution to the reduced problem corresponding 
to the limit boundary function, we carry out another iteration with a small time step, 
etc. We show that a finite number of such processes will give the solution over an arbi- 
trary time interval [0, 7]. Uniqueness of the solution follows from the contracting 
character of the transformation F. If we are content with existence alone, we can use a 
standard fixed point theorem of functional analysis and obtain a (non-constructive) 
proof which does not require subdividing the time interval.* 

We will show that the first iterative process converges to a differentiable function 





*This method was used in [8]. Uniqueness was established by a separate argument. 
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if 0 < ¢ < ¢, . Subsequent iterative processes provide solutions to the integral equations 
t R(t) As 
R(t) = A, + [ g(s) ds — [ u(x, t) dx + [ u,(x) dz, (2.5) 
Jty /0 70 


where A,(= R(t,)) and u,(x) are obtained from the previous process. The function u 
is the solution to the reduced problem 


Uae =U , 0<2< Ri, t+ ¢8< hes 
u,(0, t) = —g(b), /i_<4 < ha» (2.6) 
u(R(t),) =0, tt <t<ta, 
u(x, 0) = u,(x), o<2< A. 

In order to prove convergence, we need the following estimates: 


Lemma 1. If u is the solution to the reduced problem, then there exists a number B, 
independent of R(t), such that for all t in the interval [0, T), 


A < R(t) < Bt, 
0<uz,) <BR —2), O<2< Rd), (2.7) 
—B < u,(R(d), ) < 0. 


Lemma 2. If u and v are solutions to the reduced problem (2.6) with boundary functions 
R and S respectively, then there exists qo > 0, independent of the boundary curves and of 
the subdivision, such that0 < t — t, < qo implies that 


tin 
[ | u(x, t) — v(x, t) | de < (1/2) |R— S| *., ts = 6 + Ge.- (2.8) 


Furthermore, for all ¢ > ¢, , 


| u(x, 4) — o(z, | < B| R-— S|*,0 < xz < min (R(d), S(O). (2.9) 
The derivation of these estimates is in the appendix. 
Let 
jt) = min (F,(0, R,-,(d), (2.10) 
k(t) = max (R,(), R,-.(d), 
then if R,., = F(R,) defines the pth iterative process, and if g < q@, 
a Apr in(t) 
Risi ~ Raila < | | Au,(a, t) | dx + [ | Au,(a, t) | dx 
J0 J Ap 
pact | = s “2 (2.11) 
+ | Au,(x, t) | dr, < 1/2|R, — Rs |, + Bl iM) — Al | Ra — Rae |e 


+ B/2|R, —R,-, |% < [1/2 + 2B’q]|R,. — Ra- |. ; 
where Au, is the difference between the solutions to the reduced problems corresponding 


*|}R —S|. = max| Ri’) — S(U) |, <0 St. 
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to R, and R,_, . We have adopted the convention that the solution to the reduced 
problem is identically zero outside the original domain of definition, i.e., u(x, t) = 0, 


ifz > Rit). 
Clearly, if ¢g < min (1/4B’, q), the sequence converges to a monotonic function 


R(t). Let u(x, t; R) be the solution to the reduced problem corresponding to the limit 
function R. Then if R’ = F(R), 


|R’ —R| =| F(R) — F(R,) | +] Ran — RI, (2.12) 
|R’-—R|, < [1/2 + 2B |R—-—Rails + | Ri -—RI., S<t< aH. 


Since the right side can be made arbitrarily small, we conclude that F is invariant under 
F. We repeat this argument for each subinterval. 

The functions R(t) and u(z, t; R) are the solution to the free boundary problem if 
R is differentiable. To prove that F is differentiable, we write 


R(t+k) 


R(t+k) 
[R(t + k) — RO) + | u(r, ) de = [ g(s) ds 
d # Rit) 


. 
R(t) 


— [ | fu(x, i + k) — u(x, d))] de. 


Using the law of the mean and the fact that u is the solution to the reduced problem, 


we get 
(1/k/) [R(t + k) — R(O] = —u (RO, t) + OK). (2.14) 


This concludes the proof of the theorem. 


APPENDIX 


Proof of Lemma 1. In proving inequality (2.4), we found that u, must be non-positive 
on z = R(t) and that ~ must be non-negative in the interior of the domain. To obtain 


the lower bound on wu, , we pick a constant B so that 
gt) << B, O<t<T, 
(al) 
0 < f(x) < B(A — 2), G<2< A. 
We extend f as an even function and take v to be the solution of the heat equation taking 


on the boundary values 
v(x, 0) = f(x) + Bl zi, —-A<2<A, (a2) 
v(+R(t), t) = BR), 6 <4 < fF. 
R(t) is monotonic, so by the maximum principle, 
0 < o(z, t) < BRO, 
(a3) 
0 < v,(R(b), bt). 
Since v is an even function, v, = 0 on x = O. Hence, if we let 


w(x, t) = v(x, 2) — Br — u(z, 0), (a4) 
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then 
w,=v, -B-—u, <0 on «x =0, o-< 4 < TF. (a5) 
By construction, w = 0 on x = R(#). It follows from the maximum principle that 
0 < w(z, 2d), 0<-2z< R(d, S< <7, 
w,(R(t), t) < 0. 


(a6) 


We conclude that 
—B <v/(R(t),) —B<u(RO, db, GS<ee 7, 
u(x, t) < v(x, t) — Br < B(R() — 2), 0<-2< Rid, e<ci< Tt. 


a7) 


Proof of Lemma 2. It follows from (2.7) that (2.9) is valid on the boundary, x = R(?). 
By the maximum principle, it is valid in the interior of the domain. 
If R ¥ S, let = sup {t| R() = S(@}. Then for any r > 0, 


u(x, t) — v(z, t)| < w(z7,0|R—Sleie,l <t<l +r, 
(a8) 
0O<2z<C=R(t) = S(t), 
where w is the solution to 
WU, = Wes ; e<2< Cc, r <<, 
wz, t') = 0, o<2<t. 


(a9) 
w(C, t) = B, <= ¥, 
w(0, t) = 0, P< ¢ 
Note that 
A<C < BT, os f <€T,. (al0) 
Clearly, 
[ u(x, t) — (2, t) | dx < | R— S |e, [ w(x, t) dx. (all) 


Since 


aC 


| w(x, t) dx = 4C Y (1/(2n + 1)°(1 — exp —{(t — ¢’)(2n + 1)’/4C"}) (al2) 
approaches zero uniformly in C and ¢’ as ¢ approaches ¢’, we can restrict r so that the 
integral (a12) is less than 1/2. 

If A = 0, the estimate (2.8) is not needed for the first iterative process. The lower 
bound on C will then be R(t,). 
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—NOTES— 


ON A SPECIAL BOLZA VARIATIONAL PROBLEM, AND THE MINIMIZATION 
OF SUPERAERODYNAMIC HYPERSONIC NOSE DRAG* 


By H.S. TAN (Therm-Electric Meters Co., Ithaca, N. Y.) 
Instead of a definite integral, let the following expression be given: 
l 
D = Gyo.) + | Fy’, 2) de, (1) 
0 


where G is a known function of yo and y; , yo and y; may or may not be specified. It is 
desired to find the optimum function y(x) that minimizes expression (1). This is a special 


case of the Bolza problem [1]. 
The variation of (1) is easily obtained as follows: 


6D = D(y*) — Dy) 


I 
= Gy byo + Gu, by: + [Py bys + fF, — (Fy) by dz (2) 


Il 


G,, — (F,)o] dye + [Gy + (Fei) dyn + [ [F, — (F,,)’] ay dx. 


To insure vanishing of 6D, it is clear that: 
(i) throughout the interval 0 < x < I, Euler’s variational equation must be satisfied: 


¢g = (F,,)’ — F, = 0 (3) 
(ii) at both ends of the interval, the following end conditions must be met: 
bYo = 0, or Wo = G,, — (F,-)o => 0, (4) 
by, = 0, or w= G,, + (Fy): = 0. 
The condition of minimization is then furnished by 
(iii) the Legendre second variational inequality: 
Fy, > 0, OG<2< i. (5) 


Our generalization evidently lies in the inclusion of G term outside the integral, 
and the corresponding broadened end conditions (4). Indeed, if G is constant, or dis- 
appears from (1), then the end conditions simply reduce to the following conventional 


form for calculus of variations [2, 3, 4], i.e. 
ébyo = 0, or (F,-)o = 0, (6) 
0, or (F,.), = 0. 


II 


bY: 


* Received June 9, 1958. 
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In case 


af al 


Gyo 5 Ya) = | F(y, y’, x) dx + | F(y, y’, x) dz, 


¥ 0— Jl 


(1) reduces to following single integral: 
D = | F(y, y’, x) dx (7) 


which is, on the one hand, an immediate extension of the ordinary calculus of variations 
to include Stieltjes integral at the ends; and on the other hand, a special case of our 
present formulation. 

The present generalization of the variational problem has, in fact, arisen as a result 
of the search for optimum nose curve that minimizes the superaerodynamic hypersonic 
drag of an axially symmetric body. The problem is usually so formulated that the nose 
length / and base radius y, are both given, while the tip radius y, may or may not be 
specified, depending on the situation. By using meridianial coordinates, x axial and y» 
radial, the nose drag has been shown to be given by following integral [5]: 


D = 2npV’ {1 + cyl + yy”) |yy’ dx 


l 
= roy + ei + 2 [ yy’(1 + y”)'”? ax |} , 


0 


(3S) 


where c is a constant determined by the ratios of solid surface and gas temperature, 
and of body and molecular speed. It is easy to see that this drag formula is indeed of 
type (1), with 


G = y¥, + (y/o, F = Qyy’(1 + y”)'” (9) 


and end conditions: 


(10) 
sy, = 0, or » =G,, — (F,)o = 0. 


As already pointed out, an optimum solution must be that of Euler’s variational 
equation, which, by putting (9) into (3), takes the following form: 


y= —y"(1 + y”)/y(2 — y”). (1) 


/ 


Its solution, in parametric form, has been obtained, by simple quadrature and use of 
relations y”’ = y'dy’/dy and y” = dy’/dz, as follows [6]: 


y=e(1l+y")/y” 
(201 ys + ' 12 
ec, J201 + wy”) ( ~ ih di , : a 2 
f= x ee OY + 31n[(1 + yy? + y' |e +e. , 
3 y y ) 
On putting (9) into (5), the second variational inequality then requires: 
PF... =2—y">0, ie. y’ < (2)"” (13) 


Differential equation (11) shows there are two branch curves merging at cuspidal 
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point y’ = (2)'”*; solution (12) shows that the axis y = 0 can not be reached, and y 
grows indefinitely both as y’ approaches zero and infinity. The correct branch of curve 
is thus specified by requirement (13) which, together with (12), imposes the following 
limitations on the admissible values of y’: 


e<s < &, ww” <6. (14) 


Now the tip condition Yody, = 0 becomes: 

(i) For specified yo: 6yo = 0. In this case, in view of the limitations on the admissible 
values of y’, it is evident that with prescribed nose length / and base radius y, , there is a 
lower bound for admissible tip radius y, , below which no solution exists. Within its 
admissible range, it is easy to see that specification of xo , Yo , 2: , y: , Uniquely determines 
the four unknowns C, , C, , yi and y/ , by the four independent equations from parametric 
solution (12). 

(ii) For unspecified yo : 


y G,, — (Fy-)o = Yol[l — y(2 + yh) + yt?) *”?] = 0. (15) 


Although it is easy to see that both y, = 0 and yi = © are solutions of (15), neither of 
them can be fulfilled by our solution (12), with finite 1 and non-vanishing y, . However, 
a plot of y against y/ (Fig. 1) shows that the curve crosses axis Yo) = 0 at yj = 0.7862 
< (2)'* Thus, with unspecified y, , y; can be determined through condition (15). Actual 
construction of the optimum curve then amounts to specifying x) , x, , yj (through 
Wo 0), and y, , from which four unknowns y , yf , C,; , C2 , are determined by four 
independent equations from parametric solution (12). 

It is interesting to note that in this case y, corresponding to yi; = 0.7862 is usually 


greater than the least admissible value of yo which corresponds to yi, = (2)'*. This 
implies that with a given finite nose length 7, reducing the tip radius beyond y, corre- 
sponding to y/ 0.7862 actually has an adverse effect. To see this point, it is best to 














Fic. 1. 
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refer to Fig. 1, which shows clearly that both for yo(yj) greater and smaller than 
Yo (.7862), we will have Yodyo > 0, i.e., 5D > 0. 
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NOTE ON THE SOLUTION OF THE NEUTRON DIFFUSION PROBLEM 
BY AN IMPLICIT NUMERICAL METHOD* 


By GEORGE A. BAKER, JR. (Los Alamos Scientific Laboratory, Los Alamos, New Mexico) 
It is of interest that the implicit, numerical method of Baker and Oliphant [1, 2 
for the solution of time-dependent heat-flow problems in rectangular regions may be 
usefully extended to obtain the solution of the time-independent neutron-diffusion 
equation 

VD(nVe — xe + Sir) = 0 (1) 


for the neutron density ¢g. In the original work of [i], to obtain an accurate solution, 
even asymptotically, it was necessary to take At small. This choice was necessary be- 
cause of the occurrence of a term in a higher order space derivative multiplied by At. 
We will herein make a slight modification to remove this defect. Let us rewrite (1) as 
(for the two-dimensional case) 


ee 2 0” 1 d*o* 
Be + Vet aaa a — + fe + =— 33 
GB Ox oY ot Box oy 
‘ (2) 
o 
+ (22 — uot + ue — 30 4 VO-Ver 
© \D() 1 Dir) Dir 
where 
= ot Op _ ; 
=e — = 0, (3 
g=¢ ry, ) 
and 
uw = Min [x(r)/D(r)]. (4) 


To obtain the solution of (1) via (2) and (3) we must advance the time until the asymp- 
totic solution of (2) is obtained. In [1], we described for the special case, D, a constant, 


how to guess g*, and then use 2) to calculate ¢ - by solving the left-hand side and 
eal A > 
then how to compute a new guess, , by means of 


*Received September 30, 1958. Work performed under the auspices of the United States Atomic 


Energy Commission. 
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(F + ee 1)o* + Peale 
ee. A 
¥ ° F + @ ’ (5) 
where 
—3a~? 
A = 2 At ,B=A—-zB, w = A/B, (6) 
F = —x(r)/[D(r)6], W = —S(r)/[D(r)p]. 


In [1] this process was continued until ¢** converged within a desired accuracy to ¢. 
For the steady state, if we obtain g* by linear extrapolation from the two previous times, 
then the first y** obtained is a sufficiently good approximation for g and one can proceed 
directly to the next time step. 

If we assume that D(r) half-way between mesh points is obtained by averaging the 
reciprocals (for conservation of flux at discontinuities), then the proper difference equa- 
tion representation for @(Ax)*(Ay)* times the right-hand side of (2) would be 


B( Ax) B(Ay)* | Wa, 7) + w/5 ‘g(t, )) — 3 “oli, >] — [FG@, 7) +o — 1]e*(Z, j)} 


oa ff BE 1D = DED Vign( — oe 


Bi) = De 13 lds all ni} 
| Meds Di ry | Gd — eh — 1,9) (7) 


any| PE Lt D = DED Vion, j ) — enti, al} 
atany | Beit 1) + Di, ) [e*(t, 7 + 1) — ¢*(Z, j)] 


= —— Dt] [e*(t, as ¢*(t, J in 1)] 


+ 


DG, j) + DG, j -— 1) 
+ 4y*(7, j) — QZe*(i + 1, 7) + o*(t, 79 + 1) 
+¢e*(t-—1,) + e*%@,j3-DJ) + e@+1,5+ ) 
+e7(i+1,j-ND+e*0@-—1,j7+) + ¢*@ —1,j3-)), 
where 'y is yg at t — Atand “gis g att — 2At. 

In order to obtain the asymptotic solution as quickly as possible, we would like to 
choose At as large as possible. Hence we wish to pick — 8 as small as possible. However, 
as was pointed out in [1], the method loses about [— log, (0.2 87(Ar)*(Ay)*)] binary bits 
due to the cancellation of nearly equal numbers. Therefore, if we carry eight decimal 
places and wish to retain four figure accuracy, we must restrict 


B7(Ax)*(Ay)? 2 5 & 10°. (8) 
As a general rule, our experience has shown that about 5 to 10 iterations are required 


and the choice of a smaller 8 than that needed to assure convergence in this number of 
iterations does not speed up the convergence. We find 


B(Ax)*> and B(Ay)? ~ —(u + 50/N), (9) 


where N is the number of mesh points, to be a good rough guide to the size of 8 to be 


used, subject of course to (8). 
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When D(r) is not constant, 8 is subject to another restriction besides (8). If we 
represent ¢*, ¢.aic. , and g** as ¢ plus an error term, it is easy to show using the method 
of Sec. 3 of [1] that if e* is the error in ¢*, then the error in g**, ¢**, is to first order 
approximately 
VD-Ve* 
t# = — NE. (10) 
BD(F + w) 


From (7) it is evident that the x, and y components of (VD/D) are bounded in magnitude 
by (Az)~' and (Ay)~’ respectively. The term F + w has a minimum value of unity. At 


a discontinuity (Az or Ay) YV «* may well be of the order of «*, hence we see — 8(Ax)’ 
and — 6(Ay)* must be at least unity. We find in practice that a satisfactory choice of 


8 may be obtained from the rule that — @(Ar)* and — 6(Ay)’ are greater than 2. If the 
discontinuities in D are slight then we can reduce the value of 8 correspondingly. It 
seems desirable to provide a convergence factor of about } between e* and e** in (10). 

In the problems (NV S 2000) we have run using this method, most took about 7 to 15 
iterations to reduce the error by a factor of 100. Large problems with sharp discontinuities 
in D can take considerably longer, due to the restriction imposed by (10) on 8. 
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FUNCTIONAL EQUATIONS AND MAXIMUM RANGE* 
By RICHARD BELLMAN (The RAND Corporation) 


1. Introduction. The current interest in rockets and space travel has aroused a 
corresponding interest in the determination of maximum range, minimum time, and so 
on, for various types of trajectories. 

A variety of questions of this type have been treated by means of the theory of 
dynamic programming, see [1, 2, 4]. Here we wish to show how to use functional equa- 
tions to determine the range, the maximum elevation, and similar quantities, as func- 
tions of initial position and velocities. 

2. Vertical motion—I. Consider an object, subject only to the force of gravity and 
the resistance of the air, which is propelled straight up. In order to illustrate the technique 
we shall employ, let us treat the problem of determining the maximum altitude. 

Let the defining equation be 


47 


u’’ = —g — hiv’), (1) 


with the initial conditions u(0) = 0, u’(0) = v. Here v > 0, and h(u’) > 0 for all w’. 
Since the maximum altitude is a function of v, let us introduce the function 


f(v) = the maximum altitude attained starting with initial velocity v. (2) 


*Received October 20, 1958. 
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From the definition of the function it follows that 
fv) =v A+ fv — [g + A@)] A) + 0 (A), (3) 


for A an infinitesimal. Verbally, this states that the maximum altitude is the altitude 

gained over an initial time A, plus the maximum altitude attained starting with a velocity 

v — [g + h(v)]A, the velocity of the object at the end of time A, to within 0(A). 
Expanding both sides and letting A — 0, we see that 


te, en 
f'(e 9 + he) (4) 
Since {(0) = 0, this yields 
* _v, dv, | - 
fv) = / <a (5) 


In the particular case where h(v) = 0, we obtain the standard result v*/2g. 
3. Vertical motion—II. Consider the more general case where motion is through 
an inhomogeneous medium. Let the defining equation be 


u’’ = h(u, wu’), u(0) =c,, u’(0) = @. (1) 


Assume that h(u, u’) < 0 for all u and u’, so that c, = 0 implies no motion. 
The maximum altitude is now a function of both c, and c, . Introduce 


f(c, , ¢) = the maximum altitude attained starting with the initial position (2) 
c, and initial velocity c, . 
Then, as above, 
fle, ,¢2) = C2 A+ fle, + c, A, c. + hc, ,c2) A] + 0 (A), (3) 
which yields in the limit the partial differential equation 
Co + 2 zt + hic, , ¢2) er = (. (4) 


By virtue of our assumptions, f(c, , 0) = 0, for c, > 0. 

4. Computational aspects. One can, of course, use the method of characteristics, 
or standard difference methods, to solve (3.4). Let us present another method which 
reduces the solution to the tabulation of a sequence of functions of one variable. 

In place of (3.4), let us use the discrete approximation of (3.3), 


fle, , 2) = C2 A+ fle, + ce, A,eg + he, , c2) Al. (1) 


Since c, is monotone decreasing, it can be used to play the role of time. Let us write 
c, = N6, where 6 is a positive quantity, and f(c, , c2) = fy(c,). We consider then only 
values of c, which are multiples of 6. To overcome the fact that c. + A(c, , ¢2)A in general 
will not be a multiple of 6, we can either replace it by [(c. + A(e, , ¢2)A)/6], or use inter- 
polation. Although use of an interpolation formula slows up the computation, it greatly 
improves the accuracy. For an application of the foregoing techniques to a more compli- 
cated partial differential equation, see [3]. 
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5. Maximum altitude. Consider now the case where motion takes place in a plane. 


Let the equations be 
zc” = g(z’, y'), z(0) = 0, z’(0) = ¢,, (1) 


y'’ = h(z’, y’), y(0) = O, y'(0) = c . 


Introducing, as before, the function f(c, , c.) equal to the maximum altitude, we see that 


fle, 5 2) = (| +.3)'? A+ fle: + gli , cz) A, c. + Ale, , 2) A] + 0 (A). (2) 
Hence, 
eae a 4] 
(ci + ¢5)'”? + gle, , &) : f + hie, , ¢) of =. (). (3) 
‘ Oc, OC 


Once again, let us assume that c. = 0 implies no vertical motion. Then f(c, , 0) = 0 for 
c, > 0. It follows that we can again compute the solution by means of a sequence of 
functions of one variable. 

6. Maximum range. To tackle the problem of maximum range directly requires 
the introduction of another state variable, the initial altitude. It can also be broken up 
into two problems, corresponding to the ascent to maximum altitude, and the descent. 
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ON THE DETERMINATION OF CERTAIN 
THERMODYNAMIC AND PHYSICAL QUANTITIES* 


By A. GLEYZAL (U. S. Naval Ordnance Laboratory, White Oak, Silver Spring, Maryland) 


We consider any physical phenomenon where a quantity z is a continuous differ- 
entiable function of two independent quantities x and y. Thus: 


z= 27, y). 
Hence 
dz =F dx+Gdy, 
where 
. Vz 
FY = F(z, 4) = : 
Ou 
Az 
G = Gz, y) = —: 
oy 





*Received October 27, 1958. 
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Suppose furthermore that the quantities G and y are readily measured directly but F 
and x cannot be measured directly except that the family of curves F = const. may be 
determined and at least two curves x = const. may be found. Then the physical quantities 
F and x themselves may be determined as functions of G and y. We need merely assign 
a label F to each curve of the family of curves F = const. in such a manner that the 
areas enclosed by the two curves x = const. and two curves F = const. are proportional 
to the increment AF of F. Then additional curves x = const. may be constructed and x 
evaluated so that the area enclosed by the two curves x = x, , X = 2, + Ar, and two 
curves F = F, and F = F, + AF is equal to AF Ax. 
Proof. Given any exact differential 

dz= F dx+ Gdy, (1) 

where F and G are functions of x and y, then, if the curves F = const. and x = const. 


are plotted in the G, y plane as shown in Fig. 1, 
the area AA enclosed by the curves: 


F=F,, F=F,+4F, 


oe x= 2, + Az, 











is equal to AF Ax. lor, integrating 


[ dz = [ F dx + [ G dy = 0, 


where C is the eycle which proceeds along curves x = const. or F = const. from 2, , 


F, toz,, F; + AF, toz, + Az, F, + AF, toz, + Ar, F; tox, , F, , we find: 


« 


‘ 
0= AF Ar+ | Gay. 
JC 


Therefore: 
| G dy = AA = —AF Az. 
JC 
These statements may be generalized to three or more variables and are themselves 
generalizations of, and render obvious, certain relationships among thermodynamic 


variables. Here, in the usual notation, since 


dE = T dS — pav, (2) 
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we have 
-| pdV = AA = —AS AT, (3) 
JC 
where C is a Carnot cycle which proceeds first along an adiabatic, then along an iso- 


thermal, then an adiabatic, and then along an isothermal to the starting point. We con- 
clude that the family of isothermals 7 = 7, + nAT and the family of adiabatics S = 
S, + nAS,n = 0, + 1, +2, +3, --- drawn in the p, V plane map out equal areas AA 
in this plane. Thus, if for any gas the isothermals p = p,(V) and the adiabatics p = p,(V) 
are determined by experiment, the ‘‘labels’’ S and 7 for the curves may then be de- 
termined for it is merely necessary to label as 0 the isothermal along which water freezes 
and as 100 the isothermal along which water boils at atmospheric pressure. The labels 
of the curves 7 = const. are then determined by Eq. (3) uniquely. Any pair of curves 
may be labeled S = 1 and S = 2 along any isothermal. The unit of entropy is related 
to the unit of mechanical energy by resorting again to Eq. (3). The labels for the interven- 
ing curves S = const. are also uniquely determined by Eq. (8). Entropy and absolute 
temperature S and 7 as thus determined for one gas must be consistent with S and T 
determined for any other gas due to the principle of conservation of energy. The extent 
to which the areas AA in the p, V plane are equal indicates the extent to which the postu- 


lated equation (2) is valid. 


Correction to my paper 
ON TRANSFER FUNCTIONS AND TRANSIENTS 
Quarterly of Applied Mathematics, X VI, 273-294 (1958) 
By A. H. ZEMANIAN (New York University) 


Expression (38) should be replaced by its (m — n)th positive root and the second 
line below this expression should read, “horizontal line whose ordinate is the (m — n)th 
positive root of (m — n)!.” 


The second conclusion of Theorem 8 should be deleted and its proof adjusted such 
that zeros on the imaginary axis are counted with the zeros in the left half plane (i.e. 
the symbol g should be discarded and the expression n — p — gq should be replaced by 


n — p). 
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(Continued from p. 298) 


Some aspects of the mathematical theory of control processes. By R. FE. Bellman, I. Glicks- 


berg, O. A. Gross. The RAND Corp., Santa Monica, California, 1958. xix + 244 pp. 


The purpose of the book is to provide ‘‘a taste of the mathematical theory of control processes’, 
and this has been accomplished by a blending of techniques and considerable mixing of a variety of 
topics. The authors selected the ingredients from their own contributions while at the same time giving 
references to more complete discussions and to the work of others. The book is far superior to and much 
more interesting than a survey. The authors have purposely not gone deeply into the mathematical 
theory nor have they emphasized applications. The reader interested in either of these aspects can 
consult the references, and the authors so advise him. 

The most interesting feature of the book is the emphasis placed on the formulation of control problems 
and the techniques available for their solution and the illustration that by another formulation of the 
same problem new techniques become applicable. The techniques are those of differential equations, 
variational calculus, linear (Hilbert) spaces, dynamic programming and game theory. 

A fairly general description of the control problem—sufficiently general to encompass problems 
in control of interest to analysts and operations researchers, economists, management consultants, 
and engineers—is given. The state of a physical system is described at time ¢ by an n-dimensional vector 


a(t) = (x,(0), x(t), «++ , 2,(). 


The rate of change of the state of the system is described bv an ordinary differential equation 


dx 
dt 


where the freedom in the choice of the control function f represents our ability to control the system, 
and we wish to improve the performance of the system by a propitious choice of f. The performance 
of the system under control f is measured by a functional J(f), and an optimal choice of the control 
function f is that which minimizes (or maximizes) J(f). 

Although it is impossible here to describe all of the various control problems discussed in this book, 
we can describe the three principal problems. In each of these the differential equation is assumed to 


= G(x(0), f(b), x(0) =e, 


be linear with constant coefficients— 


dx 
dt 


and the control enters as a forcing term. The state x of the system (the output) 


=Ar+f, 20) =c. 


A is a constant matrix 
is then linearly related to the control f (the input) by 


t 


x(t) = y(t) + | K(t — s)f(s) ds, 


where y is the solution of the homogeneous equation (the uncontrolled system). The first control problem 
is one where the functional J and the constraints on the allowable control functions are linear. For 


instance, 
J(f) = | (a,x) dt, 


0s tO Ss mM... Sata fT, t=1,-++ .n, 


und 
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oe 


where 


r 


(y,2) = = Y 2, 


r=! 
is the inner product. 
In the second control problem the functional is nonlinear and there are no constraints. An example 
of a nonlinear (quadratic) functional is 


7 » 7 


J(f) = | (x —c,x—odi+a | (f, f) dt. 
0 J0 


This can be interpreted as a desire to maintain the system in the state c during the time interval [0,7']. 
The first term represents the cost of deviation from the state c, and the second term measures the cost 
of control. Optimal control minimizes the cost J(f). The third control problem is one in which the func- 
tional J is nonlinear and there are linear constraints. 

Part I of the book is a survey of fundamental results on linear functional equations (primarily, 
difference and differential equations). Part II is concerned with control problems of the first and second 
kind. A theoretical solution of a problem of the first kind is obtained using the Neyman-Pearson lemma. 
The techniques involved in computing the optimal control functions are illustrated. Problems of the 
second kind are solved by Hilbert-space techniques. Here again computation of a solution is discussed. 
Part III is a study of three problems of the third kind each attacked differently. First a variant of the 
Neyman-Pearson lemma is used, and next a problem with constraints is solved by a combination of 
classical techniques and ad hoc methods. The third problem is the “bang-bang’’ control problem and 
is solved by methods which utilize techniques of classical differential equations, linear spaces and dynamic 
programming. Their solution of the problem provides information on how to determine optimal switching 
in a bang-bang system, and this is illustrated for a two-dimensional system. In this example we see 
an aspect of the control problem that this book neglects. The practical problem—and this is certainly 
true of a servomechanism—is to determine the control function, not as a function of time, but as a 
function of the state of the system. In this case their methods give such a solution. 

In Part IV there is an introduction to the theory of dynamic programming which gives a general 
formulation of both deterministic and stochastic multistage decision processes. The variational problem 
with constraints, which was previously treated by classical techniques, is now considered to be a multi- 
stage decision process. The discrete version of the continuous process provides a means of studying 
the structure of the solution and gives a numerical method of computation. Part V corresponds to 
Part IV with game theory in place of dynamic programming. The min-max technique of game theory 
is applied to some control problems with nonanalytic functionals J. 

The authors provide, as they say, a taste of some aspects of control theory in which they have 
specialized and do so remarkably well. The book contains some misprints, some errors and a few false 
statements, none of which should cause the intelligent reader great concern. They have covered a wide 
variety of topics in a vast and exciting field of research with great skill. 

J. P. LaASau.e 


Mathematical theory of compressible fluid flow. By Richard von Mises. Completed by 


Hilda Geiringer and G. 8. 8. Ludford. Academic Press, Inc., New York, 1958. xiii + 


514 pp. $15.00. 

Hilda Geiringer (Mrs. R. von Mises) and G. 8. S. Ludford have done the scientific community a 
great service in completing for publication the last work in compressible flow of Professor von Mises. 
Unfortunately the present text is only the first part of what the author had originally intended to be 
a comprehensive work on compressible flow. However, this cannot be considered to detract from what 
has been published, but rather to leave uncovered by Professor von Mises’ unique approach various 
topics of importance in compressible flow. The book is written principally from the applied mathe- 
matical point of view and is divided into five chapters, of which the first three were written by the 
author and the last two were completed according to his plan following his papers and lecture notes. 

In the first chapter the momentum and energy equations are derived including the effect of viscous 
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stresses and heat conduction. Although the concept of a general specifying equation of state is introduced 
most of the special cases refer to a perfect gas. The chapter concludes with a discussion of the propagation 
of small disturbances in an inviscid fluid, along with a delineation of subsonic and supersonic motion. 
The second chapter presents general theorems of fluid motion and the method of characteristics, with 
principal emphasis on the mathematical aspects of problems in two independent variables. 

The third chapter considers one-dimensional flow, with the first section treating the steady case 
including viscosity and heat conduction, while succeeding sections deal with the nonsteady flow of 
an ideal fluid. The chapter concludes on a discussion of shock phenomena with the treatment restricted 
to a perfect gas of constant specific heat ratio. In the fourth chapter the author deals with plane steady 
potential flow and discusses the hodograph method, simple waves, exact solutions, and limit and branch 
lines. In the final chapter the author expands his discussion of hodograph techniques, introduces the 
oblique shock in a perfect gas, and concludes the book with a stimulating discussion on the existence 
of smooth transonic flows. For the reader who is interested, there are some forty pages of notes and 
addenda which supply interesting historical footnotes to the text. 

Certainly Mrs. von Mises and Professor Lundford are to be congratulated on the manner in which 
they completed the text, for the reader’s immediate impression is that the book was written in its en- 
tirety by Professor von Mises. Of course, as in any book as detailed as the present one, it is always 
possible to disagree with certain features of the presentation. Thus, this reviewer felt at times a lack 
of Professor von Mises’ apt physical description of fluid flow phenomena. This is particularly manifest 
in the presentation of characteristics and in the later discussion of simple waves. In addition, this reviewer 
did not always find a consistent level of treatment as evidenced, for example, by the unecessary re- 
striction to a perfect gas in treating shock waves directly after introducing a general state equation. 
Some difficulty may also be encountered as a result of the manner in which some of the material is 
presented, an example of which is the introduction of one-dimensional shock structure before any 
consideration is given to the inviscid Hugoniot shock conditions. 

In spite of these minor criticisms, to the reader who has some acquaintance with the field of com- 
pressible fluid flow and who desires a text which presents the more mathematical aspects of the subject, 
as well as to students and research workers more directly concerned with compressible flow as con- 
sidered from a mathematical point of view, the book is highly recommended. In the years to come 
this book is certain to be recognized for its fund of information. It can be considered a fitting tribute 
to the author. 

RONALD F. PRossTEIN 


Einige nichtlineare Probleme aus der Theorie der selbsttdtigen Regelung. By A. I. Lurje. 
Akademie-Verlag, Berlin, 1957. xi + 167 pp. $3.60. 


This monograph represents a revised and enlarged edition of papers written in the years 1945-1950. 
Its purpose is to bring the modern theory of nonlinear control to the engineers designing such control 
systems. Therefore, after careful exposition of methods and solutions, particular examples are treated 
in each chapter. 

The book consists of four chapters: The canonical form of the equations in the theory of automatic 
control; the stability of control systems with one controlling element; self oscillations in control systems; 
behavior of a control system at the boundary of the stability region. The first chapter serves as a base 
for the following ones, which are independent of each other. 

The state of the controlled system is described by n variables ¢,. In the direct control their deviation 
Ag: from the desired equilibrium values determines the action £ of a controlling element (e.g. motor) 
which should reestablish the equilibrium of the controlled system. In case of feedback a linear combination 
(called c) of é, its derivatives and the Ag; determines the action of the motor. The differential equations 
describing the state of the controlled system are linear; the nonlinearity enters through the characteristic 
of the controlling motor: — = f(c). Various realistic functions f(c) are considered in the examples. 

The canonical form of the differential equations describing the state of the controlled system is 
derived without the use of matrices, but with ample use of determinants and Dirac’s delta function. 

In Chapter II the stability in the large of the control system is investigated with as little restriction 
as possible concerning f(c). A Liapounoff function is constructed for the problem in question with the 
assumption that of(¢) > 0. The conditions for stability of a given system are determined by a system 
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of n quadratie equations, whose solution for n > 2 is tedious and for n > 5 extremely difficult (Lurje’s 
own statement). In addition the special case f(¢) = co + (0) with c = const. > 0 and o&(c) > 0 
or (oc) < 0 is investigated. 

In the third chapter Lurje studies the occurrence and stability of self oscillations in nonlinear 
systems by the method of Kryloff and Bogoliuboff and by the method of Poincaré, Both methods are 
explained in detail, based on a paper by Bulgakoff. Readers principally familiar with the methods 
will enjoy the discussion of their special merits and difficulties in the present problem. 

The fourth and last chapter of the book is based on a paper by Bautin which in turn used a method 
suggested by Liapounoff in 1935. (Unfortunately no translation of this paper of Liapounoff seems to 
be available; his famous book, on ‘‘the general problem of stability’’ (1892) was translated into French; 
some chapters of the latter are expected to be known to the reader of Chapter II). The boundary of 
the stability region is considered to consist of dangerous and not dangerous portions. A disturbance 
starting outside, but close to the not dangerous portion will not deviate too much from a stable motion. 
On the other hand motions with initial values close to, but outside the dangerous boundary of the 
stability domain may deviate considerably from the stable motions. 

The translation of Lurje’s monograph into German is well done and certainly helps getting acqu- 
ainted with these investigations, so important for the design of control systems. 


I. FLiiace-Lorz 


Lectures on ordinary differential equations. By Witold Hurewicz. The Technology Press 
of the Mass. Institute of Technology, Cambridge, John Wiley & Sons, Inc., New 
York, and Chapman «& Hall, Ltd., London, 1958. xvii + 122 pp. $5.00. 





This book is a reprint of lecture notes of a course given by Hurewicz in 1943. It is a very clear and 
readable introduction to the fundamentals of the theory of systems of differential equations. 

The first and second chapters deal with questions of existence and uniqueness of solution, the third 
deals with the properties of linear systems with constant or variable coefficients, the fourth with the 
theory of the singularities of second order systems, saddle points, nodes, foci, etc., leading up to a dis- 
cussion of the periodic solutions of second order nonlinear differential equations. 

For those who are interested in following the more modern developments, there is a list of twelve 
volumes dealing with various aspects of the theory of ordinary differential equations. 


RicHARD BELLMAN 


High-speed data processing. By C. C. Gotlieb and J. N. P. Hume. McGraw-Hill Book 
Co., Ine., New York, Toronto, London, 1958. xi + 338 pp. $9.50. 


The field of high-speed data processing has become a highly specialised branch of the general field 
of computing, and there has been a need for an account of its peculiarities, equipment and techniques, 
although it is still in a state of rapid development. This book goes some way to fill this need, and is 
particularly suited to the newcomer to the computing field who expects to specialise in data processing, 
but those already experienced in the use of general-purpose computers will find many of their questions 
unanswered. 

When trying to deal with a restricted aspect of the computing field which depends so much upon 
new and changing equipment there is a strong temptation to stray into the ramifications of computer 
design and programming principles and niceties, but some discussion of these aspects is unavoidable 
in the absence of books dealing with them specifically. Although such preliminaries occupy a little 
over half the book, an eye is kept on those aspects particularly applicable to data processing, and the 
effect of the special requirements of data processing of equipment compared with those of scientific 
computation are discussed in the earliest chapters, and throughout, actual machines are frequently 
referred to. 

The treatment of representation of information concentrates upon alpha-numeric data and coded 
decimal systems. This considerably assists later discussion on programming and coding. Knowledge of 
the details of binary arithmetic is not essential and is correctly relegated to an appendix. 
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Modern data processing systems are generally centered around a ‘general-purpose’ type of com- 
puting and control unit which usually possesses an instruction code specially adapted to the class of 
work to be carried out; thus there is special emphasis on input, output, store transfers and comparisons 
and so on. The organization within the central unit depends upon well established principles, and the 
discussion of the details of the central unit and various storage systems would seem to be excessive 
compared with the treatment of the various means of input and output upon which data processing 
so largely depends. 

The discussion of instruction code types and address systems is followed by a lengthy section on 
programming and coding, and the authors have adopted a hypothetical instruction code as a model. 
This code is readily understood and remembered, being single-address coded-decimal with mnemonic 
letter code for operation type. It contains special features to assist coding for data processing, including 
block transfers, input and output buffers, a number of operations to assist double length operations 
within the accumulator and storage of discriminations. It does however have the disadvantage, for 
the beginner, of packing two instructions to a word; these are addressable only in pairs, and this compli- 
cates the organisation of control transfers and requires the use of redundant ‘skip’ instructions. It has 
to be admitted that some machines are just like this, and in the words of the authors; ‘This awkward 
feature of the Hypothetical Machine has been retained to allow the programmer his occupational 
prerogative of complaining about the instruction code’. 

The main techniques of programming are illustrated via a number of simple problems of data 
processing type and occupies some 60 pages. The use of autocodes of various kinds is becoming common, 
and this method of coding is described only in the final chapter. 

A considerable simplification in the sections on coding would have been achieved, without lack 
of reality, if this tendency had been accepted, and coding been described in terms of even a simple 
autocode. The use of symbolic addresses alone would have been of help, to the beginner. The treatment 
of machine organization could have been shortened and the more interesting features of data reduction 
dealt with in the later sections dealing with the use of files, sorting, selecting and practical applications, 
could have received closer treatment, as for instance, on discussion of the various ways of marking 
file and record endings, the use of additional tracks, and forward and backward reading. 

Examples of applications have been taken from the commercial fields including insurance, accounting, 
planning and scheduling, and there is some discussion of the use of processors as simulators. This latter 
use would seem to have very great potential for the future development of data processing and would 
have justified expansion. The authors have in fact kept their feet well on the ground in discussing data 
processing as it now exists, and there is little reference to its future. Among aspects of the future of 
data processing which one would like to see discussed are the implications of the treatment of the results 
of continuously recording scientific experiments, an expanding field of data processing which calls for 
high processing speeds; the possible effects of increased processing speeds on code systems; the likely 
effects of new electronic techniques; and so on. 

The extension of data processing to wider fields is still largely dependent upon the elimination of 
the human being at the early stages, and the presentation of large amounts of treated information in 
readily assimilable forms. Future developments would seem to depend upon the design of special input 
and output devices for translating original data, as occur in banking and national census and market 
research, into directly usable media without human intervention; and high speed multi-curve plotters 
for ready assimilation by the reader. 

The book makes a good introduction to the current problems and techniques of data processing, 
it is excellently clear in style, and contains, not least in importance, an extensive bibliography. 


T. PEARCEY 


An introduction to combinatorial analysis. By John Riordan. John Wiley & Sons, Inc., 
New York, and Chapman « Hall, Ltd., London, 1958. x 244 pp. $8.50 


The author interprets combinatorial as ‘anything enumerative” and presents a fairly systematic 
survey of the subject with particular emphasis on the developments of the last fifty years and on the 
use of generating functions. The book should be of interest not only to specialists in the field but should 
serve as a useful reference as well since, between the text and the problems (of which there are about 
200), the book contains the solutions to a considerable number of problems. In addition to the classic 
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material on permutations, combinations, etc., there is a chapter entitled ‘Partitions, Compositions, 
Trees, and Networks’ plus two on “Permutations with Restricted Position’? much of the material 
in which is either new or of recent origin. 

Although the reviewer could find no specific criticism to make of the book, he did not find it very 
inspiring. This may be more a consequence of the reviewer’s attitude toward the subject than the authors 
presentation although the authors style of writing does not seem to convey much enthusiasm 


G. F. NEweELL 


Electronic digital computers. By Franz L. Alt. Academic Press Inc., New York and 
London, 1958. x 336 pp. $10.00 


The book is divided into five parts: 1. Introduction; 2. Automatic digital computers; 3. Coding 
and programming; 4. Problem analysis; 5. Matching problems and machines. Parts 1 to 3 are standard 
material but 4 and 5, which make up well over half of the book, present matter not previously collected 
in book form. In part 4, methods of numerical analysis are classified and discussed in relation to their 
usefulness and efficiency when used with computers; and in part 5 the impact of computers on scientific 
and engineering research is demonstrated by examples drawn from many fields. The exposition is clear 
and concise, and geared to a mathematical level far below that of the other books in this series. 


WALTER F. FREIBERGER 


Ordinary differential equations. By Wilfred Kaplan. Addison-Wesley Publishing Co., 
Inc., Reading, Mass., 1958. xv 534. $8.50. 


« 


This is an excellent introductory text in the field of differential equations. It is carefully written 
and carefully planned, with many important illustrative examples, illuminating graphs, and a large 
number of exercises. 

The first part, about three hundred pages, is devoted to a thorough study of linear differential 
equations with constant coefficients, n-th order equations, and linear systems. For this latter purpose, 
matrices are introduced and applied. A number of engineering applications are given and there is a 
very helpful discussion of the interconnection between various engineering terms and mathematical 
vocabulary. A long chapter, fifty pages, is devoted to the fundamental technique of power series solu- 
tions, and a brief chapter, ten pages, perhaps too brief, to numerical solution via difference techniques. 

The penultimate chapter contains a very readable introduction to the study of periodic solutions 
of nonlinear second order differential equations via phase plane analysis. In the final chapter, we find 
existence, uniqueness, and convergence theorems which the author has wisely postponed to the end 
of the volume. 

The book is heartily recommended for college classes and for those who wish to prepare themselves 
for the more advanced theory of differential equations and its applications. 


RICHARD BELLMAN 


Boundary layer research. Edited by H. Gortler. Springer-Verlag, Berlin, Gottingen, 
and Heidelberg, 1958. xii 411 pp. $16.20. 


This book differs notably from the usual symposium proceedings (so often rather aimless and full 
of ticket papers). A symposium was actually achieved, by felicitous choice of a subject of neither too 
wide, nor too narrow, a scope. And while geography played clearly a role in the selection of the partici- 
pants, the remaining freedom of choice was employed to ensure an average level of contributions well 
above that usual at such international meetings. 

The 31 papers and 22 notes range over the whole field of boundary layer research, with the keenest 
interest discernible, perhaps, in the subjects of stability and separation. 

The volume offers more than a formal record of proceedings, due to the inclusion of a part of the 
impromptu discussions. These make the meeting come partly alive again, and it was reluctance to 
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miss any of them (and especially the pungent and illuminating remarks offered by some of the English 
participants), almost as much as reluctance to miss any of the first-class papers, which made this 
reviewer read through all 411 pages. It must have been fun to be at the meeting, and it is fun to read 
its record. Apart from which, the book will be needed on the shelves of all those with a more than transi- 


tory or platonic interest in boundary layers. 


R. E. MEYER 


Elementary statistical physics. By C. Kittel. John Wiley & Sons, Ine., New York, and 
Chapman & Hall Ltd., London, 1958. ix 228 pp. $8.00. 


This book is based upon a series of lectures given to beginning graduate students in physics and 
is written in a style quite similar to the author's earlier book ‘Introduction to Solid State Physics’. 
Unlike many authors who have tried to write a book that will serve as both a text and a reference but 
succeed in doing neither, this book is only an elementary text. For the more difficult topics such as 
ergodic theory, phase transitions, ete., frequently described in books of similar titles, the author simply 
refers the reader to the more advanced works. The scope, however, is very broad and includes besides 
the usual statistical mechanics, short sections on irreversible thermodynamics, Brownian motion, 
noise theory and other topics not always found in an elementary book. 

Although the author sacrifices rigor at times, both in the mathematics and in the physics, to preserve 
simplicity, this seems to be unavoidable in a book of this type. At no place does the author embark 
on long tedious calculations or deviate very much from a fairly constant level of difficulty. Particularly 
because of this, the reviewer believes that this is the best elementary text book presently available 
on this subject. The book is primarily for physicists, however. 


G. F. NEWELL 


An introduction to multivariate statistical analysis. By T. W. Anderson. John Wiley 
& Sons, Ine., and Chapman & Hall Ltd., London, 1958. xii 374 pp. $12.50. 


This textbook gives a very systematic treatment of multivariate analysis. Throughout most of 
the book the methods of analysis are derived by the maximum likelihood method (for estimation) and 
the likelihood ratio criterion (for testing hypotheses). The first of these two methods is known to have 
certain optimality properties; the second one also has some desirable large sample features and often 
leads to reasonable tests. In this way the author has organized the material successfully and presents 
a unified treatment of the many different subjects that belong to multivariate analysis. 

After introducing the reader to the multivariate distribution the author describes how to estimate 
the mean vector and the covariance matrix and how the empirical correlation coefficients are distributed. 
One chapter is devoted to the 7°-statistic and the Behrens-Fisher problem. 

The next chapter deals with a different type of problem, namely the classification of observations 
into one or two or more normal populations. This topic is discussed from the point of view of statistical 
decision functions in terms of Bayes procedures, admissible classes, etc. 

After a detailed discussion of testing various multivariate hypotheses the author turns to principal 
components and canonical correlations, which are treated in a very lucid way. 

In the last chapter some more advanced problems are sketched briefly, and it is only to be regretted 
that they were not given more space in the book; this is true especially about factor analysis. 

An appendix on matrix theory ends the book. At the end of each chapter a number of problems 
are given. 

This book is an attractive and representative example of the modern point of view on theoretical 
questions among mathematical statisticians in this country. It should be very useful to anybody inter- 
ested in multivariate analysis. 


U.r GRENANDER 
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Elasticity and plasticity. By J. N. Goodier and P. G. Hodge, Jr. John Wiley & Sons, 
Ine., New York, and Chapman «& Hall, Ltd., London, 1958. ix + 152 pp. $6.25. 


This book is the first volume of a series of surveys in Applied Mathematics. In the preface the 
authors make a modest statement, saying that their work would be devoted only to such a group of 
problems in the domain of the theory of elasticity and plasticity which were dealt with in publications 
rather inaccessible or published in non-familiar languages (non-familiar from the English speaking peo- 
ple’s point of view). So it was not the authors’ aim to present an exhaustive and proportioned survey of all 
the branches of the subject in question, but rather to focus attention on those important latest achieve- 
ments which, according to their opinion, are little known to the readers whose first language is English. 
What is known and easily accessible has been omitted or only touched upon. Nevertheless, the book 
plays a much more important role because it represents a kind of monographical approach to some 
fields of the theory of elasticity and plasticity. 

The approach to the subject in Part I (The Mathematical Theory of Elasticity, by J. N. Goodier) 
is different from that in Part II (The Mathematical Theory of Plasticity, by P. G. Hodge, Jr.), both 
in essence and form. However, each of them has its own merits, the first part being more a survey, the 
second a survey with an attempt towards monographical approach. 

Part I (47 pages) deals successively with: two-dimensional problems, the problems of holes and 
fillets (without and with reinforcement), mixed boundary value problems, anisotropic elasticity, thermal 
problems, three-dimensional contact problems, wave propagation, seismical and vibrational problems. 
Bibliography contains 128 entries (65 Western, 63 Soviet). 

Part II (96 pages) deals successively with the foundations of the theory of perfectly plastic bodies, 
of strain-hardening bodies, piecewise linear plasticity, the minimum principles, and a number of appli- 
cations of the theory, such as bending of planes and shells, plane strain and plain stress, beams, rods, 
and miscellaneous problems. A separate chapter is devoted to Soviet, Polish, Hungarian, and Chinese 
papers. Bibliography contains: 345 entries (182 Western, 149 Soviet, 11 Polish, 2 Hungarian, 1 Chinese). 

The closing: part of the book includes Author Index and Subject Index (8 pages). 

The book is interesting and refreshing. Both authors, authorities in their specialities, made a worthy 
contribution, which is very useful and serves the purpose: for experts, the book represents a valuable 
contemporary informative source (with a detailed list of references); for non-specialists, it may be 
considered as an easily understandable introduction to the subject. 

W. Ouszak 


Games and decisions. By R. Duncan Luce and Howard Raiffa. John Wiley & Sons, Inc., 
New York, and Chapman & Hall, Ltd., London, 1957. xix -+ 509 pp. $8.75. 


Mathematical fields, as all other subject areas, are in constant need of redefinition. Applied mathe- 
matics is no exception to this rule. One of the most interesting phenomena of the past couple of decades 
is the wide application of mathematical techniques to areas that previously lacked them. The historian 
of the future may well regard the present era as one of the times when mathematics was most healthily 
challenged and reinvigorated by the posing of new problems, arising from attempts to apply mathe- 
matics to new fields of study. 

Among the new areas of application are the nonphysical (or behavioral) sciences such as anthro- 
pology, biology, economics, genetics, psychology, management science, and sociology. The book under 
review gives a wide-ranging discussion of a number of these new kinds of applications, and may properly 
be regarded as a book on applied mathematics. This is true even though the traditional tools of applied 
mathematics (such as differential and integral equations) nowhere appear in the book. Instead, the 
mathematics of convex sets and other kinds of mathematics formerly in the “‘pure’’ category are used 
to state and criticize the results of game theory. 

Game theory was put on its mathematical feet in the epic paper of John von Neumann in 1928. 
This paper was almost unnoticed for 16 years until, in 1944, the treatise Theory of Games and Economic 

Behavior, by von Neumann and the economist Oskar Morgenstern appeared in 1944. There was a brief 
period of further dormancy of the theory until about 1948 when the first of a flood of papers started 
appearing on the subject. To show how vast the subject has become in this short time, a bibliography 
of articles in game theory up to 1957, recently compiled by the reviewer and his wife, lists more than 





es 
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1000 papers and books on the subject during that period. Some of these articles are expository in nature 
and were written to introduce the subject to a new audience. However, a substantial fraction of these 
papers contain the meat of a new result. 

It is manifestly impossible for a research mathematician, say, to attempt the reading of this many 
papers unless he has a very serious professional interest in the subject. But it is completely impossible 
for a non-mathematician, say an interested behaviorial scientist, to plow through all that material and 
still maintain himself in his own profession. The authors (both of whom received their principal training 
in mathematics) have undertaken the task of surveying the field. Because many of the people who want 
to know these results are not technically trained mathematicians, the authors have merely included 
the statements and critique of the results together with a bibliography of references for those who wish 
to probe deeper into the mathematical background. 

Although both authors are mathematicians, Luce has concerned himself very much with sociology 
and Ra 


we find that they ta 






a has made extensive studies of statistics. Both these biases appear in the book. For instance, 


] 
i 


if 
k of a player’s “security level’’ instead of his mathematical expectation—clearly 
sociological terms. Also the very complete discussion of the statistical decision problem is a reflection of 
Raiffa’s special interests. These emphases are to the good, since they reflect the parts of the theory on 
which each of the authors has concentrated most fully. It might be noted that the discussion of con- 
tinuous games is fairly sketchy. 


\ brief summary of the contents of the book is that it covers the two-person zero-sum game, several 


the current solution theories of the n-person non-zero sum game, individual decision making (sta- 


tistical decisions nd group decision-making (the welfare problem). The content areas to which these 


subjects are most closely related are sociology, statistics, economics, and management science. 

Che authors are to be commended for their courage in undertaking such a formidable task as the 
survey of such a large and recent area of study. The result of their iabors is a marvelous exposition of 
most of the theory. This book is certain to be influential in the future training of students working in the 
theory as well as in stimulating work on specific research problems. It should be of interest to anyone 


who wants to find 


it about new kinds of applied mathematics from either a cultural or technical 
point ol view 

The decision as to whether or not a behavioral scientist or applied mathematician should have this 
book on his bookshelves is easy—he should. 
G. L. THompson 


An introduction to fluid dynamics. By G. Temple. Oxford University Press, New York, 
1958. xi + 195 pp. $4.00. 


Phe preface states that “the object of this book is to provide an introduction to fluid dynamics, 
primarily for students reading for honors in mathematics and theoretical physics.’’ One cannot but 
approve the object and simultaneously acknowledge that the author has achieved it with brilliant 
success. 

Keeping within the bounds of inviscid continuous fluid and with an eye firmly fixed on the physical 
problem, Temple has produced in 190 pages a fascinating account of hydrodynamics. Particular empha- 
sis has been placed on the “fluid body”’ that is a portion of fluid which always consists of the same fluid 
particles. In this connection the argument of 1-3 concerning the application of the laws of motion appears 
to be unconvincing, for it starts from the tacit assumption that the internal forces form a self-equi- 
librating system. 

The reader is led by easy stages through elementary notions to sources, doublets and vortices, to 
distributions of these singularities and the action on a body in a uniform stream. Here, in deriving the 
Kutta Joukowski theorem the author obtains a drag term due to the total source strength in Green’s 
equivalent stratum. To the reviewer this term seems to be necessarily zero, for the normal velocity on 
the surface of a body at rest in a uniform stream vanishes and it is to this normal velocity that the 
stratum is due, The text then goes on to conformal mapping, free streamlines, design of wing profiles, 
axisymmetric flow, and finally slender body theory applied to solids of revolution and checked by exact 
solutions for the ovary ellipsoid. 

The book can be heartily recommended. 

L. M. Mitne-THoomson 
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Mathematics of physics and modern engineering. By I. 8S. Sokolnikoff and R. M. Red- 
heffer. McGraw-Hill Book Company, Inc., New York, Toronto, London, 1958. 
ix + 810 pp. $9.50. 


The present volume shares with its well-known predecessor, Higher Mathematics for Engineers and 
Physicists, the aim of providing a discussion of those topics beyond the calculus which are of importance 
in engineering and physics. Thus, mathematical concepts are first presented in a precise manner and are 
then illustrated in many cases by examples and problems drawn from engineering and physics, and 
throughout the book the authors have attempted to preserve this balance between the formal aspects 
of the subject on the one hand and their application to physical problems on the other. The text is 
divided into nine chapters which deal with ordinary differential equations, infinite series, functions of 
several variables, vectors and matrices, vector field theory, partial differential equations, complex 
variable, probability theory, and numerical analysis, with short appendices on determinants, the 
Laplace transform, and the Riemann and Lebesgue integrals. Since the chapters are largely independ- 
ent, the book can easily be adapted to varying teaching requirements and will also prove useful for 


reference purposes 
W. H. Rew 


Introduction a L’ Algébre Supérieure et au Calcul Numérique Algébrique. By L. Derwidue. 
Masson et Cie, Paris, 1957. 431 pp. $15.80. 


The selection of topics in this text book on advanced algebra is especially happy from the point 
of view of engineers and physicists. The author notes in his preface that, while good texts covering these 
topics are available in English and German, the literature in French is very sparse. The principal sub- 
jects covered are: systems of linear equations and determinants, polynomials, matrices and matrix 
eigenvalue problems and the stability criteria of Routh, Hurwitz and Schur. In all of these topics great 
emphasis is placed on bringing the theoretical results down to a form which permits numerical calcu- 
lation. It is assumed that the reader has access to a desk calculator. Many sample calculations are 


given indicating appropriate layouts for recording intermediate results and pointing out methods for 


carrying along “running checks.” Ill-conditioned systems are mentioned and the advantages of select- 
ing the largest divisor in elimination methods are demonstrated but the general round-off problem is 
not touched. High speed automatic digital computers are not discussed. In the final chapter there is a 
short introduction to abstract algebra: groups, rings, fields, ete. 


All in all, the author has succeeded in his attempt to provide a transition between elementary 


algebra and the specialized fields of numerical analysis and modern algebra. 


S. H. CRANDALL 


Space -charge waves and slow electromagnetic waves. By A. H. W. Beck. Pergamon Press, 
New York, London, Paris, Los Angeles, 1958. xi + 596 pp. $15.00. 


This book is the eighth volume in the International Series of Monographs on Electronics and 
Instrumentation. It covers a great deal of ground and on the whole the material is well presented, with 
more or less complete reference to pertinent published papers. There is a fair amount of mathematics 
used, although this entails not much more than a knowledge of the simpler properties of Bessel functions 
and some matrix algebra. There are some places where the logical structure is rather poor and this is 


complicated further by numerous misprints in some of the more mathematical sections. Also, the 





author occasionally quotes results from published papers in a highly condensed form without complete 
definitions of symbols used. This is again complicated by misprints and it was not always easy for the 
reviewer to follow the deductions without reference to the paper in question. However, these difficulties 
do not detract from the general worth of the book and the reviewer believes that the applied mathe- 
matician who wishes to become thoroughly acquainted with the theory underlying the various micro- 
wave valves in use today may do so by reading this book which is obviously written by one with an 


extensive knowledge of the subject. A condensed summary of the contents follows: 














1959] BOOK REVIEWS 331 


A short general introduction is given in Chapter I in which various types of amplifiers are briefly 
mentioned. Chapter II gives a condensed account of Maxwell’s electromagnetic theory. Propagation 
modes are discussed for rectangular and circular waveguides and expressions for the power flows are 
obtained. Vector and scalar potentials and the Hertzian vectors are introduced. Chapter III deals 
with various slow wave structures and in particular with disk-loaded waveguides, interdigital delay 
lines and both sheath and tape helices. Their properties are discussed in relation to the Brillouin diagram 
(frequency versus phase constant). Chapter IV is devoted to space-charge wave theory and topics in- 
clude the determination of the space-charge reduction factors for cylindrical and annular beams in 
cylindrical tunnels; space-charge waves on Brillouin beams and on confined beams; space-charge waves 
in accelerating fields and in crossed fields; the distribution function approach for multivelocity electron 
beams and plasma oscillations. Chapter V utilizes some of these results in discussing the matching of 
input conditions and space-charge waves on cylindrical, annular, Brillouin and confined beams, and 
on beams in helices. 

Chapter VI is concerned with space-charge waves in klystrons while Chapter VII is devoted to 
travelling-wave tubes and backward-wave oscillators. Some crossed field devices are discussed in Chap- 
ter VIII while some special space-charge wave devices, such as the two-beam tube and the transverse 
current T.W.A. (travelling-wave amplifier), are mentioned in Chapter IX. The final Chapter deals with 
noise phenomena in space-charge wave devices and after discussing shot noise in diodes, the noise factor 
for valves with input resonators and the noise figure of T.W.A.’s, the general theory of noise in beams, 
based on the theorems of power flow in space-charge waves, is presented in the matrix formulation 
adopted by Haus and Robinson, although the proof of the noise invariants is omitted. The smoothing 
of current fluctuations near the potential minimum is also briefly discussed. There are twelve Appen- 
dices, mostly mathematical in content, and these are followed by a few questions on each chapter, a list 
of recent. references, a list of the major symbols used throughout the book and, finally, the index. 


J. A. Morrison 


Selected papers on quantum electrodynamics. Edited by Julian Schwinger. Dover Publi- 
cations, Ine., New York, 1958. xvii + 424 pp. $2.45. 


This book is a collection of thirty-four papers taken from the literature of quantum electrody- 
namics—papers which together with a preface summarize the present state of the theory. In summing 
up the problems still facing quantum electrodynamics the editor points out that “the real significance 
of the work of the past decade lies in the recognition of the ultimate problems facing quantum electro- 
dynamics.”’ 

Papers are included in this collection by the following: Dirac, Fermi, Fock, Poldolsky, Jordan, 
Wigner, Heisenberg, Weisskopf, Block, Nordsieck, Foley, Kusch, Lamb, Retherford, Bethe, Schwinger, 
Oppenhemier, Tomonga, Pauli, Villars, Feynman, Dyson, Karplus, Klein, Kallen, and Kroll. 

Of these papers twenty-nine are in English, three in German, and one each in French and Italian. 


R. TRvELL 


Computability and unsolvability. By Martin Davis. MeGraw-Hill Book Company, Inc., 
New York, Toronto, London, 1958. xxv + 210 pp. $7.50 


This book is an introduction to the theory of computability and noncomputability, usually referred 
to as the theory of recursive functions, a branch of pure mathematics. In the light of recent develop- 
ments in computers, decision problems, i.e., problems “‘which inquire as to the existence of an algorithm 
for deciding the truth or falsity of a class of statements,”’ are potentially of interest to other than pure 
mathematicians. The work is an outgrowth of a graduate course at the University of Illinois and series 
of lectures given at the University of Illinois and the Bell Telephone Laboratories. 

The major part of the book is self-contained and assumes no particular mathematical training on 
the part of the reader. A degree of mathematical maturity, in particular the ability to follow abstract 
proofs, is, however, necessary. Acquaintance with elementary mathematical logic is desirable. 
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The concept of Turing machine is made central to the development. Thus Turing’s approach is 
combined with the methods of Gédel and Kleene to present the various aspects of the theory of com- 
putability in a unified manner. The general theory is applied to combinatorial problems, problems re- 
lated to Hilbert’s Tenth Problem, and systems of symbolic logic. The general theory is extended to 
include the Kleene hierarchy of arithmetical predicates, computable functionals, and the classifications 
of unsolvable decision problems. 

A. A. GRAI 


Modern geometrical optics. By Max Herzberger. Interscience Publishers, Inc., New York, 
London, 1958. xii + 504 pp. $15.00. 


This work is both a treatise and a text on the methods of geometrical optics. As far as the reviewer 
is aware it is the only book of this kind in English (other than the notes of Lumeberg which were never 
published). 

There are essentially two distinct sections of the book. The first part is concerned with the ray 
tracing and the calculation of optical systems; this consists of about the first thirteen chapters. Much 
of the second section of the book is concerned with the general laws of geometrical optics. The charac- 
teristic function methods of Hamilton, the general laws of image formation, the properties of concentric 
systems, and the properties of rotation symmetric systems are naturally an important part of the book. 

The contents of the book are listed as follows: 

Part I, Ray Tracing; Part II, Precalculation of Optical Systems; Part III, General Laws; Part IV, 
Concentric Systems; Part V, Rotation-Symmetric Systems; Part VI, Approximation Theory for Normal 
Systems; Part VII, Third and Fifth Order Image-Error Theory; Part VIII, Interpolation Theory of the 
Optical Image; Part IX, Optics in General Media; Part X, Appendix. 

The only topics that the reviewer felt were noticeably missing, topics that might have been included, 
are those of integral invariants and perhaps something about electron optics; it may be that this is 
asking for too much. 

The book represents a tremendous amount of work. It is well written, and it should be the standard 
work in this field for a long time. 

R. TRvELI 


Principles of quantum electrodynamics. Translated from the German by J. Bernstein, 
with additions and corrections by Walter E. Thirring. Academic Press Inc., New 
York, London, 1958. xv + 234 pp. $8.00. 


This book is a text on the quantum electrodynamic part of field theory. It is more advanced in 
nature than the book of Wentzel entitled “Quantum Theory of Fields.’’ The text is the result of an 
effort by the author to present mainly what cannot be found elsewhere in books on this subject. There 
are four parts of the text divided as follows: 

I. General Introduction: Units and Orders of Magnitude; Classical Electrodynamics; General 
Formalism of Quantum Theory of Fields. II. Free Fields: General Discussion; Special Fields; Matrix 
Elements; Fluctuation Phenomena. III. Fields With External Sources: General Formulae; Emission 
of Light; The Dirac Field in an External Electric Field; The Limitations of Measurability. IV. Interact- 
ing Fields: General Orientation; Scattering Processes; Renormalization Theory; Higher Order Cor- 
rection; Outlook. 

Appendix I, Dirac Matrices. Appendix II, Green’s Functions (Relativistic Wave Equation) 

While this text is intended for the advanced student, there are parts such as section 16 of Part IV 
which will be of interest even to those who only want to know what the problems are at the present 


time. 
R. TRUELL 
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Theorie schallnaher Strémungen, By Kk. G. Guderley. Springer-Verlag, Berlin, Gottingen, 
Heidelberg, 1957. xv + 376 pp. $10.07. 


The study of transonie flow as a special field was stimulated by the failure of classical treatments 
of compressible flow, both by theoretical and experimental means, to describe flow phenomena in the 
Mach number range close to critical conditions. This gap arises because, near Mach number one, the 
linearized theory of compressible flow is invalid and normal High Speed Wind Tunnels become choked. 

Much of the fundamental theory of transonic flow was developed by Dr. Guderley himself and 
it is appropriate that he should now present a connected account of his own and other contributions 
to the subject. Although his book is primarily concerned with transonic small disturbance theory based 
on Tricomi’s equation it is self contained and the theory is carefully introduced with the necessary 
background of fundamental Gas Dynamics. The book will be of considerable interest to Applied Mathe- 
maticians and Theoretical Aerodynamicists. A great. deal of space is given to special solutions of Tricomi’s 
equation and their use in solving the problem of flow past a double wedge profile in a sonic stream. 
This is a field to which the author devoted many years of research, previously published in somewhat 
inaccessible form. It is certainly valuable to have this work in a single volume. The book would, in 
the reviewer's opinion, have been better balanced if more space had been given to other approaches 
to transonic small disturbance theory, notably by Oswatitsch, and to a fuller description of the purely 
numerical solution of transonic flow problems. The book is welcome as the first, and much needed, 
book on transonic flow and is strongly recommended to Fluid Dynamicists. There is still room, however, 
for a more complete account of the theory and also for an account of the difficult experimental work 
in this field. Further, since the book was written a number of transonic flow problems have been solved 
in the U. S. 8. R. by Chushkin, using the successful numerical technique of Dorodnitsin for second 
order partial differential equations of mixed type. This approach certainly competes with the classical 
treatment based on Tricomi’s equation. It is hoped that future editions of Dr. Guderley’s book will 
include an account of this recent work. 

There are 11 chapters in the book. Following a first chapter devoted to the basic principles of 
Steady Gas Dynamics the transonie similarity concept is introduced in Chapter II. The third chapter 
contains a brief account of Linearized Transonic Flow. In Chapter IV exact solutions of the simplified 
transonic potential equation are given for flow through a Laval nozzle and flow of a sonic jet. Chapter 
V gives a fairly complete account of the hodograph method leading to Tricomi’s equation when transonic 
similarity is introduced. Chapter VI is largely a physical discussion of transonic flow phenomena in 
the hodograph plane. The author's special solutions of Tricomi’s equation based on the singularity at 
the sonic point are described in full in Chapter VII. This is followed, naturally, by the solution of problems 
of flow with Mach number 1 in Chapter VIII. Chapter IX presents methods for determining flow fields 
with free stream Mach numbers slightly different from 1, based essentially on a linear perturbation 
of the solution given in Chapter VIII. In Chapter X some special points are discussed such as the re- 
flection of disturbances at a sonic line. The final chapter is concerned with the extension of Guderley’s 
basic theory to axially symmetrical flow. 


M. Hott 


Advances in applied mechanics. Editors: H. L. Dryden and Th. von Karman. Academic 
Press, Inc., New York, 1958. x + 459 pp. $12.00. 


The appearance of a new volume in this series, with its outstanding surveys of selected topics in 
Applied Mechanics, is always most welcome. For it is only by means of review articles of the type pre- 
sented here that one can even attempt to keep abreast of developments in so inclusive a field as that 
of Applied Mechanics, an inclusiveness well demonstrated by the titles of the present articles: ‘Supersonic 
air ejectors” by J. Fabri and R. Siestrunck, “Unsteady airfoil theory’’ by A. I. Van De Vooren, ‘‘The 
theory of distributions’ by Charles Saltzer, “Stress wave propagation in rods and beams”’ by H. N. 
Abramson, H. J. Plass and E. A. Ripperger, “Problems in hydromagnetics’’ by Edward A. Frieman 
and Russell M. Kulsrud, “Mechanics of granular media” by H. Deresiewiez, and ‘Condensation in 
in supersonic and hypersonic wind tunnels” by P. P. Wegener and L. M. Mack. 


W. H. Rerp 
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Notes on analog-digital conversion techniques. Edited by Alfred K. Susskind. The Tech- 
nology Press of Mass. Institute of Technology, Cambridge, and John Wiley & Sons, 
Inc., New York, 1957. x + 877 pp. $10.00. 

The articles, by various authors, in this book had their origin in a special summer program on 
the subject at M. I. T. in 1957. The subject matter is divided into three parts: systems analysis, engine- 
ering analysis of devices and a case study, of which the first bears most interest to mathematicians, 
in particular Chapter 2 on the theory of sampling and the theory of quantizing, and Chapter 3 on codes. 


W. F. FREIBERGER 


Elementary mathematical programming. By Robert W. Metzger. John Wiley & Sons, Inc., 
New York, and Chapman & Hall, Ltd., London, 1958. ix + 246 pp. $5.95. 


This book is addressed to industrial engineers and others with little mathematical background 
who wish to know the mechanics of programming without being burdened with proofs or conceptual 
discussions. It is lucid and provides many practical illustrations but does not discuss the limitations 


of the very simplified models. 
WALTER FREIBERGER 


Notes on Operations Research 1959. Assembled by the Operations Research Center, M.1.T. 
The Technology Press, Massachusetts Institute of Technology, Cambridge, Mass., 
1959. viii + 256 pp. $4.00. 

These are the lecture notes for a special program in operations research that members of the Opera- 
tions Research Center at M.I.T. conducted for persons from European NATO countries at the Center 
for Experimental Aerodynamics in Brussels in August 1959. The scope of the book is indicated by the 
following table of contents: Introduction (P. M. Morse)—Probability (G. P. Wadsworth)—Search 


Markov Processes (P. M. Morse)—Queuing (H. P. Galliher)—Control Processes 


(B. O. Koopman 
Sequential Decision 


(R. A. Howard)—Organization of Operations Research Groups (P. M. Morse) 
Processes (G. E. Kimball and R. A. Howard )—Reliability and Maintenance (G. E. Kimball)—Informa- 
tion theory (B. O. Koopman and G. E. Kimball )—Production Scheduling (H. P. Galliher)—Simulation 
of Random Processes (H. P. Galliher)—Bibliography. 


The measurement of power spectra. By R. B. Blackman and J. W. Tukey. Dover Publica- 
tions, Ine., New York, 1959. x + 190 pp. $1.85. 
This is a republication of “Power Spectra from the Point of View of Communications Engineering’’, 
which originally appeared in Volume 37 of the Bell System Technical Journal. 


Les Principes de la Théorie Electromagnétique et de la Relativite. By Marie-Antoinette 
Tonnelat. Masson et Cie, Editeurs, Paris VI*, 1959. 394 pp. $10.25. 


i 

This book is concerned with the study of the principles of classical and relativistic electromagnetic 
field theory and gravitational field theory. The discussion of the principles of electromagnetic theory 
includes among other things, discussion of the electron theory of Lorentz and the electrodynamics of 
bodies in motion. The fifth chapter begins the discussion of relativity and the remaining three quarters 
of the book does an excellent job of this. The only feature of the book that detracts from its value is 
the rather poor printing job, especially the tensor indices which are sometimes illegible. 


Roun TRUELL 
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Theory of functionals and of integral and integro-differential equations. By Vito Volterra. 
Dover Publications, New York, 1959. xiv + 226 pp. $1.75. 


This is an unabridged republication of the English translation of Volterra’s work published by 
Blackie & Son Ltd. in 1930. A preface by Professor G. C. Evans and a biography and bibliography by 
Sir Edmund Whittaker (originally published as an obituary note) have been added. 


Plastic analysts of structures. By Philip G. Hodge, Jr. McGraw-Hill Book Co., Ine., New 
York, Toronto, London, 1959. xiv + 364 pp. $10.50 


The point of view that the rational analysis and design of structures composed of elastic-ductile 
materials requires inclusion of plastic action in the theory of structures has been attracting more interest 
recently, as evidenced by the increasing number of books published on the subject. The present book 
constitutes a systematic and up-to-date account of the theory of plastic analysis of structures and it 
is broader in coverage than its predecessors. Structures treated in the book include beams and frames 
in bending, beams under combined stresses, plates and shells and slabs with cut-outs. The subject 
matter is approached from the view-point of limit analysis. Emphasized are two basic theorems which 
enable lower and upper bounds on the load carrying capacity of structures to be determined. Part I, 
Bending of Beams and Frames, treats the application of plastic methods to frame-type structures. It 
also deals with such topics as elastic plastic deformations, variable and repeated loading and procedures 
of design for economy in materials. Part II, Structures Under Combined Stresses, is primarily concerned 
with combined stresses in beams, circular plates, and circular cylinderical shells. More general plate 
and shell problems and problems in plane stress have also been covered. Finally the book closes with a 
brief introduction to some of the problems encountered in the dynamic loading of plastic structures. 
References and Problems are given at the end of each chapter. 

The clear presentation and reasonably elementary level of mathematical tools employed will un- 
doubtedly make this text attractive to the interested students (both graduate and undergraduate) and 
also to the practicing engineer. However the structural engineer is probably too much elasticity-minded 
and may need additional convincing before accepting such concepts of plasticity theory as rigid-plastic 
action, plastic hinges and plastic collapse. The critical presentation and discussion of some experimental 
work, in addition to the test results given in the text, may prove helpful in this respect. 

It may also be noted that questions concerning the effects of geometry changes and strain-hardening 
do not receive extensive attention in the book. The question of plastic stability is hardly considered in 
the book and beneficial effects of geometry changes and strain-hardening on the load carrying capacity 
of structures are mentioned only briefly. Some systematic discussion of these topics in the future editions 
of the book would increase the value of this interesting work. 


E. T. Onat 





Notice of 


SYMPOSIUM ON PLASTICITY 


The Second Symposium on Naval Structural Mechanics will be held April 
5, 6, and 7, 1960 at Brown University, Providence, R. I., under the joint spon- 
sorship of the Office of Naval Research, Department of the Navy, and of Brown 
University. The Symposium will be devoted exclusively to the field of plasticity. 
The program will consist of critical surveys in selected areas and of reports on 
original research, with ample time for discussion. The organizing committee con- 
sists of Professors E. H. Lee and P. 8. Symonds (co-chairmen), D. C. Drucker 
and W. Prager. 
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HANDBOOK OF AUTOMATION, 
COMPUTATION, AND CONTROL 


Edited by EUGENE GRABBE, SIMON RAMO, and DEAN E. WOOLDRIDGE, all of 
Thompson Ramo Wooldridge Inc., with sections contributed by a staff of 104 specialists. Written 
and edited with an emphasis on systems engineering, the handbook covers material of direct use 
to all levels of students interested in the associated fields of automatic control and computers. The 
major objective is to provide design data for research, development and design in feedback con- 
trol, computers, data processing, control components, and control systems. The stress throughout 


is on new techniques and components for designing, and developing control systems. 
Volume I—Control Fundamentals. Includes specific coverage of aspects of mathematics as applied 


to control; a compilation of the mathematics of digital computers; the latest techniques and com- 
parisons of different techniques involving computers; and all necessary material on information 
theory and transmission. 1958. 1020 pages. $17.00. 

Volume II—Computers and Data Processing. This volume brings together in one place all avail- 
able techniques for design and use of digital and analog computers, 1959. 1096 pages $17.50. 


iD 


TESTING STATISTICAL HYPOTHESES 


By E. L. LEHMANN, University of California, Berkeley. Gives a systematic account of the mathe- 
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matical theory of hypothesis testing and theory of confidence sets. 1959, 396 pages. S$ 


MEASUREMENTS: Definitions and Theories 
ed by ( W. CHURCHMAN, University of California: and P. RATOOSH, The Ohio State 
" Presents various approaches to the problems of measurement, 1959. 274 pages. $7.95 


METHODS OF CORRELATION AND REGRESSION ANALYSIS, Third Edition 


By M. EZEKIEL. United Nations; and K. A. FOX, Iowa State College. Revised to reflect advances 


1959. 548 pages. $10.95, 


in the field. 


MATHEMATICAL PROGRAMMING AND ELECTRICAL NETWORKS* 


. A new approach to mathematical pro 


By J. B. DENNIS, Massachusetts Institute of Technology 
1959, 186 pages. $4.50. 


gramming based on an analogy with electrical networks 


STUDIES IN DISCRETE DYNAMIC PROGRAMMING* 


B RONALD A. HOWARD, Massachusetts Institute of 


niques in dynamic programming. 1959, In press. 


Technology. Develops specialized tech- 


REGRESSION ANALYSIS 


By E. J. WILLIAMS, University of North Carolina. Deals with regression analysis from theoretical 


and applied viewpoints. 1959. 214 pages. Prob. $7.50. 


CLASSICAL DYNAMICS 


By R. H. ATKINS, Unischol Tutorial College, London. 1959. In Press. 


* A Technology Press Research Monograph Send for examination copies. 


JOHN WILEY & SONS, Inc. 440 Fourth Ave. New York 16, N. Y. 























